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Abstract 

Two  years  of  moored  oceanographic  and  automatic  weather  station  data  which 
span  the  winter  ice  seasons  of  2001-2003  within  Marguerite  Bay  on  the  western  Antarctic 
Peninsula  (wAP)  shelf  were  collected  as  part  of  the  Southern  Ocean  Global  Ocean 
Ecosystems  Dynamics  program.  In  order  to  characterize  the  ice  environment  in  the 
region,  a  novel  methodology  is  developed  for  determining  ice  coverage,  draft  and 
velocity  from  moored  upward- looking  acoustic  Doppler  current  profiler  data.  A  linear 
momentum  balance  shows  the  importance  of  internal  ice  stresses  in  the  observed  motion 
of  the  ice  pack.  Strong  inertial,  not  tidal,  motions  were  observed  in  both  the  sea  ice  and 
upper  ocean. 

Estimates  of  upward  diapycnal  fluxes  of  heat  and  salt  from  the  Upper 
Circumpolar  Deep  Water  to  the  surface  mixed  layer  indicate  almost  no  contribution  from 
double  diffusive  convection.  A  one-dimensional  vertical  mixed  layer  model  adapted  for 
investigation  of  mixing  beneath  an  ice-covered  ocean  indicates  that  the  initial  w  ind  event, 
rather  than  subsequent  inertial  shear,  causes  the  majority  of  the  mixing.  This  work  points 
towards  episodic  wind-forced  shear  at  the  base  of  the  mixed  layer  coupled  with  static 
instability  from  brine  rejection  due  to  ice  production  as  a  major  factor  in  mixing  on  the 
wAP  shelf. 

Thesis  Supervisors:  Robert  C.  Beardsley  and  W.  Brechner  Owens 
Senior  Scientists,  Department  of  Physical  Oceanography 
Woods  Hole  Oceanographic  Institution 
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1.  Chapter  1 


Introduction 

1.1.  Overview 

Icc  in  the  world  oceans  plays  an  important  role  in  regulating  the  Earth’s  climate. 
Sea  icc  acts  to  insulate  the  relatively  warm  ocean  from  the  cold  polar  atmosphere. 
Furthermore,  movement  of  icc  from  the  formation  to  the  melt  region  acts  as  a  transport  of 
latent  heat  and  freshwater.  This  movement  depends  on  a  complicated  balance  between 
the  winds,  the  surface  ocean  currents  and  the  dynamics  of  the  sea  icc  itself.  Sea  ice  is 
also  believed  to  play  an  important  role  in  biological  processes.  The  importance  of  the 
polar  oceans  in  climate  studies  combined  with  a  rapid  increase  in  available  technology  to 
access  and  measure  ice-covered  seas  have  prompted  a  renewed  interest  in  understanding 
sea  icc  and  the  ocean  below  it. 

This  interest  in  high  latitude  physical-biological  processes  led  to  the  Southern 
Ocean  Global  Ocean  Ecosystems  Dynamics  (SO  GLOBEC)  program,  which  chose  the 
continental  shelf  of  the  western  Antarctic  Peninsula  (wAP)  for  study  due  to  its  unusually 
high  concentration  of  Antarctic  Krill  {Euphasia  superba)  and  krill  predators  such  as 
penguins,  seals  and  whales  (Figure  1 . 1. 1 ).  SO  GLOBEC  aims  to  determine  what 
conditions  favor  the  reproduction,  recruitment,  and  retention  of  krill  populations, 
including  icc  cover  (high  icc  versus  low  icc  years),  on-shelf  intrusion  of  warmer  nutrient 
rich  waters  and  subsequent  vertical  mixing. 
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Figure  1.1.1  Left  panel  shows  the  Marguerite  Bay  region  of  the  western  Antarctic  Peninsula  shelf, 
depth  is  in  meters.  White  squares  and  circles  show  the  locations  of  the  moorings  and  weather 
stations,  respectively.  The  thick  black  line  represents  the  coast  and  approximate  position  of  the 
permanent  ice  shelf.  The  depth  contours  end  at  the  ltt00-m  isobath,  the  approximate  location  of  the 
shelf  break.  The  right  panel  shows  the  location  of  the  study  region  in  (he  larger  context,  w  ith  South 
America  to  the  north. 

The  physical  oceanographic  environment  of  the  wAP  shelf  has  a  number  of 
features  which  differentiate  it  from  a  typical  mid-latitude  shelf,  c.g.,  the  North  American 
Middle  Atlantic  Bight.  The  wAP  shelf  is  nearly  fully  ice-covered  in  winter,  and  almost 
totally  ice-free  in  summer.  Most  of  the  shelf  is  quite  deep  (-500  m)  with  a  -800  m  deep 
trough  running  from  Marguerite  Bay  nearly  to  the  shelf  break.  The  Antarctic 
Circumpolar  Current  (ACC)  flows  northeastwards  along  the  shelf  break  carrying 
relatively  warm  and  salty  Circumpolar  Deep  Water  (CDW). 

The  surface  waters  above  about  100  m,  known  as  Antarctic  Surface  Water 
(AASW),  show  a  strong  seasonal  signal,  becoming  warmer  in  summer  as  a  result  of  the 
positive  surface  heat  flux  (Figure  1.1.2).  In  addition,  the  wAP  shelf  surface  waters 
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freshen  slightly  as  a  result  of  glacier  melt,  coastal  runoff  and  to  a  lesser  extent  ice  melt 
(Dierssen  ct  al.,  2002;  Meredith  ct  al.,  2006).  In  fall,  negative  surface  heat  flux  causes 
rapid  cooling  of  the  surface  waters  that  become  saltier  as  they  mix  down  into  the  remnant 
Winter  Water  (WW)  and  brine  is  rejected  during  the  formation  of  sea  ice  (Smith  ct  al, 
1999).  On  intcrannual  timescales,  ice  production  and  subsequent  brine  rejection  has  been 
found  to  be  the  main  control  on  upper  ocean  stratification  (WW  production  and 
maintenance)  in  northern  Marguerite  Bay  off  the  British  Antarctic  Survey  station  Rothcra 
(Meredith  ct  al.,  2004).  Relatively  warm  and  salty  Circumpolar  Deep  Water  (CDW) 
flows  onto  the  shelf  below  the  AASW  and  WW,  and  is  divided  by  its  properties  into  an 
Upper  and  Lower  Circumpolar  Deep  Water  (UCDW  and  LCDW).  Despite  strong 
atmospheric  forcing,  waters  from  the  wAP  shelf  do  not  become  dense  enough  to  form 
Antarctic  Bottom  Water  (Klinck,  1998;  Hofmann  and  Klinck,  1998). 
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Figure  1*1.2  Four  CTD  casts  on  the  wAP  shelf  showing  water  masses  and  the  seasonal  progression  of 
the  temperature  and  salinity  in  (a)  a  Temperature-Salinity  plot  and  <b)  vertical  profiles  of  potential 
temperature,  salinity  and  potential  density. 
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Air-sea  heat  flux  is  one  to  two  orders  of  magnitude  larger  over  open  water  than 
over  thick  (~  1  m)  sea  ice.  Wind  stress  on  sea  ice  is  frequently  balanced  in  large  part  by 
internal  ice  stresses  rather  than  imparting  momentum  to  the  sea  ice  and  upper  ocean  (e.g. 
Chapter  4  on  wAP  shelf;  Uotila  et  al.,  2000  in  the  Weddell  Sea;  Steele  et  al.,  1997  in  the 
Arctic).  The  strong  influence  sea  ice  has  on  the  vertical  transport  of  heat  and  momentum 
between  the  atmosphere  and  the  ocean  plays  a  major  role  in  the  hydrography  of  the  wAP 
shelf  and  is  a  central  theme  of  this  thesis. 

The  data  used  in  this  thesis  represent  the  most  comprehensive  set  of  physical 
oceanographic  measurements  made  in  Marguerite  Bay  to  date.  It  consists  primarily  of 
two  years  of  moored  oceanographic  and  Automatic  Weather  Station  (AWS)  data  which 
cover  the  winter  ice  seasons  of  2001-2003.  During  the  ice-free  austral  summer, 
subsurface  moorings  were  deployed  and  remained  beneath  the  ice  for  the  duration  of  the 
ice-covered  austral  winter,  including  the  onset  and  breakup  of  sea  ice.  They  had  an 
upward-looking  acoustic  Doppler  current  profiler  (ADCP)  and  measured  temperature, 
salinity,  pressure  and  velocity  at  various  depths  (Figure  1.1.3).  The  locations  ol  the 
moorings,  as  well  as  the  two  AWSs  are  shown  in  Figure  1.1.1.  Details  of  the  mooring 
and  weather  station  data  can  be  found  in  Moffat  et  al.  (2005).  Satellite  data  and 
Conductivity-Temperature-Depth  (CTD)  profile  data  collected  on  mooring  and 
hydrographic  cruises  were  used  as  well. 
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Figure  LI *3  Schematic  of  the  C3  mooring  deployed  during  SO  GLOBKC  in  Marguerite  Trough  in 
the  mouth  of  Marguerite  Bay  {68"  6.0'  S,  70°  31.8*  W)  (Moffat  et  aL  2005). 

1.2.  Outline 

It  is  noted  that  Chapter  2  has  been  published  (Howard  et  al.,  2004)  and  Chapters  3 
and  4  have  been  submitted  for  publication  (Hyatt  et  al.,  2006a  and  2006b), 

In  Chapter  2,  the  small-scale  processes  that  contribute  to  upward  diapycnal  iluxcs 
of  heat,  salt  and  nutrients  from  the  UCDW  to  the  surface  mixed  layer  are  discussed  and 
initial  flux  estimates  based  on  published  parameterizations  are  made.  Near-inertial  shear 
across  the  sharp  pycnocline  at  the  base  of  the  mixed  layer  is  observed  and  considered  in 
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the  context  of  vertical  mixing.  A  previous  suggestion  that  double-diffusive  convection 
(DDC)  provides  significant  upward  heat  fluxes  (of  order  10  W  m' )  in  a  nearby  region 
(Smith  and  Klinck,  2002)  is  not  supported  by  the  analyses  in  Chapter  2,  which  indicate 
almost  no  contribution  (0.2-0. 4  W  m~2)  to  diapycnal  fluxes  from  DDC. 

In  Chapter  3,  a  technique  for  the  analysis  of  data  from  a  moored  upward-looking 
ADCP  to  determine  ice  coverage,  draft  and  velocity  is  presented  and  applied  to  data 
collected  on  the  wAP  shelf.  This  method  provides  sea  ice  information  when  no  dedicated 
upward- looking  sonar  (ULS)  data  are  available. 

In  Chapter  4,  the  onset,  draft,  motion  and  dynamics  of  sea  ice  within  Marguerite 
Bay  arc  characterized  for  the  two  ice  seasons  spanning  2001-2003  with  data  derived  in 
part  using  the  methods  described  in  Chapter  3.  Both  years  have  roughly  seven  months  of 
nearly  complete  ice  cover,  but  the  timing  is  quite  different  between  the  two  years.  The 
ice  draft  is  considerable,  over  3  m  lasting  for  weeks  during  both  winters.  Internal  ice 
stresses  are  shown  to  be  important  in  explaining  the  observed  motion  of  the  ice  pack. 
Tidal  and  near-inertial  band  motions  are  discussed  as  well. 

Vertical  fluxes  of  heat,  salt  and  nutrients  in  a  seasonally  ice-covered,  high- latitude 
coastal  ocean  have  become  a  central  question  to  the  SO  GLOB  EC  program.  To  address 
this,  in  Chapter  5  a  one-dimensional  model  of  vertical  mixing  is  used  to  investigate  two 
potentially  important  mixing  processes,  namely  shear  caused  by  a  wind  event  and 
subsequent  inertial  oscillations;  and  convection  due  to  brine  rejection  associated  with 
freezing  as  leads  open  and  close  with  divergent  and  convergent  inertial  motions.  The 
model  results  indicate  that  the  wind  event  is  responsible  for  the  majority  of  heat  flux 
across  the  base  of  the  mixed  layer,  not  the  subsequent  inertial  oscillations.  When  forced 
by  observed  atmospheric  and  ice  conditions  on  the  wAP  shelf,  large  (one-day  average  of 
-  300  W  m 2)  fluxes  across  the  base  of  the  mixed  layer  are  produced  during  a  strong 
forcing  event.  However,  averaging  on  an  annual  timescale  predicts  a  smaller  flux,  on  the 
order  of  1 0  W  m'2. 

Finally,  Chapter  6  provides  a  summary  of  the  thesis  and  discusses  its  implications 
and  possible  future  work. 
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2.  Chapter  2 


Mixing  in  the  pycnocline  over  the  western  Antarctic 
Peninsula  shelf  during  Southern  Ocean  GLOBEC 

2.1.  Introduction 

The  western  Antarctic  Peninsula  (WAP)  shelf  region  (Figure  2. 1 . 1 )  is  an  area  of 
unusually  high  krill  production  and  thus  a  favorable  habitat  for  krill  predators  such  as 
Adelie  penguins  and  Crabeater  seals.  The  Southern  Ocean  Globat  Ecosystem  Dynamics 
(SO  GLOBEC)  program  examined  the  unique  physical  and  biological  factors  that 
contribute  to  enhanced  krill  growth,  reproduction,  recruitment  and  survivorship  in  this 
region  throughout  the  year  (Hofmann  et  al.,  2004). 
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Longitude 


Figure  2.1.1.  Map  of  the  study  area  surrounding  Marguerite  Bay  on  the  western  side  of  the 
Antarctic  Peninsula.  Thin  black  lines  are  the  500  and  2000  m  isobaths  based  on  the  Smith  and 
Sandwcll  ET0P02  dataset  updated  with  GLOBEC  bathymetry  data.  Bold  black  lines  indicate 
fronts  of  major  ice  shelves.  The  location  of  Surv  ey  Station  3  near  66.2°S,  69.1°\V  (see  Fig.  2.1.2)  is 
indicated  by  a  star,  while  CTD  76  near  69"S,  75°VV  from  the  fall  2001  survey  cruise  (see  Fig.  2.3.3)  is 
indicated  by  a  filled  triangle.  The  drift  track  for  the  section  of  the  winter  process  cruise  on  L.  M. 
Could  used  in  Fig.  2.3.1  is  also  indicated  (centered  near  67.2“S,  70.5°\V). 

In  this  chapter  we  identify,  and  attempt  to  quantify  the  fluxes  associated  with,  the 
small-scale  physical  oceanographic  processes  that  facilitate  the  supply  of  heat,  salt,  and 
nutrients  to  the  upper  ocean  over  the  WAP  shelf.  The  principal  pathway  for  these  fluxes 
to  the  shelf  is  via  large-scale  lateral  subsurface  intrusions  of  warm,  salty,  nutrient-rich 
Upper  Circumpolar  Deep  Water  (UCDW).  The  UCDW  is  derived  from  the  Antarctic 
Circumpolar  Current  (ACC)  (Prezelin  et  al.,  2000;  Smith  and  Klinck,  2002)  and  flows 
shoreward  along  bathymetric  features  such  as  the  Marguerite  Trough  (Klinck  et  al., 
2004).  Over  the  shelf,  the  UCDW  typically  lies  below  about  200  m,  well  below  the  less 


18 


dense  surface  waters  that  are  formed  from  ocean/atmosphere  buoyancy  exchange  and 
sea-ice  processes.  In  a  plot  (not  shown)  of  potential  temperature  (0)  against  salinity  (5), 
the  transition  region  between  the  UCDW  and  the  surface  waters  (roughly  the  depth  range 
80-200  m  in  Figure  2.1.2)  is  fairly  linear,  suggesting  that  the  0-S  relationship  in  this  depth 
range  is  set  by  diapycnal  mixing  between  the  end-member  water  types.  Our  goal  is  to 
determine  the  processes  that  drive  this  mixing. 


Figure  2.1.2.  Profiles  of  temperature  (left)  and  salinity  (right)  at  survey  station  3  from  the  fall  (gray 
lint1)  and  winter  (black  line)  survey  cruises  (see  Fig.  2.1.1  for  location).  The  profiles  were  separated 
by  -H9  days.  Upper  Circumpolar  Deep  Water  (UCDW)  that  has  intruded  onto  the  shelf  occupies  the 
water  column  below  about  200  m*  On  the  assumption  that  lateral  advert  ion  is  negligible  (Smith  and 
Klinck,  2002),  the  heat  and  salt  flux  divergences  associated  with  the  change  in  heat  and  salt  content 
of  the  profiles  above  110  m  (indicated  by  the  horizontal  dashed  line)  are  -72  W  m  2  and  0,7  mg  m"2  s'1, 
respectively. 

The  most  likely  candidates  for  pyenocline  mixing  processes  in  this  environment 
are  shear-driven  instability  (SDI)  and  double-diffusive  convection  (DDC).  In  turn,  the 
most  likely  sources  of  energy  for  SDI  are  from  wind  stress  and  tides  (see,  c.g.,  Munk  and 
Wunsch,  1998).  The  mixing  rates  associated  with  these  processes  can  be  estimated  in 
three  ways:  a  box-model  of  heat  and  salt  budget  calculations  for  the  surface  mixed  layer 
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(e.g..  Smith  and  Klinck,  2002);  measurements  of  scalars  or  velocity  shear  that  fully 
resolve  the  turbulent  microscales  of  order  1*10  cm  (see  Gregg,  1987);  or  by  invoking 
parameterizations  that  relate  mixing  rates  to  measured  “large-scale”  quantities. 

The  difficulty  of  a  box-model  approach  for  the  WAP  region  is  illustrated  in 
Figure  2. 1 .2.  For  a  representative  site  (survey  station  3;  see  Figure  2.1.1  for  location) 
over  the  central  WAP  shelf,  the  change  in  upper-ocean  heat  and  salt  content  between  the 
austral  fall  and  winter  2001  survey  cruises  is  consistent  with  net  vertical  flux  divergences 
of  heat  and  salt  across  the  mixed  layer  of -72  W  m 2  and  0.7  mg  m 2  s',  respectively.  In 
contrast,  the  mean  upward  fluxes  from  below  the  mixed  layer  are  likely  of  order 
1-10  W  m 2  and  0.01  mg  m 2  s'1  (Smith  and  Klinck,  2002),  implying  that  the  upper  ocean 
heat  and  salt  budgets  are  dominated  by  the  seasonal  cycles  of  heat  and  fresh  water 
exchange  at  the  ocean  surface.  Attempts  to  estimate  the  diffusivities  in  the  pycnocline 
from  box-model  studies  of  surface  mixed  layer  heat  and  fresh  water  budgets  are  therefore 
subject  to  the  significant  errors  associated  with  parameterizations  of  the  large  surface 
fluxes. 

No  microstructure  data  are  available  from  SO  GLOBEC  year  I  cruises,  hence  the 
direct  measurement  of  diffusivity  in  the  pycnocline  is  not  possible.  Instead,  for  the 
present  chapter  we  concentrate  on  the  third  approach;  using  fine-scale  hydrographic  and 
velocity  data  to  identify  mixing  processes  in  the  pycnocline.  The  data  we  analyze  consist 
of  hydrographic  profiles  from  conductivity-temperature-depth  (CTD)  recorders  and 
acoustic  Doppler  current  profilers  (ADCPs).  The  CTD  resolution  is  -1  m  when  sea-state 
is  not  severe,  and  the  ADCP  collects  velocity  profiles  at  a  vertical  resolution  of  ~10  m. 

We  first  describe  the  relevant  data  sets  (Section  2.2),  then  review  the  processes  of 
SDI  and  DDC  and  estimate  their  contribution  to  pycnocline  mixing  through 
parameterizations  that  are  commonly  employed  in  numerical  models  (Section  2.3). 
Section  2.4  provides  a  summary  and  discussion  of  our  results. 

2.2.  Data 

The  SO  GLOBEC  program  collected  data  in  the  WAP  area  during  the  austral  fall 
and  winter  seasons  of  200 1  and  2002  using  4  survey  cruises  and  4  process  cruises.  Three 
additional  cruises  to  deploy,  service  and  recover  moorings  took  place  in  the  austral  late 
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summers  of  200 1 ,  2002,  and  2003.  For  this  study  we  concentrate  on  data  front  CTD 
profiles  and  the  vessel-mounted  ADCP  during  the  fall  and  winter  2001  survey  cruises  on 
the  RVIB  N.  B.  Palmer  (U.S.  SO  GLOBEC,  2001  b,c)  and  the  winter  2001  process  cruise 
on  the  ARSV  L.  M.  Gould  (U.S.  SO  GLOBEC,  2001c)  (Table  1).  Vertical  profiles  of 
temperature  (7)  and  conductivity  (C)  as  functions  of  pressure  ( P )  were  obtained  on  a  grid 
pattern  (Figure  2.2.1)  on  the  survey  cruises  using  a  SeaBird  SBE  91  V  CTD  system. 
Additional  CTD  profiles  were  obtained  from  the  L.  M.  Gould  process  cruise  using  a 
similar  CTD  system.  Salinity  (,$)  and  potential  density  (oH)  were  derived  front  C,  T,  and 
P.  Accuracies  of  T and  5  are  -1 0”1  °C  and  —1 0  1  practical  salinity  units  (psu), 
respectively,  after  post-calibration.  The  data  that  we  use  here  were  averaged  initially 
over  1-nt  vertical  bins:  additional  averaging  used  in  our  analysis  is  described  where 
necessary.  The  cruise  reports  (U.S.  SO  GLOBEC  2001a.b,c)  and  Klinck  ct  al.  (2004) 
provide  more  details  about  the  CTD  data  set  and  the  general  hydrographic  structure  and 
circulation  that  it  reveals. 
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Figure  2*2*1.  Left:  Locutions  of  CTD  stations  during  the  2001  SO  GLOBEC  survey  cruises*  Fall 
stations  are  shown  as  plus  symbols  (+).  Winter  stations  are  shown  as  circles  (o)*  Right:  Location  of 
useable  vessel-mounted  ADCP  data  collected  by  the  RV  IB  /V*#*  Palmer  for  the  same  cruises*  Kail 
survey  data  are  shown  as  gray  dots.  Winter  survey  data  are  shown  as  black  dots*  Large  symbols 
indicate  the  locations  of  the  same  sites  shown  in  Fig*  2*1. L  The  cross-shelf  transect  I  (Fig,  2.3.2),  is 
indicated  by  an  arrow. 


Vertical  profiles  of  upper-ocean  horizontal  current  (U(z)=u{z)i+v(z)j)  were 
obtained  from  each  vessel  with  RDI  150  khlz  vessel-mounted  ADCPs.  The  ADCP 
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systems  were  set  for  8-m  vertical  bins,  8-m  pulse  length,  and  an  ensemble  averaging 
period  of  5.0  minutes  for  all  SO  GLOBEC  cruises  except  the  2001  mooring  deployment 
and  2001  fall  process  cruises  on  the  L  M.  Gould.  These  two  cruises  used  a  16-m  pulse 
length.  The  resultant  nominal  accuracy  of  each  estimate  of  U(z)  is  -  I  cm  s'1,  per  the 
manufacturer's  design  error  estimate.  The  ADCP  on  the  N.  B.  Palmer  measures  currents 
from  3 1  m  to  ~300-350  m  depths.  Quoted  depths  refer  to  the  center  of  the  8-m  vertical 
bins.  The  ADCP  on  the  L.  M.  Gould  measures  currents  from  26  m  to  -300-350  m  depths. 
Navigation  data  were  obtained  with  an  Ashtech  3D  GPS  system  and  were  used  in  the 
processing  of  the  ADCP  data.  Bottom  tracking  was  used  when  the  ship  was  in  water  less 
than  500  m  deep.  Figure  2.2.1  (right  panel)  shows  the  location  of  usable  ADCP  data 
from  the  2001  N.  B.  Palmer  survey  cruises.  Sea  state  and  ice  cover  both  interfered  with 
collection  of  ADCP  data.  Sea-state  effects  can  be  seen  as  gaps  in  the  grid  for  the  fall 
survey  cruise.  Sea  ice  interferes  with  ADCP  data  collection  while  the  ship  is  underway. 
During  the  winter  cruise,  therefore,  data  were  acquired  only  as  time  series  at  fixed  sites 
while  the  vessel  was  stopped  during  CTD  and  other  station  operations.  These 
measurements  lasted  for  as  long  as  the  on-station  work,  which  typically  took  a  few  hours 
depending  on  water  depth  and  sampling  activities. 


Year  1  Cruises 

Type 

Dale 

NBP0I-03 

Fall  survey 

Apr.  23  -  June  6 

NBP0I-04 

Winter  survey 

July  22  -  Aug.  31 

LMG01-06 

Winter  process 

July  21  -  Aug.  31 

Table  2.1.  Summary  of  201)  1  SO  GI.GBEC  cruises  discussed  in  this  chapter.  Dates  and  types  of 
cruises  are  included.  LMC  -  ARSY'  L.  M.  Gould,  NBP  -  RVIB  N.  B.  Palmer 

2.3.  Vertical  mixing  processes 

2.3.1.  Shear-Driven  Instability 

Shear-driven  instability  (SDI)  is  the  primary  source  of  mixing  in  the  World  Ocean 
pycnocline.  The  source  of  energy  for  SDI  is  a  combination  of  atmospheric  forcing 
(primarily  wind  stress,  but  also  surface  buoyancy  fluxes)  and  tides:  Munk  and  Wunsch 
(1998)  and  Wunsch  (2000)  suggest  that  about  50%  of  the  globally  integrated  mixing  is 
associated  with  each  source.  This  ratio  is,  however,  strongly  dependent  on  local 
conditions.  Tides  can  be  the  dominant  source  of  mixing  via  production  of  baroclinic 
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tides  and  higher  frequency  internal  waves  in  regions  where  strong  barotropic  tidal 
currents  flow  across  steep  topography;  see,  e.g.,  Padman  and  Dillon  ( 1991)  and 
Wijesekcra  et  al.  (1993)  for  high-latitude  examples.  The  spatial  distribution  of  shear  is 
complicated  by  the  propagation  of  super-inertial  waves  away  from  the  energy  source 
(see,  e.g.,  Dushaw  et  al.,  1995;  Muench  et  al.,  2002;  Lozovatsky  et  al.,  2003).  Thus  for 
any  particular  location,  determining  the  likely  strength  of  mixing  through  SDI  requires 
knowledge  not  only  of  the  local  strength  of  forcing  of  baroclinic  waves  but  also  the 
potential  for  wave  energy  to  propagate  from  elsewhere  to  the  observation  site. 

For  reasons  explained  below,  we  believe  that  the  dominant  source  of  shear  in  the 
WAP  pycnocline  is  associated  with  wind-forced  near-inertial  waves.  Baroclinic  tides  may 
provide  additional  shear  in  certain  localized  regions.  We  discuss  each  of  these  shear 
sources  separately. 

2.3.1. 1.  Near-inertial  waves 

At  process  site  3  from  the  winter  process  cruise  while  the  L.  M.  Gould  drifted  in 
pack  ice,  we  observed  persistent  shear  in  the  ADCP  data,  with  the  shear  being 
concentrated  across  the  sharp  pycnocline  at  the  base  of  the  mixed  layer  (Figure  2.3.1). 
The  dominant  signal  in  this  record  is  a  semidiurnal  oscillation,  where  the  current  vectors 
in  the  layers  above  and  below  the  sharp  pycnocline  rotate  anticlockwise  but  are  roughly 
180°  out  of  phase  with  each  other.  Current  speeds  are  up  to  24  cm  s'1  in  the  upper  mixed 
layer,  and  weaker  but  still  significant  in  the  lower  layer.  Although  the  ADCP  record  does 
not  reach  full  depth  (-500-650  m)  in  this  region,  the  impression  is  that  the  net  depth- 
integrated  flow  is  small,  with  the  volume  flux  in  the  faster  and  thinner  surface  layer  being 
compensated  for  by  the  opposing  flux  in  the  thicker  but  slower-moving  lower  layer. 
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Figure  2.3.1.  Density  and  velocity  fields  for  process  site  3  front  the  winter  process  cruise  on  L  M. 

(it) u Id  (see  Fig.  2.1. 1  for  location).  Vertical  profile  of  potential  density  for  CTD  227  (upper  left),  time 
series  of  shear  magnitude  <[L.,|>  (s  ')  calculated  front  ADCP  data  (upper  right),  time  series  of  current 
speed  (cm  s'1)  at  50  in  depth  (black  line)  and  the  negative  of  the  speed  at  154  m  depth  (gray  line)  from 
the  vessel-mounted  ADCP  (central  right),  and  time  series  of  current  direction  for  50  m  depth  (black 
line)  and  154  m  depth  (gray  line)  (tower  right).  The  solid  vertical  lines  on  the  right-hand  panels  show- 
the  time  of  CTD  227.  The  fundamental  periodicity  of  the  phase  lines  is  -12,9  h  (-/),  and  the  surface 
layer  is  roughly  180  degrees  out  of  phase  with  the  deeper  layer  for  the  entire  time  period. 

The  duration  of  this  record  at  a  fairly  constant  location  allows  us  to  identify  the 
semidiurnal  signal,  on  the  basis  of  its  dominant  periodicity,  ellipticity  and  rotation 
direction  (counter-clockwise),  as  being  near-inertial  (period  -2 n/f  s;l 3.0  h  at  67°S)  rather 
than  tidal:  the  principal  semidiurnal  tidal  constituents,  M;  and  S2,  have  periods  of 
-12.42  h  and  12.00  h,  respectively.  Analyses  of  satellite-tracked  drifter  paths  (Beardsley 
et  al.,  2004)  indicate  that  near-inertial  surface  currents  are  frequently  energetic  over  the 
WAP  shelf  and  into  the  mouth  of  Marguerite  Bay.  Near-inertial  oscillations  are 
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commonly  associated  with  rapid  changes  in  wind  stress,  i.e.,  the  passage  of  storm  fronts, 
and  such  events  are  frequent  in  this  region  (U.S.  SO  GLOBEC,  200la,b,c).  The  shear 
across  the  sharp  pycnocline  associated  with  the  near-inertial  currents  is  strong  (f  igure 
2.3.1),  suggesting  that  these  waves  may  provide  sufficient  energy  for  SDI  mixing  to  be 
significant. 


2.3.1. 2.  Tides 

A  second  plausible  source  of  shear  and  turbulence  for  the  pycnocline  in  the  study 
area  is  from  tides,  specifically  baroclinic  tides  generated  as  the  barotropic  tide  flows 
across  sloping  topography;  see,  e.g.,  Levine  et  al.  (1997)  and  Muench  et  al.  (2002)  for 
two  Antarctic  examples.  The  two  principal  tidal  species,  semidiurnal  and  diurnal,  behave 
very  differently  in  the  Antarctic  (Foldvik  et  al.,  1990):  in  the  WAP  region,  u>tj &>/  for  al! 
semidiurnal  tidal  harmonics,  whereas  o></for  the  diurnal  tides.  Super-inertial  baroclinic 
semidiurnal  tides  can  be  easily  generated  and  propagate  freely  away  from  the  generation 
sites,  whereas  the  sub-inertial  diurnal  tides  are  much  more  restricted  in  their  ability  to 
support  shear  and  will  be  constrained  to  topographic  wave  guides  such  as  the  shelf  break. 

There  are  few  existing  current  meter  data  sets  for  the  WAP  region,  hence  we  base 
our  knowledge  of  tides  on  a  numerical  model,  the  Circum- Antarctic  Tidal  Simulation, 
version  02.01  (CATS02.01)  (Padman  et  al.,  2002).  This  model  predicts  that  the 
barotropic  semidiurnal  tidal  kinetic  energy  on  the  WAP  shelf  is  very  small  Indeed, 
Beardsley  et  al.  (2004)  show  that,  over  the  mid-shelf,  predicted  tidal  current  speeds  are 
only  -1/10  of  the  surface  near-inertial  currents.  Thus,  it  is  unlikely  that  significant  shear 
can  be  generated  by  semidiurnal  tides,  and  this  hypothesis  is  confirmed  by  a  primitive- 
equation  numerical  model  that  we  have  run  of  the  three-dimensional  structure  of  M:  tides 
for  the  region  offshore  of  Marguerite  Bay.  In  contrast,  CATS02.01  predicts  that  the  time- 
averaged  diurnal  tidal  current  speed  exceeds  10  cm  s"1  at  a  few  locations  along  the  WAP 
shelf  break.  Maximum  tidal  current  speeds  (“spring”  tides)  usually  exceed  mean  speeds 
by  a  factor  of  ~2.  The  strong  currents  are  associated  with  diurnal-band  topographically 
trapped  vorticity  waves  (DTVWs),  which  are  known  to  have  a  significant  cross-slope 
component  over  the  upper  slope  and  shelf  break  in  many  locations  around  Antarctica 
(Middleton  et  al„  1987;  Padman  and  Kottmeier,  2000;  Padman  et  al.,  2003).  For  a 
variety  of  reasons  the  energy  of  DTVWs  are  difficult  to  accurately  model  (Foldvik  et  al.. 
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1990;  Padman  and  Kottmcier,  2000;  Cummins  et  al..  2000)  and  so  we  only  use 
CATS02.01  to  indicate  the  approximate  magnitude  of  currents.  Some  support  for  the 
model  is  provided,  however,  by  observations  of  the  diurnal -band  variability  of  velocity 
evaluated  from  surface  drifters  (Beardsley  et  al.,  2004).  The  diumal-band  velocity  of  one 
of  these  drifters  that  spent  significant  time  over  the  outer  continental  shelf  was  up  to 
-15  cm  s"1,  broadly  consistent  with  model  predictions.  These  data,  and  studies  of 
DTVWs  in  the  Arctic  Ocean  (Padman  and  Dillon,  1991 ;  Padman  et  al.,  1991;  Padman  et 
al.,  1992),  support  our  hypothesis  that  barotropic  tidal  energy  in  the  form  of  DTVWs  may 
result  in  mixing  near  the  shelf  break.  Further  work  is  required  to  determine  whether  a 
significant  fraction  of  the  energy  in  the  primarily  barotropic  DTVWs  can  be  transferred 
to  higher  frequency  baroclinic  motions  capable  of  inciting  the  growth  of  shear 
instabilities  in  the  pycnocline.  While  the  previously  cited  studies  support  this  energy 
pathway  to  mixing  for  the  Arctic  Ocean,  the  weaker  stratification  of  the  WAP  region  may 
preclude  significant  baroclinic  wave  generation  associated  with  DTVWs. 

2.3. 1.3.  Mixing  Rates  Associated  with  Shear 

In  the  absence  of  microstructure  measurements,  estimates  of  pycnocline  mixing 

rates  due  to  SDI  can  only  be  made  through  parameterizations  based  on  larger-scale 
features  of  the  hydrographic  and  velocity  structure.  Such  parameterizations  fall  into  two 
broad  categories;  schemes  based  on  Richardson  number  Ri  (e.g.,  Pacanowski  and 
Philander,  1981;  hereinafter  called  PP-81),  and  those  based  on  energy  flux  to  smaller 
scales  via  wave- wave  interactions  in  an  internal  wave  field  frequency-wavenumber 
spectrum  (e.g.,  Gregg,  1989;  Wijesekera  et  al.,  1993;  Levine  et  al.,  1997). 
Parameterizations  of  vertical  diffusion  in  numerical  models  tend  to  be  based  on  Ri ,  e.g.. 
Large  et  al.  (1994).  However,  the  relationship  between  Ri  and  the  turbulent  diffusivity  A'v 
in  the  ocean  is  not  well  constrained  (Peters  et  al.,  1988),  and  any  A/'-based  estimate  of  Kv 
will  have  large  errors  associated  with  it.  Nevertheless,  it  is  instructive  to  see  what  value 
of  A'v  would  be  predicted  by  physical  or  coupled  physical-ecological  models  that  attempt 
to  reproduce  the  SO  GLOBEC  data  sets. 

The  Richardson  number  is  defined  as  N2l\d\ildzf,  where  V  is  the  buoyancy 
frequency  and  U(z)=i/(z)i+v(z)j  is  the  vertical  profile  of  horizontal  velocity.  We  define  a 

“mean”  gradient  Richardson  number  Ri ,  as 
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(2.1) 


Ri  = 


where  angle  brackets  (...)  indicate  averaging  in  space  and/or  time  and 

=  {duldzf  +  ( dv/dz ): .  The  overbar  on  Ri  reminds  the  reader  that  this  parameter  is 

based  on  averaged  (jV2)  and  (u^ .  Values  of  Ri  based  on  fully  resolved  density  and 


velocity  fields  can  be  used  to  identify  when  the  necessary  condition  for  dynamic 
instability  of  the  shear  flow  is  met.  Values  of  Ri< 0.25  are  generally  taken  to  indicate  that 
shear  instabilities  are  likely  to  occur  according  to  linear  theory  for  two-dimensional 
stratified  shear  flow,  while  Ri<\  indicates  that  the  potential  exists  for  advective 

instabilities  (Orlanski  and  Bryan,  1969).  The  interpretation  of  Ri  is,  however,  intimately 
related  to  the  vertical  and  temporal  averaging  scales  for  density  and  velocity  gradients 

(Padman  and  Jones,  1985).  In  this  view,  Ri  is  some  measure  of  the  probability  of  shear 
instabilities  occurring  at  time  and/or  space  scales  smaller  than  the  averaging  scales,  not  a 
threshold  criterion  separating  fully  turbulent  from  laminar  flow. 


Velocity  shear  (Ur=<3U/&)  is  approximated  in  the  ADCP  ensemble  data  as 

'A(ii)  A(v)^ 

A 2  5  Az 


<«.)- 


(2.2) 


where  A(h)  and  A(v)  are  the  differences  in  east  and  north  velocity  components  between 
adjacent  ADCP  depth  bins  separated  by  Az  m.  As  before,  angle  brackets  (...)  denote 
averaging  over  space  and/or  time  scales.  Here,  u  and  v  are  averaged  over  an  ADCP  depth 
bin  that  is  Az  (8  m)  thick  and  over  times  as  described  herein,  typically  of  order  1  h. 

Large  values  of  shear  magnitude,  |(UZ)|,  were  frequently  found  in  the  upper  water 


column  on  the  WAP  shelf.  Figure  2.3.2  shows  a  cross-shelf  plot  of  |(U*)|  for  transect  line 
1  from  the  fall  survey  cruise  on  N.  B.  Palmer  (see  Figure  2.2. 1  for  location).  The  “on- 
station"  data  were  removed  from  the  total  transect  record  prior  to  making  this  figure. 
There  is  a  strong  shear  layer  (|Uj-|>0.01  s  ')  in  the  70-80  m  depth  range,  coincident  with 
the  strong  stratification  at  the  base  of  the  surface  mixed  layer.  Throughout  the  study 
region,  the  strongest  shears  tended  to  occur  in  this  portion  of  the  water  column,  although 
there  arc  also  several  occurrences  of  relatively  strong  shear  in  the  deeper  pycnocline. 
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The  (ransect  shown  in  Fig.  2.3.2  took  ~26  h  to  complete  (~2  inertial,  and  ~l  or  ~2  tidal 
periods)  and  so  the  plotted  variability  is  a  combination  of  both  spatial  and  temporal 
changes  in  currents. 
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Figure  2,3*2,  Transect  of  shear  magnitude  (|U*[)  (s’1)  calculated  from  ADCP  data  collected  on  the  fall 
survey  cruise  Transect  l  (see  Fig,  2,2*  1  for  location).  Station  locations  are  shown  at  the  top  of  the 
plots*  Station  I  is  at  the  shelf  break  and  Station  5  is  nearest  to  the  coast*  The  transect  took  -2h  h  to 
complete. 


We  can  only  calculate  Ri  from  (I )  with  on -station  data,  where  measurements  of 
(N2)  and  (U:  ^  arc  simultaneously  available*  The  value  of  {N2)  in  (1)  is  based  on 

potential  density  a^(z)  from  a  single  CTD  profile.  For  consistency  with  the  estimation  of 
we  first  average  oh  in  the  same  8-m  depth  bins  used  for  ADCP-based  velocity, 

-j 

then  calculate  (N~)  from 


Nz)  = 


g 


A z 


(2.3) 


where  A(ae)  is  the  potential  density  difference  between  adjacent  8-m  ADCP-based  depth 
bins  and  p„  is  a  characteristic  density.  The  averaging  time  for  (if ;  )  is  taken  as  the 

duration  of  the  CTD  cast,  up  to  a  maximum  of  2.7  h  centered  on  the  mean  time  of  the 
cast.  With  Az=8  m  and  A/=0.5-2.7  h,  we  cannot  fully  resolve  the  expected  shear  and 
density  gradient  joint  frequency- wavenumber  spectrum  {see,  e.g.,  Wijesekera  et  al„ 
1993).  These  averaging  scales  are,  however,  of  the  same  order  of  magnitude  as  the  time 
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step  and  vertical  resolution  in  regional  numerical  models,  and  thus  our  interpretation  of 
mixing  based  on  Ri  emulates  mixing  parameterizations  in  those  models. 

In  Figure  2.3.3  we  show  an  example  of  the  steps  in  calculating  Ri  for  one  station; 
station  76  during  the  fall  2001  survey  cruise  on  N.  B.  Palmer.  Well-defined  shear  layers 
were  present  at  the  base  of  the  mixed  layer  and  between  -150  m  and  -210  m  within  the 

main  pycnocline  (Figure  2.3.3c).  Values  of  Ri  within  the  latter  shear  layer  are  near  unity 

(Figure  2.3.30-  We  have  applied  these  methods  to  create  profiles  of  Ri  for  all  on-station 
data  obtained  during  the  2001  fall  and  winter  N.  B.  Palmer  survey  cruises.  Very  few 

profiles  show  Ri< 0.25  and  the  fraction  of  the  total  data  set  with  Ri<\  is  also  low, 
indicating  that  parameterized  occurrences  of  energetic  mixing  by  SD I  is  rare. 


Figure  2.3.3*  Example  showing  results  of  mixing  parameter  calculations  from  the  CTD  and  ADCP 
for  Station  76  from  the  fall  survey  cruise  (see  Fig,  2.2.1  for  location):  (a)  east  component  of  time- 
averaged  velocity  (u):  (h)  north  component  of  time-averaged  velocity  (v);  (c)  time-averaged  shear 
magnitude  <jf/z|>;  (d)  potential  density  o»;  (e)  Buoyancy  frequency,  (TV);  (f)  Richardson  number 

(vertical  dashed  lines  indicate  /?/=  1/4  and  /?/—!). 
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The  mean  diffusivity  that  would  be  predicted  by  the  turbulence 
parameterizations  in  general  circulation  models  can  be  estimated  by  applying  the  same 

algorithms  to  our  profiles  of  Ri .  We  report  this  calculation  for  just  one  such  scheme 
(PP-81).  The  PP-81  scheme  is  based  on  matching  numerical  model  predictions  with 
observations  of  the  Pacific  Ocean  Equatorial  Undercurrent,  and  so  may  have  little 
relevance  to  a  high-latitude  coastal-ocean  pycnocline.  Nevertheless,  Muench  et  al. 

(2002)  used  this  model  to  obtain  predictions  of  mixing  rates  that  were  consistent  with  the 
expected  dissipation  of  internal  tides  generated  along  the  South  Scotia  Ridge.  That  paper 
describes  in  detail  how  CTD  and  ADCP  data  can  be  used  with  PP-81  to  make  predictions 

of  Ky  ,  although  note  that  eq.  (2)  (Kv=f(Ri))  in  Muench  et  at.  (2002)  is  incorrectly  printed. 

Most  estimates  of  Ri,  in  both  the  sharp  pycnocline  and  the  more  weakly  stratified 
“deep'’  part  of  the  pycnocline,  are  sufficiently  large  that  PP-81  predicts  values  of  Kv 
which  are  close  to  the  asymptotic  lower  limit  of  lx  I  O'5  nr  s'1.  That  is,  there  is  no 
significant  difference  between  PP-81  and  using  a  constant  diffusivity  of  1  x  1 0  5  nr  s'1,  this 
constant  being  consistent  with  values  found  in  the  mid-latitude  main  thcrtnocline 
(Ledwell  et  al.,  1993;  Ruddick  et  al.,  1997).  The  heat  fluxes  F»  associated  with  this  value 
of  Kv  in  the  sharp  pycnocline  and  the  deep  pycnocline  are  about  1-2  W  m  ~  and 
<0.5  W  m  \  respectively.  A  more  recent,  “modified”  form  of  PP-81  (R.  Houghton, 

personal  communication,  2003)  has  a  lower  asymptotic  limit  for  Kv  at  high  Ri,  thus  the 
mean  value  of  F\\  may  be  <1  W  m  r 

There  are,  however,  a  few  areas  where  Ri  was  low  enough  to  indicate  mixing 

(e.g.,  the  2  depth  ranges  of  relatively  low  Ri  in  Figure  2.3.3f).  Shear  associated  with  the 
near-inertial  currents  at  process  site  3  (Figure  2.3. 1 )  of  the  winter  process  cruise  leads  to 

low  Ri  at  the  time  of  the  CTD  cast  shown  on  the  figure.  When  the  heat  flux  is  averaged 
vertically  over  the  entire  sharp  pycnocline  (the  region  in  which  N exceeds  3  eph),  the 
estimated  average  pycnocline  flux  at  this  time  is  ~5  W  m  \  The  shear  is  clearly  stronger 
later  in  the  time  series,  but  no  CTD  data  are  available  for  further  calculations  of  Kv  and 
Fh-  Intermittent  but  relatively  energetic  events  such  as  this  may  be  the  dominant 
contributors  to  mean  turbulent  diapycnal  fluxes. 
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Our  focus  in  this  chapter  is  on  mixing  in  the  main  pycnocline  as  a  mechanism  for 
transporting  the  scalar  properties  of  UCDW  into  the  surface  layer  through  the  intervening 
pycnocline.  Nevertheless,  as  we  have  previously  noted  (see  Figure  2.3.1 ),  the  variability 
of  mixed  layer  properties  is  dominated  by  the  seasonal  cycle  associated  with  surface 
buoyancy  fluxes.  Surface  stress  also  contributes  to  the  seasonal  cycling  of  the  upper 
ocean.  M.  McPhee  (personal  communication,  2004)  has  run  a  one-dimensional  local 
turbulence  closure  (LTC)  model  (McPhee,  1999)  for  the  winter  process  site  3  data  set 
reported  herein  (see  Figure  2.3.1).  His  model  successfully  reproduces  the  ice  drift 
velocity,  which  provides  some  support  for  other  predictions  of  the  LTC  model.  Strong 
upward  heat  flux  (~10  W  m 2  averaged  over  the  4  days)  is  predicted  in  the  center  of  the 
surface  layer  near  55  m  depth  from  ?«227.5  to  f*23l  .5  days.  This  average  is,  however, 
dominated  by  a  short  period  of  intense  mixing  from  /»230.4  to  /*230.8  days  in  which  the 
modeled  Fh  exceeds  50  W  m‘\  The  modeled  4-day  average  flux  of  ~  1 0  W  m':  for  the 
center  of  the  surface  mixing  layer  is  close  to  our  estimate  of  the  flux  through  the  sharp 
pycnocline  based  on  ADCP  and  CTD  data  at  the  time  of  the  single  CTD  profile  for  this 
site.  Wind  stress  clearly  contributes  to  the  transfer  of  heat  from  the  upper  pycnocline  to 
the  ocean  surface.  Importantly,  the  model  results  also  confirm  our  data-based  view  that 
rare  but  energetic  events  may  dominate  the  mean  mixing  rates  and  associated  fluxes  in 
the  upper  ocean.  That  is,  it  is  quite  possible  that  our  mean  flux  values  based  on  field  data 
are  lower  than  the  true  mean  values  because  of  inadequate  sampling  of  these  important 
events. 


2.3.2.Double  Diffusive  Convection 

Double-diffusive  convection  (DDC)  is  a  process  that  is  driven  by  the  different 
molecular  diffusivities  of  T and  S  in  a  stratified  fluid  where  both  properties  increase  with 
depth  (Turner,  1973).  Such  stratification  is  frequently  found  in  high-latitude  oceans 
where  a  cold,  fresh  surface  layer  formed  by  local  air/sea/ice  interactions  overlies  warmer 
and  saltier  water  that  has  been  advected  into  these  regions  from  lower  latitudes  (Padman 
and  Dillon,  1987;  Robertson  et  al.,  1995).  Where  DDC  is  active,  the  process  forms  a 
“thermoha  line  staircase",  a  vertical  sequence  of  quasi -homogeneous  layers  separated  by 
sharp  interfaces  in  which  vertical  gradients  of  T and  S  are  large.  This  feature  allows  us  to 
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identify,  from  high-resolution  CTD  or  scalar  microstructure  profiles,  where  DDC  is  likely 
to  be  a  significant  contributor  to  total  vertical  scalar  fluxes. 

Smith  and  Klinck  (2002)  estimated  that,  for  a  region  of  the  west  Antarctic 
Peninsula  continental  shelf  just  north  of  the  SO  GLOBEC  study  region,  ~5-10  W  in  ~  of 
upward  vertical  heat  flux  across  the  pycnocline  was  required  to  balance  the  total  heat 
budget  of  the  upper  ocean  in  their  one-dimensional  vertical  numerical  model  of  this 
region.  They  also  concluded  that  they  needed  a  higher  effective  pycnocline  diffusivity 
for  T  than  for  S,  which  is  consistent  with  DDC  being  the  dominant  source  of  mixing. 
These  authors  noted  the  presence  of  steppiness  in  many  of  the  CTD  profiles  from  their 
study  region:  the  vertical  scale  of  this  finestructurc  was  consistent  with  DDC  layer 
heights  that  would  be  required  to  reconcile  laboratory-based  DDC  flux  formulas  with  the 
heat  and  salt  fluxes  determined  by  their  model -based  analyses. 

Our  inspection  of  the  “raw”  (24  Hz)  SO  GLOBEC  CTD  profiles  from  the  winter 
2001  survey  cruise  identified  a  few  profiles  containing  finestructurc  that  was  reminiscent 
of  thcrmohaline  steps,  suggesting  that  DDC  may  occasionally  be  contribute  to  vertical 
fluxes  in  the  study  area  and  that  further  study  of  DDC  is  warranted.  However,  in  general 
the  observed  “interfaces”  are  much  thicker  than  we  would  expect  for  DDC  (—1-10  cm: 
Padman  and  Dillon  (1987)).  We  therefore  believe  that  the  sleppy  hydrographic  structure 
is  generally  associated  with  non-DDC  processes  such  as  intrusive  activity  near  the  coast 
and  ice  fronts,  and  the  fi  restructure  remnants  of  shear-driven  mixing. 

In  keeping,  however,  with  our  goal  of  investigating  fluxes  that  would  arise  from 
the  application  of  1-D  turbulence  parameterizations  such  as  Large  et  al.  (1994),  we  apply 
models  of  DDC  fluxes  to  the  CTD  data  collected  in  2001  during  the  winter  survey  cruise 
on  N.  B.  Palmer.  The  fundamental  parameter  describing  DDC  activity  is  the  density  ratio 
RP (Turner,  1973),  given  by: 


_  pdS  /  dz 
adO/dz' 


(2.4) 


where  p  (psu1)  is  the  haline  contraction  coefficient  and  a  (°C  ')  is  the  thermal 
expansion  coefficient.  To  create  profiles  of  Rp,  7\z)  and  S(z)  were  first  filtered  over  a 
30-m  vertical  scale,  which  was  chosen  to  be  much  larger  than  likely  DDC  layer  heights  of 
—  1  - 1 0  m.  The  mean  gradients  of  the  filtered  T  and  S  were  then  calculated  over  the  high- 
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gradient  core  of  the  pycnocline  separating  the  surface  mixed  layer  from  the  depth  of  the 
UCDW  temperature  maximum.  The  DDC  instability  cannot  occur  below  the  depth  of 
maximum  temperature  due  to  the  change  in  sign  of  dOfdz  .  Values  of  Rp  slightly  above  I 
indicate  strong  double  diffusion.  As  Rp  increases  above  2,  DDC  fluxes  become  quite 
small.  Estimates  of  DDC  heat  fluxes  in  polar  oceans  range  from  <0. 1  W  m 2  in  the 
Canada  Basin  in  the  Arctic  (Rp>4;  Padman  and  Dillon,  1987)  to  ~20  W  m 2  in  the  central 
Weddell  Sea  (Rp«l;  Mueneh  et  al.,  1990). 

We  calculated  Rf,  over  the  winter  survey  cruise  to  estimate  the  possible  strength 
of  DDC  processes  on  the  WAP  shelf.  We  concentrated  on  this  survey  because  the 
seasonal  upper-ocean  stratification  that  is  present  throughout  much  of  the  fall  survey 
cruise  complicates  automated  detection  of  the  core  of  the  pycnocline  below  the  Winter 
Water,  the  subsurface  temperature  minimum  layer  representing  the  remnant  of  the 
previous  winter's  near-freezing  surface  mixed  layer.  Furthermore,  CTD  profiles  taken 
while  the  vessel  is  becalmed  in  ice  are  of  higher  quality  than  the  wave- influenced  “open 
ocean"  profiles  obtained  during  the  fall  survey  cruise.  Some  steppiness  is  seen  in  profiles 
from  the  winter  survey  cruise  when  Rp  is  low,  ~2.  This  relatively  low  value  of  Rp 
indicates  the  potential  for  DDC  to  occur;  however,  inspection  of  the  raw  (24  Hz)  CTD 
profiles  of  T  and  C  confirm  that  sharp  steps  are  rare.  Most  of  the  profiles  in  the  WAP 
region  have  higher  Rr  values  (typically  -3.6)  in  the  pycnocline,  indicating  only  weak 
DDC. 

Several  models  have  been  developed  to  estimate  DDC  heat  fluxes.  Robertson  ct 
al.  (1995)  provide  an  overview  and  comparison  (their  Fig.  7)  of  these  models.  In  the 
present  study,  we  use  the  Kelley  (1990)  heat  flux  parameterization,  defined  as 

Fh-k  ~ 0.0032pi(£*/jtr”1  exp[4.8/?^'72](g/r,2L>_1  )I/3(£?A0)4/J ,  (2.5) 

where  p„  is  the  mean  density,  cp  is  the  specific  heat  capacity  of  seawater  at  constant 
pressure  (  -4.0  x  10?  J  kg  '  K  '),  k,  is  the  molecular  diffusivity  for  heat,  o  is  the  kinematic 
viscosity  of  seawater,  and  A0  is  the  thermal  step  between  adjacent  well-mixed  layers. 
The  thermal  step  is  estimated  by  first  obtaining  the  DDC  layer  height,  H ,  using  a 
relationship  based  on  large-scale  stratification  parameters  proposed  by  Kelley  (1984); 

H  =[0.25jcl0,«;iw<l(i?p-l))M5(ir((Ar1»M.  (2.6) 
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The  thermal  step  is  then  obtained  using 

AO  =  (dO/dz)H ,  (2.7) 

where  the  temperature  gradient  (dO/dz)  is  the  average  gradient  over  the  depth  range  in 
which  evidence  of  DDC  is  sought. 

Heat  fluxes  were  estimated  by  applying  (5)-(7)  to  all  stations  from  the  winter 
2001  survey  cruise.  The  area-averaged  heat  flux  due  to  DDC  was  ~0. 2-0.4  W  m  ‘  over 
the  central  and  inner  continental  shelf,  and  -0.4-0. 9  W  m  ‘  over  the  continental  slope 
where  the  average  value  of  /?p  was  smaller  than  for  the  central  shelf.  The  error  bars 
associated  with  these  values  are  uncertain.  Another  DDC  heat  flux  model  developed  by 
Marmorino  and  Caldwell  (1976)  predicts  values  up  to  a  factor  of  2  greater  than  Fh-k  for 
Rvalues  between  1.4  and  1.5  (Robertson  et  al.,  1995).  However,  as  /(^increases,  the 
differences  between  the  estimates  from  the  various  analytical  models  decrease  (see  also. 
Fig.  4b  in  Large  et  al.,  1994).  Note,  however,  that  all  models  of  DDC  assume  that  the 
hydrography  has  organized  into  a  thermohaline  staircase,  in  which  a  continuous  upward 
flux  of  heat  is  maintained  by  molecular  diffusion  through  extremely  sharp  interfaces 
(“sheets”)  -1-10  cm  thick  and  buoyancy-driven  convection  in  the  adjacent  layers 
-1-10  m  thick  (see  Padman  and  Dillon,  1987).  When  these  sheets  and  layers  are  not 
present,  as  in  the  present  data  set,  expected  heat  fluxes  will  be  less  than  the  models 
predict. 

Assuming,  however,  that  these  fluxes  do  at  least  represent  some  potential  level  of 
DDC  activity,  we  can  estimate  the  thermal  diffusivities  associated  with  these  heat  fluxes 
from 


k!  = 


F„ 


P„cp(dO/dz) 


(m;  s'1). 


(2.8) 


where  ( dQ/dz )  is  again  averaged  over  the  pycnocline.  Using  this  method,  the  average 
DDC  thermal  diffusivity  in  the  main  pycnocline  on  the  continental  slope  is 
-0.5  x  10  s  m:  s'1,  and  -0.3  x  1 0'5  m2  s'1  over  the  central  and  inner  continental  shelf. 
The  diffusivity  for  salt,  Ff ,  is  lower  than  K[  for  DDC  (see  Fig.  2  in  Kelley, 


1984).  The  ratio  of  the  two  diffusivities  is  given  by 


K>IKl  =  R,.Rp 


(2.9) 
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(2.10) 


{Kelley,  1984),  where  the  buoyancy  flux  ratio  R/.  is  given  by 

Rp+\A(Rp-])m 
F  l  +  14(/f/, -1),/2  ’ 

(Kelley,  1990).  Based  on  (9)  and  (10),  AT;S  is  -0.5  x  10'<>  m2  s'1  over  the  continental  slope 

and  -0. 1  x  10  m“  s*  over  the  central  and  inner  continental  shelf.  The  salt  fluxes 
associated  with  these  diffusivities  are  -0.002  mg  tn 2  s'1  over  the  continental  slope  and 
-0.001  mg  m'2  s'1  over  the  shelf. 

It  remains  an  open  question  whether  DDC  fluxes  in  the  ocean  can  be  correctly 
parameterized  by  formulas  developed  in  laboratory  studies.  For  example,  Padman  { 1 994) 
suggested  that  the  addition  of  mean  shear  to  thermohaline  steps  can  increase  the  vertical 
momentum  and  scalar  fluxes  relative  to  unsheared  steps.  A  fundamental  assumption  in 
the  Padman  (1994)  model  was  that  the  applied  shear  was  “relatively”  weak  (incapable  of 
causing  mixing  in  the  absence  of  the  finestructure  established  by  DDC),  and  the  sharp 
DDC  interfaces  were,  therefore,  not  significantly  disrupted  by  the  applied  shear.  We 
reiterate  that  very  few  cases  of  sharp  interfaces  separating  truly  homogeneous  layers  were 
found  in  the  SO  GLOBEC  data,  which  implies  that  DDC  was  not  active  at  the  time  the 
observations  were  made.  There  are  several  plausible  explanations  for  the  observed 
steppy  fine-structure  in  some  CTD  profiles,  however  it  is  possible  that  some  of  this 
finestructure  points  to  prior  double-diffusive  activity  that  has  since  been  partially  erased 
by  more  recent  shear-driven  turbulence. 

2.4.  Conclusions 

We  use  CTD  and  ADCP  data  from  the  first  year  (2001)  of  Southern  Ocean 
GLOBEC  field  work  to  study  diapycnal  mixing  processes  in  the  pycnocline  separating 
UCDW  from  the  surface  layer  over  the  western  Antarctic  Peninsula  (WAP)  continental 
shelf  near  Marguerite  Bay.  Based  on  a  recent  box-model  study  of  a  region  of  the  WAP 
shelf  just  north  of  the  SO  GLOBEC  study  area  (Smith  and  Klinck,  2002),  our  a  priori 
expectation  was  that  double-diffusive  convection  (DDC)  would  be  the  dominant  process 
driving  the  upward  flux  of  heat,  salt  and  nutrients  from  the  UCDW.  with  an  additional 
contribution  from  shear-driven  instability  (SDI). 
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Detailed  examination  of  high-resolution  CTD  profiles  from  the  winter  2001 
survey  cruise  on  N.  B.  Palmer,  combined  with  the  application  of  common 
paramcterizations  of  DDC  fluxes  to  the  observed  hydrographic  profiles,  suggest  that 
DDC  is  either  not  active  or  is  very  weak  during  this  cruise.  The  maximum  likely  DDC 
heat  flux  is  <1  W  m'2,  compared  with  Smith  and  Klinck’s  estimate  of  order  10  W  m  \ 

We  attribute  most  of  the  velocity  shear  observed  in  the  pycnocline  with  hull- 
mounted  ADCP  to  wind  forcing.  With  the  exception  of  small  regions  near  the  shelf 
break,  the  tidal  currents  predicted  from  an  Antarctic  tide  model  (CATS02.01 )  are  an 
order  of  magnitude  weaker  than  typical  wind-forced,  near-inertial  oscillations  of  the 
upper  ocean.  A  simple  parameterization  of  vertical  mixing  due  to  SDI,  following 
Pacanowski  and  Philander  (1981),  suggests  that  the  mean  vertical  diffusivity  (Kv)  in  the 
pycnocline  in  winter  2001  is  ~1  x  10  ^  nr  s  '  with  a  corresponding  upward  heat  flux  from 
the  UCDW  of  ~1  W  m 2.  We  caution,  however,  that  this  parameterization  was  developed 
for  oceanic  conditions  (the  mean-sheared  Equatorial  Undercurrent)  that  are  quite  different 
from  the  WAP  region. 

The  pycnocline  fluxes  that  we  estimate  are  two  orders  of  magnitude  less  than 
those  required  to  explain  the  loss  of  heat  from  the  surface  mixed  layer  between  the  fall 
and  winter  survey  cruises:  as  expected,  surface  buoyancy  fluxes  are  the  primary  causes  of 
the  strong  seasonal  cycle  in  the  upper  ocean.  The  destabilizing  buoyancy  flux  in  fall  and 
winter,  as  a  source  of  entrainment  of  the  pycnocline  at  the  base  of  the  mixed  layer,  is 
augmented  by  the  turbulent  stress  penetrating  the  ice  pack. 

We  conclude  that  the  downward  flux  of  momentum  originating  in  the  wind  stress 
field  is  the  dominant  process  driving  diapycnal  fluxes  in  the  pycnocline  through 
downward  propagation  of  near-inertial  internal  waves,  and  is  also  a  critical  component  of 
upper-ocean  seasonal  variability.  There  are  several  practical  consequences  to  this 
observation.  First,  since  the  actual  stress  transmitted  to  the  upper  ocean  through  the  sea 
ice  (when  present)  depends  on  the  roughness,  thickness  and  mobility  of  the  sea  ice, 
mixing  in  the  pycnocline  cannot  be  evaluated  in  isolation  from  either  accurately 
prescribed  or  modeled  sea  ice  conditions.  Second,  since  the  wind  stress  in  this  region  is 
extremely  variable  (U.S.  SO  GLOBEC,  2001a,b,c),  we  risk  undersampling  the  rare  but 
energetic  events  that  might  dominate  mean  mixing  rates,  and  so  may  significantly 
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underestimate  mean  pycnoclinc  diffusivity.  Third,  generating  realistic  near-inertial 
behavior  in  a  model  of  the  upper  ocean  requires  atmospheric  model  forcing  fields  with 
high  resolution  in  both  time  and  space.  The  output  from  global  atmospheric  models  is 
typically  every  6-12  h  and  on  a  coarse  spatial  grid,  thus  a  new  high -resolution  regional 
atmospheric  model  will  be  required  before  WAP  shelf  circulation  models  can  hope  to 
capture  the  physics  of  mixing. 

We  hypothesize  that,  given  the  importance  of  wind  forcing,  coastal  upwelling 
could  be  a  significant  factor  contributing  to  the  upward  flux  of  UCDW  properties. 
Upwelling  may  raise  UCDW  into  the  depth  range  affected  directly  by  surface  stress. 
Given  the  short  Rossby  radius  (~5  km)  in  high-latitude  coastal  seas,  sampling  very  close 
to  the  coast  under  upwelling- favorable  wind  conditions  would  be  required  to  test  this 
hypothesis. 

Finally,  we  note  that  the  study  period  covers  only  the  fall  and  winter  of  one  year, 
200 1 ,  whereas  we  expect  that  both  SD1  and  DDC  could  be  sensitive  to  interannual 
variability  of  stratification.  The  Amundsen-Bellingshausen  Sea  sector  of  Antarctica 
experiences  significant  El  Nino-Southern  Oscillation  (ENSO)  atmospheric  and  sea-ice 
variability  (Kwok  and  Comiso,  2002),  while  the  Antarctic  Circumpolar  Wave  (ACW) 
provides  an  additional  multi-year  signal  (White  and  Peterson,  1996).  It  is  possible  that 
quite  different  annual-mean  fluxes  in  the  pycnocline  occur  during  appropriate  phases  of 
the  ENSO  and  ACW  cycles. 
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3.  Chapter  3 


Measurements  of  sea  ice  properties  using  a  moored 
upward-looking  acoustic  Doppler  current  profiler 
(ADCP). 

3.1.  Introduction 

Upward-looking  Acoustic  Doppler  Current  Profilers  (ADCPs)  have  been 
successfully  deployed  in  the  world  oceans  for  many  years  on  bottom-anchored  subsurface 
moorings.  The  primary  purpose  for  the  ADCP  has  been  to  measure  vertical  profiles  of 
ocean  currents.  When  an  upward-looking  ADCP  is  deployed  in  a  location  that 
experiences  seasonal  sea  ice  cover,  it  can  also  be  used  to  collect  ice  data  (Belliveau  et  al 
1990;  Strass,  1998).  This  has  become  more  important  as  the  effort  to  observe  polar  seas 
has  increased  due  to  their  significant  role  in  the  climate  system.  The  remoteness  of  the 
polar  seas  demands  the  use  of  every  piece  of  information  obtained. 

The  presence  of  sea  ice  has  a  strong  influence  on  the  exchange  of  heat  and 
momentum  between  the  atmosphere  and  the  ocean.  The  lateral  motion  of  sea  ice  moves 
both  fresh  water  and  heat.  Furthermore,  the  presence  of  seasonal  sea  ice  is  believed  to 
play  a  vital  role  in  krill  ecosystems  throughout  Antarctica.  Combined  in  situ  observations 
of  sea  ice  coverage,  velocity  and  draft  provide  data  are  of  great  benefit  to  many 
investigations,  particularly  when  measured  at  the  same  time  and  location. 

Ice  draft  has  been  successfully  measured  with  a  dedicated  upward-looking  sonar 
(ULS)  (Hudson,  1990;  Strass,  1998).  An  ULS  typically  samples  at  high  frequency  (e.g., 
0.5  Hz)  using  one  narrow  vertical  beam.  This  single  beam  avoids  errors  due  to  lateral 
scattering  between  multiple  slanted  beams  that  can  occur  with  the  four-beam  ADCP,  Ice 
drift  velocities  can  also  be  measured  using  ADCPs  in  bottom-track  mode  (Belliveau  et  al, 
1990).  These  bottom-track  pings  require  additional  power  which  can  be  problematic  in 
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long-term  moored  deployments  and  may  be  unavailable  in  post- processing  of  archived 
ADCP  data.  Visbeck  and  Fischer  (1995)  estimated  the  presence  of  ice  and  drift  velocity 
using  a  narrow-band  ADCP  in  water-track  mode,  the  mode  normally  used  for  water 
velocity  profiling,  but  did  not  estimate  ice  draft.  The  combination  of  ULS  and  ADCP 
provides  the  best  measurements  of  ice  draft  and  motion  and  allows  proper  spatial 
averaging  of  ice  properties  to  obtain  statistical  descriptions  of  ice  (Mel ling  et  a/,  1995). 
This  configuration  still  has  problems  estimating  ice  draft  due  to  errors  in  the  estimation  of 
range-averaged  sound  speed  since  one  cannot  make  in-situ  measurements  between  the 
ULS  and  the  sea  surface  (the  sensors  would  block  the  single  vertical  acoustic  beam).  The 
methodology  presented  here  provides  an  improvement  that  could  be  used  even  in  the 
ULS/ADCP  configuration. 

In  sections  3.2  and  3.3,  we  describe  a  method  for  processing  water- track  mode 
ADCP  data  from  a  standard  subsurface  mooring  (as  in  Figure  3. 1 . 1 )  to  estimate  ice 
coverage  and  draft  that  can  be  applied  to  both  broad-  and  narrow-band  ADCP  data.  This 
method  provides  estimates  of  sea  ice  draft  in  the  absence  of  ULS  data  by  taking 
advantage  of  auxiliary  data  from  other  instruments  on  the  mooring.  We  estimate  the 
uncertainty  in  the  ice  draft  estimate,  and  compare  this  with  the  uncertainty  of  an  ULS 
estimate.  The  information  provided  by  this  method,  while  by  no  means  a  substitute  for 
actual  ULS  data,  can  improve  the  interpretation  of  the  ocean  current  data  recorded  by  the 
ADCP.  This  method  can  also  be  used  with  archived  ADCP  data  to  provide  historical 
estimates  of  ice  draft. 
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Figure  3.1.1  Schematic  of  mooring  and  names  of  variables  used  in  ice  draft  calculations. 

The  following  information  is  needed  to  estimate  sea  ice  coverage,  draft  and 
velocity: 

1 .  The  data  from  an  upward-looking  ADCP  on  a  subsurface  mooring. 

2.  The  data  from  a  high-precision  pressure  gauge  on  the  mooring,  preferably 
close  to  the  ADCP. 

3.  Atmospheric  pressure  and  some  estimate  of  the  lime- variable  sound  speed 
profile  above  the  ADCP.  In  the  case  shown  here,  the  atmospheric  pressure  was  measured 
nearby  and  the  sound  speed  profile  is  estimated  based  on  moored  temperature  and  salinity 
measurements  made  near  the  ADCP  in  conjunction  with  shipboard  CTD  profiles  made  in 
the  area. 

The  ice  drift  velocity  can  simply  be  taken  as  the  velocity  in  the  bin  just  above  the 
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surface  bin  (maximum  backscatter)  during  times  of  ice  coverage  (Visbeck  and  Fischer, 
1995)  using  standard  water-track  pings. 

Data  from  a  mooring  in  Marguerite  Bay  on  the  western  Antarctic  Peninsula  shelf 
are  analyzed  to  demonstrate  that  the  proposed  method  can  provide  useful  estimates  of  ice 
draft.  These  data  are  described  in  greater  detail  in  the  appendix.  The  final  two  sections 
of  this  chapter  make  recommendations  for  future  deployments  of  ADCPs  in  ice  covered 
regions  and  present  conclusions  of  our  study. 

3.2.  Ice  detection 

The  detection  of  ice  has  been  performed  using  ADCP  bottom-track  pings 
(Belliveau  el  al.  1990),  upward -loo king  sonars  (ULS)  (Strass,  1998),  and  narrow-band 
ADCP  water-track  pings  (Visbeck  and  Fischer,  1995).  The  method  described  here 
essentially  duplicates  the  method  for  a  narrow-band  ADCP  in  water-track  mode,  with 
some  noted  differences  for  a  broad-band  ADCP. 

We  employ  a  set  of  criteria  which  discriminate,  from  time  averages  of  the  ADCP 
data,  in  a  binary  way  the  presence  of  ice.  First,  the  ADCP  bin  that  samples  the  sea 
surface  or  bottom  of  the  sea  ice  is  identified  as  the  bin  above  the  one  with  maximum 
backscatter  intensity.  The  bin  above  maximum  backscatter  is  chosen  to  ensure  ice-only 
data  due  to  overlapping  cell  information  as  a  result  of  the  range  gating  in  an  ADCP  (RDI- 
Primer,  1996).  While  it  is  possible  to  use  a  weighted  average  of  bins,  we  chose  the 
simplest  option  of  using  only  one  bin. 

Possible  ice  detection  criteria  are  vertical  and  error  velocity  windowed 
variances,  surface  backscatter  intensity,  the  horizontal  surface  speed,  and  the  signal 
correlation.  The  signal  correlation  is  the  only  additional  piece  of  information  provided  by 
the  broad-band  ADCP  which  is  not  provided  by  a  narrow-band  ADCP.  These  properties 
arc  all  greatly  influenced  by  the  presence  of  sea  ice  (Visbeck  and  Fischer,  1995).  A 
successful  criterion  for  the  presence  of  sea  ice  is  one  where  the  variable  has  a  strongly 
bimodal  distribution  with  the  two  maxima  in  the  probability  distribution  function 
significantly  separated  and  identifiable  with  times  when  the  location  is  ice-covered  or  ice- 
free.  One  then  defines  a  “cut-off’  value  that  delineates  the  boundary  between  the  two 
regions. 


45 


The  top  pane!  of  Figure  3.2. 1  shows  a  time  series  of  the  high-pass  boxcar  filtered 
error  velocity  with  a  filter  period  of  24  hrs  and  clearly  reveals  approximate  times  of  ice. 
This  is  used  to  make  a  first  approximation  of  the  presence  of  ice  (marked  in  black)  based 
on  inspection  of  the  time  series  and  identification  of  periods  of  low  variance  by  eye.  (In 
this  example,  we  based  the  first  approximation  of  ice  cover  on  the  filtered  error  velocity, 
but  one  could  choose  any  one  of  the  criteria,  such  as  vertical  velocity,  and  get  nearly 
identical  results).  Next,  cumulative  histograms  are  made  of  each  of  the  criteria  to  find  a 
suitable  cutoff  value  (Figure  3.2. 1 ).  Criteria  which  show  a  clear  separation  between  ice 
and  no  ice  are,  in  our  case,  the  windowed  variance  of  the  error  velocity,  the  squared 
windowed  variance  of  vertical  velocity,  and  the  surface  speed.  We  squared  the 
windowed  variance  of  vertical  velocity  in  order  to  expand  the  scale  of  the  data. 
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Figure  3.2.1  (a)  Time  series  of  surface  error  velocity  (It)'1  cm/s),  with  black  showing  initial  selection  of 
ice  periods.  The  other  panels  show  cumulative  histograms,  with  the  percentage  on  the  y  axis  and  the 
value  of  the  data  on  the  x  axis  of  (b)  variance  of  surface  error  velocity  (cmV),  (c)  squared  variance 
of  surface  vertical  velocity  (ciriVs4),  (d)  surface  speed  (cni/s),  (c)  surface  backscatter  intensity  ( Rl)l 
relative  backscatter  units,  or  BSU),  and  (f)  surface  signal  correlation  (no  units).  In  the  histograms, 
the  black  and  gray  lines  contain  data  from  the  time  periods  shown  in  black  and  grey  in  panel  (a) 
respectively,  and  the  dashed  lines  contain  data  from  the  entire  dataset.  The  vertical  black  lines  in 
panels  (b)  -  (d)  show  the  threshold  values  as  described  in  the  text. 


The  cutoff  values  for  the  presence  of  ice  were  first  determined  by  examining  these 
histograms.  As  an  ad  hoc  rule,  we  have  taken  the  cutoff  values  to  be  when  the 
cumulative  histograms  for  the  ice  period  reaches  90%  and  the  open  water  period  is  less 
than  10%.  These  values  were  chosen  using  a  subjective  trial  and  error  process.  They  not 
only  work  best  at  identification  of  the  best  criteria  for  discerning  ice  in  the  record,  but 
provide  the  best  cutoff  values  for  these  criteria  when  compared  with  satellite-based 
estimates  of  ice  coverage,  as  described  below. 

Then  a  time  series  of  ice  presence  is  made  by  tagging  each  half  hourly  data  point 
as  ‘ice’  or  'no  ice’.  One  can  choose  to  simply  use  one  of  the  criteria,  or  combine  them  by 
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taking  the  average  or  cross-correlation  of  two  or  more  indicators.  Next,  a  daily  average 
ice  concentration  in  the  area  of  the  mooring  is  calculated  as  the  percentage  of 'ice' 
measurements  in  one  day  (Visbeck  and  Fischer,  1995).  This  assumes  that  the  ice  is 
moving  overhead  in  a  statistically  random  way  on  daily  time  scales. 

Finally,  the  time  series  of  ice  concentration  can  be  compared  with  remotely- 
sensed  ice  estimates  and  cutoff  values  for  criteria  for  ice  coverage  can  be  fine  tuned  if 
necessary.  One  might  need  to  identify  and  possibly  omit  times  of  calm  winds  which  can 
produce  low  vertical  velocity  variance  and  thus  false  ice  detection  (Visbeck  and  Fischer, 
1995).  However,  in  the  example  shown  here,  surface  wind  data  collected  nearby  had  no 
periods  of  time  where  winds  were  light  enough  for  a  long  enough  period  of  time  (less 
than  5  ms'1  for  more  than  2  days)  to  warrant  this  extra  step. 

The  ice  concentrations  derived  using  the  initial  choice  of  cut-off  values  from 
Figure  3.2.1  are  indicated  by  the  black  line  in  Figure  3.2.2.  In  this  case,  we  applied  the 
cutoff  criteria  of  windowed  variance  of  error  velocity  (0.0006  cm2  s'2),  the  squared 
windowed  variance  of  vertical  velocity  (0.4  cm4s4),  and  the  surface  speed  (26.5  cm  s’1), 
and  then  averaged  the  ice  concentration  derived  from  these  three.  The  two  SSMI  satellite 
estimates  in  Figure  2.2  employ  different  algorithms  and  show  the  range  in  estimates  of 
remotely-sensed  sea  ice  coverage  by  the  shaded  gray  area.  The  ADCP  estimate  derived 
from  the  initial  selection  falls  within  the  range  of  the  satellite  estimates  for  most  of  the 
time  series  without  further  correction.  That  is,  the  satellite  data  was  not  used  to  correct 
the  cutoff  criteria. 
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Figure  3.2.2  Estimated  sea  ice  concentrations  at  the  mooring  site  using  the  A  DCF  data  and  the 
methods  described  in  the  text  (black  line),  and  estimated  by  two  different  algorithms  applied  to  SS1VII 
satellite  data  whose  range  are  show  n  by  the  gray  shaded  area. 

The  greater  variability  in  the  ADCP  estimate  can  be  attributed  to  the  spatial 
average  for  the  moored  ADCP  versus  that  of  the  SSMI  satellite.  The  nominal  pixel  size 
of  the  satellite  products  is  12.5  km  and  25  km  (see  appendix  for  details),  making  156.25 
km2  and  625  krrf  footprints.  The  total  area  of  the  footprints  of  the  four  ADCP  beams  is 

4^[o.5/?[sin{0  +  <pi2)  -  sin(6>  -  (p! 2)Jj2 ,  where  0  is  the  beam  angle,  ip  is  the  beam  spread 
and  hti  is  the  nominal  depth  of  the  ADCP  (Figure  3.1.1).  In  our  example,  the  300  kHz 

broad-band  ADCP  had  a  beam  angle  of  20°,  a  beam  spread  of  2.2°  at  -3dB  (RDI-Primer, 
1996),  and  a  nominal  depth  of  108  m,  giving  an  instantaneous  footprint  of  48  m\  One 
must  take  into  account  ice  motion  when  making  averages  of  moored  ice  observations 
using  the  width  of  the  beam  times  the  distance  the  ice  has  traveled  for  the  footprint. 

Here,  the  total  width  of  the  four  beams,  4[/t[sin((?+  (pi  2)-  sin{0-(p/2)]],  is  16  m  and  the 
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average  distance  traveled  by  the  ice  in  a  day  is  I .  I  km,  with  95%  of  the  values  less  than 
28  km.  This  gives  an  average  footprint  of  1 7.600  m\  or  four  orders  of  magnitude  smaller 
than  the  satellite  footprints.  For  a  thorough  treatment  of  the  subject  of  spatial  and 
temporal  averaging  of  ice  data  from  a  moored  platform,  the  reader  is  referred  to  Melling 
etal(  1995). 

3.3.  Ice  draft 

Ice  draft,  hjie  in  Figure  3.1 . 1 ,  is  defined  as  the  depth  of  ice  below  waterline.  The 
estimate  of  ice  draft  is  made  using  ocean  pressure,  atmospheric  pressure,  sound  speed 
derived  from  temperature  and  salinity,  and  backscatter  and  tilt  data  from  a  moored 
upward-looking  ADCP.  The  ice  draft  is  calculated  as  the  difference  between  the  depth  of 
the  ADCP  and  the  acoustically  measured  distance  to  the  underside  of  the  ice: 

K,  =  K-K  (3-D 

Using  hydrostatics,  the  ADCP  depth,  or  distance  from  the  waterline  to  the  ADCP, 
can  be  calculated  as: 

p  -P 

/?„  =  ^  ~h,  (3.2) 

Pg 

with  the  variables  defined  in  Figure  3.1.1.  The  pressure  P„  is  measured  anywhere  near 
the  ADCP  and  h2  is  determined  by  the  mooring  design.  The  time  varying  average 
density  between  the  ADCP  and  the  sea  surface,  p ,  is  calculated  from  CTD  casts  taken  at 

deployment  and  recovery  and  moored  temperature  and  salinity  time  series  data. 
Substituting  equation  (3.2)  into  the  ice  draft  expression,  equation  (3. 1 ),  we  have  an 
expression  for  the  ice  draft: 


P  -P 

h  -  -1L. _ JUL.  _  /,  _  /, 

n,ce  ~  ~  nl  nl 

Pg 


(3.3) 


The  cumulative  uncertainty  in  the  first  two  terms  is  less  than  1  cm.  The  acoustic  range 
estimate  from  the  ADCP  to  the  underside  of  the  ice,  h\  requires  careful  consideration  and 
a  series  of  specific  corrections. 
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The  initial  estimate  of  /»/  is  obtained  for  each  of  the  four  beams  using  the 
backscatter  profile.  To  obtain  an  estimate  of  hi  at  a  resolution  greater  than  the  bin 
intervals,  the  centered  first-difference  of  the  observed  backscatter  profile  dBS  at  zf ,  the 

midpoint  between  bins  i  and  /+ 1 ,  is  modeled  using  a  slightly  modified  first-difference 
of  a  standard  normal  distribution: 


dBS{zi)  =  A 
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This  functional  form  includes  a  constant  multiplier  A  and  an  added  constant  C  to 
improve  the  fit  with  the  data  and  have  the  same  units  as  dBS .  a ,  ft  and  z,  have  units  of 


meters.  The  four  parameters  A ,  C ,  a  and  //  are  obtained  by  a  non-linear  least-squares 
fit  to  the  observed  dBS  in  3  bins  above  and  below  the  bin  with  the  maximum 
backscatter,  for  7  total  bins,  or  6  differences.  This  gives  two  degrees  of  freedom  for  this 

fit.  The  zero-crossing  of  dBS ,  defined  as  dBS{z„)  =  0 ,  is  the  peak  of  the  fitted 

backscatter  profile,  which  is  taken  as  the  ice- water  or  air- water  interface  (Visbeck  and 
Fischer,  1995).  We  believe  that  this  function  improves  upon  the  second-order  polynomial 
fit  used  by  Visbeck  and  Fischer  (1995)  since  it  gives  a  range  to  the  surface  which 
achieves  a  better  match  to  the  hydrostatic- based  range  ho  during  times  of  open  water. 

The  initial  attempt  to  estimate  hi  using  the  parabolic  fit  did  not  track  variations  of 
pressure  at  tidal  frequencies,  while  equation  (3.4)  did  accomplish  this  (section  3.6). 
Furthermore,  there  is  no  acoustic  process  which  favors  or  suggests  one  empirical  fitting 
function  over  another. 

In  the  example  here,  a  20°-beam  angle,  a  4.4°-beani  aperture  and  a  nominal 
instrument  depth  of  108  rn  yields  a  minimal  width  of  the  surface  energy  peak  of  2.9  m 
(equation  1  in  Visbeck  and  Fisher,  1995),  which  is  resolved  by  the  2~m  bin  averaging. 
This  confirms  that  the  surface  energy  peak  affects  more  than  one  of  the  ADCP  bins. 

A  set  of  corrections  to  this  raw  estimate  must  be  made.  These  corrections  take 
into  account  changes  which  would  not  greatly  affect  the  velocity  profile  and  arc  therefore 
not  included  in  standard  ADCP  processing,  but  are  of  the  order  of  the  signal  when  ice 
draft  is  estimated.  Since  they  are  multiplicative,  they  become  more  important  as  the 
depth  of  the  ADCP  increases.  These  corrections  for  sound  speed  and  instrument  tilt  are 
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described  in  detail  in  separate  sections  below.  As  a  final  step,  bad  data  are  removed  and 
a  footprint  error  correction  is  made. 

Finally,  a  best  estimate  for  vertical  range  za  is  made  from  the  4  beams  by 
averaging  the  2  most  similar  corrected  estimates  4(2,,,  +  z()2).  Local  atmospheric 

pressure  is  used  to  estimate  the  depth  of  the  ADCP  (equation  3.3).  In  our  case,  Pair  was 
measured  at  a  nearby  automatic  weather  station.  If  this  data  were  not  available,  one 
could  use  meteorological  reanalysis  products  to  estimate  Pair-  However,  in  a  location 
with  little  actual  data  these  reanalysis  products  are  much  less  reliable  than  actual 
measurements  (e.g.  Marshall,  2002). 


3. 3.1.  Sound  speed  correction 

A  site-specific  method  for  estimating  the  sound  speed  profile  is  required  to 
accurately  convert  the  round-trip  travel  time  to  range  (Strass,  1998).  Standard  RDI  ADCP 
processing  uses  a  depth- indepen dent  sound  speed  profile  computed  from  the  ADCP 
temperature  measurement  and  a  constant  salinity  specified  by  the  user,  since  sound  speed 
variation  is  much  more  sensitive  to  changes  in  temperature  than  salinity  (RDI-Primer, 
1996).  The  multiplicative  factor  to  correct  for  the  sound  speed  profile  is 


r  „  „ 

AIK'P 

c 

^  real 


(3.5) 


where  Cadcp  is  the  sound  speed  used  by  the  ADCP,  and  Creai  is  the  harmonic  mean  of  the 
true  sound  speed  profile. 

In  most  cases,  one  does  not  have  a  full  time  series  of  sound  speed  profile.  In  fact, 
deep  keeled  icebergs  often  occur  in  regions  of  sea  ice,  jeopardizing  any  attempts  at  upper 
water  column  density  measurements.  Therefore,  one  must  use  other  information  and  a 
model  of  the  evolving  mixed  layer  to  estimate  the  sound  speed  profile  (Strass,  1998).  We 
use  the  UNESCO  formulas  for  calculating  seawater  properties,  such  as  sound  speed 
(Fofonoff  and  Millard,  1983)  in  our  example. 

The  first  correction  uses  salinity  measured  at  a  sensor  very  close  (in  this  case,  10 
m  above)  to  the  ADCP  to  replace  the  constant  salinity  preprogrammed  into  the  ADCP. 
Another  simple  correction  takes  into  account  the  effect  of  pressure  on  the  sound  speed 
profile.  This  assumes,  for  lack  of  better  information,  a  well-mixed  water  column  from 


52 


the  surface  down  to  the  ADCP.  which  means  a  constant  profile  of  potential  temperature 
Since  sound  speed  depends  on  in  situ  temperature,  the  sound  speed  will  be  greater  and 
actual  range  will  decrease.  Figure  3.3.1  shows  the  resulting  corrections  for  a  nominal 
range  of  108  m.  Salt  and  pressure  have  a  combined  effect  of  about  5- 10  cm  for  most  of 
the  time  series. 
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Figure  3*3.1  Time  series  of  sound  speed  correction  to  nominal  1 08-m  range  estimate  based  on 
different  methods  described  in  the  text,  where  Ml.  profile  refers  to  the  mixed  laver  profile*  The  thick 
black  line  at  the  bottom  shows  the  approximate  period  of  nearly  complete  ice  cover* 
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The  next  correction  is  site-specific  and  requires  a  simple  model  of  the  upper  layer 
hydrographic  structure.  In  the  presence  of  sea  ice,  the  surface  temperature  value  can  be 
assumed  to  be  at  the  freezing  point  at  the  ambient  surface  salinity.  Using  all  available 
CTD  profiles,  one  constructs  a  functional  form  to  estimate  the  mixed  layer  depth  and  the 
underlying  temperature  and  salinity  profile  that  can  be  empirically  related  to  the 
temperature  and  salinity  measurements  on  the  mooring  (Strass,  1998)*  The  exact 
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structure  of  the  profile  is  not  critical,  but  rather  the  harmonic  mean  of  the  resulting  sound 
speed  profile.  This  correction  must  also  be  applied  to  ULS  sea  ice  draft  estimates,  so  it  is 
a  general  problem  for  ice  draft  estimates  based  on  acoustic  measurements. 

In  the  case  of  Marguerite  Bay,  the  general  features  of  the  wintertime  density 
profile  are  an  upper  mixed  layer  of  water  at  the  surface  freezing  point  overlying  a 
pycnocline  which  warms  towards  the  relatively  warm  Upper  Circumpolar  Deep  Water 
(Smith  et «/.,  1 999).  Therefore,  under  these  circumstances  we  wish  to  find  the  time- 
varying  depth  of  the  mixed  layer  which  is  at  freezing,  and  then  linearly  interpolate  down 
to  the  uppermost  temperature  and  salinity  measurement. 

The  base  of  this  mixed  layer  is  usually  associated  with  high  rates  of  shear 
(Howard  et  at. ,  2004)  so  one  could  use  the  ADCP  velocity  data  to  find  the  mixed  layer 
depth.  We  found  the  ADCP  current  data  too  noisy  to  derive  a  reliable  mixed  layer  depth 
based  on  shear,  and  instead  used  the  temperature  and  salinity  measured  at  99  m.  just 
above  the  ADCP.  Mixed  layer  depths  from  87  CTD  stations  during  a  previous  winter 
hydrographic  survey  were  examined.  A  linear  regression  was  made  between  the 
temperature  at  99  m  and  the  mixed  layer  depth.  The  results  of  this  regression  can  be  seen 
in  Figure  3.3.2,  with  the  fit  to  the  temperature  data  having  an  R:  value  of  0.59.  This 
linear  fit  was  then  applied  to  the  time  series  of  temperature  measured  at  99  m  on  the 
ADCP  mooring,  providing  a  reasonable  time  series  of  estimated  mixed  layer  depth.  The 
salinity  of  the  mixed  layer  was  held  at  33.95  PSU  based  on  examinations  of  temperature- 
salinity  plots  of  the  moored  sensors,  the  hydrographic  surveys,  and  historical  data 
(Hofmann  et  a/.,  1996).  It  should  be  noted  that  this  method  is  applied  to  the  entire  year, 
while  it  is  strictly  only  valid  for  winter  stratification.  During  times  of  open  water  there 
was  often  a  surface  layer  of  warm  fresh  meltwater  that  could  not  be  represented  with  a 
simplistic  empirical  approach.  Therefore,  the  zero  ice  draft  estimates  are  more  uncertain 
because  the  sound  speed  correction  is  most  likely  incorrect  during  times  of  extended  open 
water.  Application  of  the  sound  speed  profile  measured  by  CTD  during  the  mooring 
recovery  reveals  an  overestimation  of  the  ice  draft  by  0.31  m  on  February  21, 2003, 
which  brings  the  estimated  ice  draft  from  0.61  m  to  0.30  ±  0.40  m  and  has  zero-draft 
within  the  uncertainty. 
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Figure  3.3.2  Linear  regression  of  lemperature  measured  at  99  m  with  a  CTI)  and  the  mixed  layer 
depth  i>f  the  cast. 


As  seen  in  Figure  3.3.1,  these  mixed  layer  depth  corrections  have  a  larger  effect 
than  the  salt  and  pressure  corrections  described  above.  The  correction  is  as  much  as  50 
cm  for  May  and  June,  when  sea  ice  is  forming.  During  September  and  parts  of  August 
and  December,  the  mixed  layer  extended  all  the  way  down  to  the  ADCP.  so  the 
correction  is  unnecessary. 

As  a  final  check,  we  have  compared  the  modeled  sound  speed  profile  with  that 
observed  during  the  only  winter-time  CTD  cast  taken  9.3  kilometers  away  on  August  26, 
2002.  As  shown  in  Figure  3.3.3,  the  salinity,  temperature,  density  and  sound  speed 
profiles  match  well.  The  overall  range  correction  is  nearly  identical,  with  the  modeled 
correction  adding  19.9  cm  onto  the  range  and  the  sound  speed  correction  using  the 
observed  profile  adding  20.6  cm. 
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Figure  3.33  Comparison  of  observed  and  modeled  sound  speed  correction  on  August  26,  2002. 
Panels  are  (left  to  right):  temperature,  salinity  and  sound  speed.  The  black  line  is  the  modeled 
profile,  and  the  grey  line  is  the  observed  profile.  Units:  °C,  psu,  and  in/s. 


3.3.2.  Tilt  correction 

The  backscatter  profile  must  be  corrected  for  instrument  pitch  and  roll,  as 
standard  ADCP  processing  uses  the  instrument  tilt  to  binmap  the  velocities  but  not  the 
backscatter  (RD I -Primer,  1996).  Due  to  the  beam  geometry,  this  correction  affects  each 
beam  differently,  and  inherently  collapses  the  separate  estimates  from  each  of  the  four 
beams  towards  one  value.  The  multiplicative  factor  applied  to  the  range  estimate  to 
correct  for  instrument  tilt  can  be  expressed  as 


cos^,±^,(/)  +  ^;),(^(r)  +  f?:))) 
cos(0„) 


(3.6) 


where  0o  is  the  beam  angle  of  the  ADCP,  20°  in  the  case  of  the  RD1  Workhorse,  and  (0r, 
0P)  is  the  roll  or  pitch  angle,  with  the  ±  indicating  that  one  should  add  or  subtract 
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depending  on  the  sign  convention  dictated  by  the  configuration  of  the  sensors  in  the 
instrument.  For  an  upward- looking  RDI  ADCP,  one  should  add  roll  for  beam  1,  subtract 
roll  for  beam  2,  add  pitch  for  beam  3  and  subtract  pitch  for  beam  4. 

The  standard  tilt  sensors  on  an  RDI  Workhorse  ADCP  have  a  resolution  of  0.01°, 
yet  have  an  accuracy  and  precision  of  only  ±  0.5°.  The  latter  represents  about  32  cm  of 
error  for  every  100  m  of  instrument  depth.  To  improve  the  accuracy  of  these 
measurements,  a  constant  offset  (&r\ffp')  is  added  to  each  of  the  time  series  of  roll  and 

pitch.  These  values  are  chosen  to  minimize,  in  a  least  squares  sense,  the  difference 
between  the  range  estimates  from  the  two  beams  affected  by  that  axis  of  the  tilt  sensor  for 
the  entire  time  series.  The  range  of  acceptable  values  for  {0r\0f')\s  constrained  by  the 

accuracy  and  precision  of  the  sensor.  In  this  case,  Bt '  =-0.27°  and  6^ '=0.04°, 

corresponding  to  19  cm  and  3  cm  corrections,  respectively,  on  a  total  nominal  range  of 
108  m. 


3.3.3. Removing  obviously  erroneous  data 

A  series  of  criteria  are  applied  to  find  erroneous  data  and  remove  them  (Strass, 

1998).  Since  many  data  sources  contribute  to  the  ice  draft  estimate,  the  first  step  involves 
acceptable  range  checking  of  any  auxiliary'  data  (temperature,  salinity,  air  or  water 
pressure).  Next,  data  were  eliminated  for  times  when  the  ADCP  ensemble  averaging  in 
space  and  time  and  the  fitting  function  (equation  3.4)  to  the  backscatter  profile  for  a 
single  ADCP  beam  failed  to  accurately  determine  the  range  to  the  surface.  When  the 
standard  deviation  of  the  residuals  of  the  fitted  function  cr.  exceeded  a  chosen  threshold, 

the  range  estimate  was  omitted.  The  threshold  chosen  here  was  6.7  backscatter  units 
(BSU),  which  eliminated  about  10%  of  the  values.  There  was  no  significant  bias  in  the 
residuals. 

Finally  the  difference  between  the  estimates  from  the  two  most  similar  beams, 
dtfif,  as  described  in  section  3.3,  was  examined.  A  maximum  limit,  0.25  m,  for  this 
difference  excluded  8%  of  the  remaining  data.  Overall,  about  18%  of  the  data  was 
removed  due  to  these  combined  tests. 
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3.3.4.  Footprint  error 

A  footprint  error  is  associated  with  the  wide  beam  of  the  ADCP.  This  error  is 
defined  as  the  difference  between  the  measured  draft  and  the  mean  ice  draft  and  does  not 
scale  with  the  ADCP  depth  (Vinje  ei  ai,  1998,  appendix).  The  diameter  of  an  individual 
beam,/j[sin(#  +  ^/2)-sin(0-0?/2)],  is  3.9  m  for  this  ADCP,  which  gives  a  footprint 
error  of  0.4  meters  which  has  been  subtracted  from  the  measured  ice  draft  to  get  the 
expected  mean  draft.  We  have  assumed  the  underwater  roughness  of  ice  in  Marguerite 
Bay  resembles  Arctic  winter  ice  for  purposes  of  the  footprint  error.  Readers  are  referred 
to  Vinje  el  a /,  1998,  (appendix)  to  find  an  estimate  of  footprint  error  given  the  footprint 
diameter  and  winter  or  summer  ice. 

3.3.5.  Uncertainty  Estimate 

The  uncertainty  in  the  range  estimate  for  each  ensemble  can  be  divided  into  that 
due  to  the  empirical  fitting  function  (equation  3.4)  and  that  due  to  the  sound  speed 
correction.  First,  we  describe  each  separately  here,  since  the  former  is  specific  to  the 
ADCP  and  the  latter  is  universal  to  all  acoustic  ice  draft  estimates.  Then  we  discuss  the 
combined  uncertainty. 

3. 3.5.1.  Uncertainty  due  to  the  empirical  fitting  function 

Estimation  of  the  error  in  the  ice  draft  estimate  is  complicated  because  the  range 

is  obtained  from  a  non-linear  fitting  procedure  for  in.  Therefore,  we  performed  a  Monte 
Carlo  simulation  of  the  range  estimate  zoMC  obtained  from  fits  of  equation  (3.4),  in  order 

to  estimate  range  uncertainties  due  to  random  errors  in  the  backscattcr  intensity.  The 
median  values  for  the  one-year  time  series  of  the  coefficients  in  equation  (3.4)  (  A  =  - 
1450  BSU,  C  =  - 1  BSU,  cr=  2,8  m  and  //=  107.5  m)  were  used  to  calculate  a  median 
range  and  six-point  synthetic  profile.  Normally  distributed  random  noise  was  then  added 
to  each  JBS  point  in  the  synthetic  profile  to  test  the  sensitivity  of  the  range  estimate  to 
noise  in  the  backscatter  profile.  The  synthetic  noise  had  zero  mean,  and  a  standard 
deviation  crn  which  we  varied  from  z.ero  to  ten  BSU.  For  each  value  of crn,  the  mean  of 

zoliK  was  indistinguishable  from  that  of  the  median  value  profile.  A  Kolmogorov- 
Smimov  test  confirmed  that  z0  in  the  Monte  Carlo  simulation  had  a  normal  distribution. 
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The  standard  deviation  ofz,,.„.  ,<r;  ,  was  calculated  as  a  function  of  <x,  and 

shown  in  Figure  3.3.4.  u;  b  increases  linearly  with  <r„  ( a.*  ~  0.03  la„  -  0.0019  ,  R2  = 
0.99).  We  then  applied  this  relationship  to  each  of  the  two  most  similar  beams 
g i  ( ,<T;  and  then  estimated  the  total  uncertainty  o-  -  4^o\:  ^  +  <J.  .  in  the  final  range 

estimate  zu  due  to  the  empirical  fitting  in  our  example  data,  a  time  series  and  histogram 

of  which  is  shown  in  the  light  gray  lines  in  Figure  3.3.5.  The  peak  in  the  distribution  at 
0.065  m  gives  an  estimate  of  the  most  likely  uncertainty  due  to  the  empirical  fitting  in  our 
example  dataset.  The  upper  bound  on  uncertainty  lies  at  a  ay  of  0. 1 5  m.  This  appears  to 

be  relatively  constant  during  the  one-year  time  series. 


Figure  3,3*4  Standard  deviation  of  zero-crossings  versus  standard  deviation  of  residuals  ^  of 
synthetic  profiles  created  in  the  Monte  Carlo  simulation  (circles)  and  linear  fit  (line). 
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Figure  3,3,5  Uncertainty  in  ice  draft  estimate.  Upper  panel  shows  time  series,  lower  panel  shows 
histogram  of  time  series.  In  both  panels,  light  gray  line  is  empirical  fitting  uncertainty,  dark  gray 
line  is  sound  speed  uncertainty,  and  black  line  is  total  uncertainty. 


3.3.S.2.  Uncertainty  in  the  Sound  Speed  Correction 

The  range  correction  due  to  sound  speed  carries  an  uncertainty  associated  with  the 

modeling  of  the  sound  speed  profile.  In  our  example,  this  depends  on  the  uncertainty  in 
the  estimate  of  the  mixed  layer  depth.  To  quantify  this,  we  express  equation  (3.5)  as: 


C  C  CT 

-  ADCP  +  ^  ADCF 


"  real 


c  c 

retjl  real 


(3.7) 


where  the  first  term  is  the  ratio  in  equation  (3.5)  and  the  second  term  quantifies  the 
uncertainty.  To  calculate  cr(.  ,  we  added  and  subtracted  the  standard  deviation  <rM/  of 


the  misfit  in  the  linear  regression  shown  in  Figure  3.3.2  (21.8  ni)  to  the  estimated  mixed 
layer  depth  for  each  ensemble.  We  used  these  upper  and  lower  bounds  and  the  method  in 
section  3.1  to  calculate  bounds  on  Creal,  the  difference  being  err  .  A  time  series  and 

histogram  of  the  resulting  uncertainty  in  the  range  1 08  ■  CADCP  ac  JC*ml  is  shown  by  the 
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dark  grey  lines  in  Figure  3.3.5.  The  time  series  shows  the  uncertainty  to  be  about  0.15  m 
from  April  to  July  when  the  mixed  layer  is  well  above  the  ADCP,  and  about  0.025  m 
from  August  to  February  (except  November)  when  the  mixed  layer  extends  down  to  the 
ADCP  and  is  thus  well-resolved.  These  two  periods  correspond  to  the  two  peaks  in  the 
histogram.  The  upper  bound  on  uncertainty  due  to  the  sound  speed  correction  is  0. 1 X  m. 


3.3. 5.3.  Combined  total  uncertainty 

The  combined  total  uncertainty  due  to  the  empirical  fitting  function  and  the  sound 

speed  correction  can  be  expressed  as: 
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where  the  first  term  on  the  right  side  is  the  best  estimate  of  the  range  and  the  second  term 
in  brackets  is  the  total  uncertainty.  The  total  uncertainty  is  comprised  of  two  parts,  the 
first  is  associated  with  the  empirical  fitting  function  and  the  second  with  the  sound  speed 
profile.  The  time  series  and  histogram  of  the  total  uncertainty  and  its  two  parts  are  shown 
by  the  black  line  in  Figure  3.3.5. 

The  time  series  and  histogram  show  that  during  periods  of  low  uncertainty  in  the 
sound  speed  correction  (i.e.,  deep  mixed  layer),  the  uncertainty  associated  with  the 
empirical  fitting  dominates  and  we  calculate  a  total  uncertainty  of  about  0. 1 0  m  ( ±  0.20 
m  95%  confidence  interval).  In  April  to  July,  the  sound  speed  correction  uncertainty  was 
greater  than  that  associated  with  the  empirical  fitting,  making  the  total  uncertainty  around 
0.22  m  (±  0.43  m  95%  confidence  interval).  This  demonstrates  that  the  uncertainty  in  the 
sound  speed  profile  is  of  the  same  importance  as  that  of  the  acoustic  ranging  for  ice  draft 
estimation. 


3.3.6.  Example  results:  Marguerite  Bay 

As  a  first  test  of  the  method,  range  estimates  at  high  frequencies  are  examined. 
The  surface  height  tidal  amplitude  for  Marguerite  Bay  is  1-2  m  which  is  the  same  order 
as  the  ice  draft  and  smaller  than  the  ADCP  bin  size.  One  would  expect  this  signal  to  be 
seen  in  both  the  bottom  pressure  record  and  in  the  range  estimates  from  the  ADCP  to  the 
surface  so  that  the  effect  of  tides  is  removed  when  equation  (3.3)  is  evaluated.  The 
surface  detection  method  was  accurate  enough  to  discern  the  tidal  fluctuations  in  the 
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acoustic  estimate  of  range  to  surface  even  when  the  surface  remains  in  one  ADCP  bin  for 
nearly  the  entire  year.  Figure  3.3.6  shows  the  ADCP  range  estimate  A/  and  the  hydrostatic 
range  h„  for  a  period  of  10  days  of  open  water  in  April,  2002.  The  average  difference  of 
the  two  time  series  is  0. 1 1  m.  The  average  of  the  absolute  value  of  the  difference  is  0. 14 
m,  with  95%  of  those  values  less  than  0.32  m. 


1096 


1094 


109.2 


E 


O 


XI 


109 


E 


&  108.8 
c 


CL 


g  108.6 


< 


108.4 


108,2 


108 


6 


7 


8 


9 


10 


31 


2 


3 


4 


5 


1 


Apr 


Figure  3.3.6  Detail  of  nearly  ice  free  period  in  April,  2002  showing  the  strong  presence  of  tides  in  the 
corrected  ADCP  range,  /r, (black)  when  compared  with  the  tides  derived  from  hydrostatics  Mgray). 

Two  phenomena  are  responsible  for  the  non-zero  ice  draft  estimate  during  times 
of  open  water.  First,  the  ambiguity  of  the  sound  speed  profile  prevents  a  truly  accurate 
range  estimate  during  the  spring  melt,  around  December  and  January  in  Marguerite  Bay, 
when  a  fresh  mixed  layer  with  temperatures  well  above  freezing  forms  which  is 
unresolved  by  our  sound  speed  correction.  Second,  it  is  unclear  how  the  ensemble 
averaging  of  the  ADCP  resolves  the  surface  waves  during  times  of  open  water.  The 
acoustic  return  from  wave  troughs  versus  peaks  may  not  be  identical  and  could  introduce 
some  bias.  Furthermore,  air  bubbles  entrained  during  wave  breaking  could  have  an  effect 
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on  the  acoustic  return.  Because  these  errors  are  based  on  sea  state,  they  arc  difficult  to 
correct.  However,  as  shown  in  Figure  3.3.6  the  method  does  work  welt  in  the  fall  when 
the  spring  melt  layer  is  absent. 

The  final  ice  draft  time  series  is  shown  in  Figure  3.3.7.  The  maximum  drafts  were 
quite  large,  with  extended  periods  of  average  daily  values  of  2  to  4  m  in  August  and 
September.  2002,  This  period  of  thick  ice  ended  abruptly,  but  ice  was  still  present  until 
January,  as  seen  in  Figure  3.2.2.  This  early  arrival  of  multimeter-thick  ice  in  Marguerite 
Bay  in  2002  was  observed  by  the  R  VIB  Nathaniel  B.  Palmer  as  it  was  trapped  in  the  ice 
for  several  days  during  this  period.  The  indication  of  50-100  cm  thick  ice  in  February 
2003  is  most  likely  due  to  a  poor  estimate  of  the  sound  speed  profile  as  a  result  of  an 
unresolved  surface  melt  layer.  It  is  worth  noting  that  a  warm  surface  layer  increases 
sound  speed,  which  causes  the  ADCP  range  to  the  surface  to  decrease  and  apparent  ice 
draft  to  increase,  making  this  explanation  of  the  error  physically  consistent.  Accurate  ice 
draft  estimates  for  springtime  would  be  possible  provided  that  the  mixed  layer 
temperatures  were  measured.  Furthermore,  one  could  potentially  use  short  windows  of 
open  water  to  improve  the  uncertainty  in  the  sound  speed  estimate  during  this  period. 


63 


1 


-9  L-1 


Mar 


Apr  May  Jun 


Jul 


Aug  Sep  Oct  Nov  Dec 


Jan 


Feb 


Figure  3.3.7  lee  draft  observed  in  Marguerite  Bay.  2002.  The  thick  black  line  is  a  24  hour  running 
average  and  the  shaded  area  indicates  the  daily  maximum  and  minimum  values  within  this  average. 


The  time  series  of  ice  draft  allows  one  to  compute  an  ice  draft  distribution  as  seen 
in  Figure  3.3.8.  The  distribution  is  binned  in  1-m  bins  and  shows  the  initial  condition  of 
nearly  ice-free  waters,  which  freeze  up  and  by  September  2002  have  1 8%  2-3  m  ice,  and 
12%  3-4  m  ice. 
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Figure  3.3*H  Monthly  average  draft  distributions  of  sea  ice  observed  in  Marguerite  Bay  in  20(12-2003 
in  l-m  bins. 


3. 3.7. Averaging 

The  temporal/spatial  averaging  of  the  ADCP  also  affects  ice  draft  estimation. 
Figure  3.3.9  shows  a  transect  of  ice  draft  and  thickness  measured  by  drilling  through  the 
ice  every  2  m  along  a  100-m  linear  transect  at  an  ice  camp  9  km  from  the  mooring  site 
during  August  2001  (Perovich  el  a/.,  2004).  This  figure  clearly  illustrates  the  great 
variability  in  ice  draft  on  very  short  horizontal  scales.  The  ensemble  averaging  in  time  of 
the  ADCP  results  in  further  spatial  averaging  as  the  ice  moves  beyond  the  area  averaged 
by  the  combination  of  the  data  from  the  four  beams  into  a  single  estimate  of  draft  for 
each  observation.  Thus,  since  the  ADCP-based  estimates  are  an  areal  average,  they  can 
not  resolve  the  short  spatial  variability  found  in  Marguerite  Bay.  An  ULS  is  needed  to 
sample  this  small-scale  variability.  If  the  ULS  were  moored  at  108  m,  its  beam  width  at 
the  sea  surface  would  be  3.4  m  (1.8°  spread  vertical  beam). 
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Figure  3,3.9  Ice  thickness  transect  measured  at  ice  station  Robert  in  Marguerite  Bay  by  the  Perovich 
group  on  August  5  through  9,  2001,  Ice  thickness  mean  =  71.4  cm;  standard  deviation  -  50.3  cm; 
minimum  =  19  cm;  and  maximum  =  235  cm.  (This  figure  kindly  supplied  by  D.  Perovich.) 


3.4.  Recommendations  for  future  deployments 

Use  of  an  upward-looking  ADCP  to  observe  ice  can  be  improved  by  carefully 
including  additional  measurements.  A  precise  pressure  gauge  is  required  for  the  estimate 
of  ice  draft.  First,  measuring  pressure  P0  close  to  the  ADCP  would  reduce  or  eliminate 

any  possible  error  due  to  mooring  motion.  Second,  locally  measured  atmospheric 
pressure  data,  such  as  that  used  in  this  example,  removes  the  uncertainty  introduced  by 
using  distant  observations  or  atmospheric  reanalysis  products.  Third,  some  plan  for 
estimating  the  sound  speed  profile  would  be  of  great  benefit  and  would  improve  not  only 
ADCP  estimates  but  also  ULS  estimates  of  ice  draft.  Finally,  the  shallower  the  ADCP  is 
deployed,  the  more  accurate  the  ice  draft  estimate  will  be  both  for  decreased  geometric 
errors  and  smaller  errors  in  the  extrapolation  of  sound  speed.  However,  ADCPs  are 
usually  deployed  to  provide  velocity  profiles  over  a  significant  depth  range.  Sound  speed 
depends  more  on  temperature  than  salinity  so  a  simple  thermistor  chain  or  a  winched 
system  such  as  the  HOMER  (HOMing  Environmental  Recorder)  could  provide  an 
improved  sound  speed  correction.  The  latter  would  provide  coverage  throughout  the 
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upper  water  column  and  also  avoid  damage  due  to  deep  ice  keels  or  ice  bergs.  An 
improved  tilt  sensor  in  the  ADCP  would  also  avoid  the  need  for  the  constant  tilt  offsets 
as  described  in  section  3.2.  Finally,  there  is  firmware  available  for  RDI  ADCPs 

(LADCP  bottom-track  mode)  which  examines  the  entire  backscatter  profile  and  locates 
the  range  of  maximum  backscatter  near  the  surface  and  its  peak  intensity. 

3.5.  Conclusions 

Upward- looking  ADCPs  have  been  shown  to  be  viable  tools  for  determination  of 
ice  coverage,  ice  draft  and  drift  velocity  as  shown  in  the  example  case  presented  here 
from  the  western  Antarctic  Peninsula  shelf.  The  techniques  described  provide  additional 
information  about  the  sea  ice  at  a  mooring  location  without,  in  most  cases,  any  additional 
instrumentation. 

Ice  detection  with  the  ADCP  using  windowed  variances  of  vertical  velocity  and 
error  velocity,  and  surface  horizontal  speed  can  be  made  independent  of  additional 
satellite  information.  In  the  example  case,  we  found  signal  correlation  and  backscatter 
intensity  to  be  poor  indicators  of  the  presence  of  ice  at  this  site,  but  this  may  not  hold 
elsewhere. 

Estimation  of  the  ice  draft  using  the  ADCP  requires  corrections  to  the  range  for 
instrument  tilt  and  sound  speed  profile.  The  instrument  tilt  correction  is  greatly  improved 
by  a  constant  offset  which  minimizes  the  difference  between  beams.  This  offset 
introduces  a  3  cm  and  19  cm  correction  to  pitch  and  roll,  respectively.  The  sound  speed 
correction  introduced  an  offset  of  between  about  0  and  55  cm,  with  an  average  magnitude 
of  about  30  cm. 

There  are  two  separate  ways  to  calculate  the  uncertainty  in  the  ice  draft  estimate. 
First,  the  range  estimate  during  times  of  open  water,  as  described  in  section  3.6,  resolves 
the  surface  tides  to  within  ±  0.32  m  95%  of  the  time.  Second,  a  Monte  Carlo  simulation 
detenu ines  the  sensitivity  of  the  zero-crossing  of  the  fitting  function  in  equation  3.4  to 
backscatter  noise.  The  most  likely  uncertainty  associated  with  the  empirical  fitting 
function  is  ±  0. 1 3  m  (95%  confidence  interval).  The  uncertainty  associated  with  the 
sound  speed  correction  had  two  regimes.  One  regime  occurred  during  midwinter  when 
the  mixed  layer  extended  down  to  the  ADCP  and  the  sound  speed  correction  was  nearly 
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zero  with  a  small  uncertainty,  most  likely  0.025  m.  The  other  regime  carries  a  larger 
uncertainty,  most  likely  0.15  m,  due  to  the  ambiguity  in  the  location  of  the  mixed  layer 
depth.  We  calculated  a  most  likely  total  combined  uncertainty  of  ±  0.20  m  (95% 
confidence  interval)  during  midwinter  and  ±  0.43  m  (95%  confidence  interval)  when  we 
had  less  confidence  in  our  assessment  of  the  mixed  layer  depth. 

The  two  independent  estimates  of  the  uncertainty  during  times  of  open  water 
( ±  0.32  and  ±  0.43  m)  are  consistent  and  suggest  an  uncertainty  of  ±  0.40  m  on  the 
estimate  of  the  ice  draft  for  this  example,  or  about  0.4%  of  the  nominal  range  of  the 
ADCP.  This  is  larger  than  that  reported  for  an  optimized  ice  draft  estimate  (4  cm)  using  a 
ULS  deployed  at  a  nominal  depth  of  150  m  in  the  Weddell  Sea  (Strass  1998),  or  0.026% 
of  the  nominal  depth.  In  other  words,  on  a  percentage  basis,  the  ULS  has  one  order  of 
magnitude  less  uncertainty  than  the  ADCP.  The  uncertainty  in  both  cases  depends  highly 
upon  the  modeling  of  the  sound  speed  profile. 

Our  study  indicates  that  the  ADCP  is  not  a  substitute  for  a  dedicated  ULS  for  ice 
detailed  ice  draft  measurements,  but  it  does  give  additional  information  about  both  ice 
and  ocean  velocities.  An  ADCP  is  more  than  adequate  to  observe  the  onset  and  breakup 
of  sea  ice  and  can  provide  estimates  of  ice  draft  and  velocity,  valuable  information  with 
little  extra  effort. 

3.A  APPENDIX:  Data  from  Marguerite  Bay,  western 
Antarctic  Peninsula  shelf 

The  methods  described  in  this  chapter  were  applied  to  a  data  set  from  the  US. 
Global  Ocean  Ecosystems  Dynamics  Southern  Ocean  (GLOBEC-SO)  program,  which 
collected  data  using  a  wide  variety  of  instruments  on  the  western  Antarctic  Peninsula 
(WAP)  shelf,  in  and  around  Marguerite  Bay.  This  work  specifically  utilized  data 
collected  by  an  Automatic  Weather  Station  (AWS),  CTD  casts  from  hydrographic 
surveys,  SSMI  satellite  data  and  one  year  of  data  from  one  subsurface  mooring  as  an 
example. 

The  mooring  (C3  in  GLOBEC-SO  array)  was  deployed  during  the  ice-free  season, 
recording  data  during  the  formation  and  subsequent  breakup  of  the  seasonal  sea  ice.  The 
instrumentation  was  relatively  successful  in  enduring  the  ice  in  Marguerite  Bay. 
Upward-looking  acoustic  ice  profilers  were  deployed  on  some  other  moorings.  However, 


68 


regionally  concurrent  ULS  and  ADCP  ice  draft  measurements  were  not  obtained.  At  the 
mooring  used  in  this  example,  conductivity-temperature  sensors  at  various  depths  gave 
some  data  of  the  evolution  of  the  stratification,  with  the  highest  sensor  at  a  depth  of  99  m. 
The  reader  is  referred  to  the  moored  array  technical  report  (Moffat  el  a}.,  2005)  for  more 
details  about  this  and  other  moorings  and  the  moored  data  collected  in  the  GLOBEC-SO 
program. 

The  ADCP  was  a  standard  300  kHz  broad-band  RD!  Workhorse  deployed  at  a 
nominal  depth  of  108  m,  bin  spacing  of  2  m  and  120  pings  per  ensemble  over  a  30-min 
period  (15  sec  between  pings).  The  VACM  pressure  sensor  (Paine  0-4400  psi)  was  part 
of  a  Vector  Averaging  Current  Meter  (VACM)  deployed  at  a  nominal  depth  of  250  m. 

Two  separate  Automatic  Weather  Stations  (AWSs)  were  deployed  in  Marguerite 
Bay  in  close  proximity  to  provide  redundancy  in  case  of  failure,  as  well  as  the  ability  to 
check  for  failure  by  intercomparison.  The  stations  were  located  on  small  rocky  islets 
named  Kirkwood  Island  and  Dismal  Island  and  had  sensor  heights  of  12  and  25  m  above 
mean  sea  level,  respectively,  and  were  separated  by  29  km.  They  transmitted  their  data 
via  ARGOS  satellite  to  the  Antarctic  Meteorological  Research  Center  at  the  University  of 
Wisconsin,  where  the  data  were  processed  into  a  time  series  of  atmospheric  pressure, 
temperature,  relative  humidity,  and  wind  speed  and  direction  with  an  averaging  interval 
of  10  minutes.  The  data  were  then  edited  and  averaged  to  construct  time  series  with 
hourly  resolution.  The  Kirkwood  and  Dismal  atmospheric  pressure  data  agreed  to  within 
±  1  mb  after  adjustment  for  the  difference  in  sensor  height.  The  Kirkwood  station  was 
slightly  closer  to  the  mooring  (68  km)  so  it  was  used  in  this  study. 

Time  series  of  ice  concentration  derived  from  the  passive  microwave  sensor 
SSM/I  microwave  imaging  satellite  data  were  provided  by  the  Arctic  and  Antarctic 
Research  Center  at  Scripps  Institution  of  Oceanography  in  San  Diego.  Two  separate  time 
series  of  estimates  of  ice  concentration  are  used  because  they  feature  different  spatial 
averaging  and  distance  from  the  center  of  the  estimate  to  the  mooring  location.  In 
addition,  they  use  different  bands  and  algorithms,  and  since  neither  is  clearly  superior, 
both  are  utilized.  The  first  estimate  is  based  on  the  85  GHz  channel,  with  a  pixel 
resolution  of  12.5  km.  The  second  estimate  uses  the  NASA  team  algorithm,  which 
combines  the  37  GHz  channel  with  a  variety  of  other  data  and  has  a  resolution  of  25  km 
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(Cavalieri,  1997).  The  icc  concentration  estimated  at  the  pixel  center  closest  to  the 
mooring  location  was  extracted  about  twice  per  day  and  then  a  daily  average  was  made. 
The  distance  from  the  pixel  center  to  the  mooring  location  is  only  170  m  for  the  85  GHz 
channel,  while  the  distance  for  the  NASA  algorithm  varies,  averaging  about  10  km,  and 
never  exceeding  1 7  km. 
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4.  Chapter  4 


Characterization  of  sea  ice  cover,  motion  and  dynamics 
in  Marguerite  Bay,  Antarctic  Peninsula 

4.1.  Introduction 

The  Southern  Ocean  Global  Ocean  Ecosystems  Dynamics  Southern  Ocean  (SO 
GLOBEC)  project  has  chosen  the  continental  shelf  of  the  western  Antarctic  Peninsula 
(wAP)  for  study  due  to  its  unusually  high  concentration  of  Antarctic  Krill  ( Euphasia 
superha)  and  the  predator  populations  supported  such  as  penguins,  seals,  and  whales.  SO 
GLOBEC  aims  to  detenu ine  what  conditions  favor  the  reproduction,  recruitment,  and 
retention  of  krill  populations,  including  ice  cover  (high  ice  versus  low  ice  years),  on-shelf 
intrusion  of  warmer  nutrient  rich  waters  and  subsequent  vertical  mixing,  and  a  sluggish 
cyclonic  circulation.  A  number  of  features  differentiate  the  wAP  shelf  from  a  typical 
mid-latitude  shelf,  e.g.,  the  North  American  Middle  Atlantic  Bight,  in  addition  to  those 
mentioned  above.  The  wAP  shelf  is  about  twice  as  deep  (500  m,  see  Figure  4.1.1),  and 
has  the  zonal ly  uninterrupted  Antarctic  Circumpolar  Current  (ACC)  flowing  along  the 
shelf  break.  The  surface  waters  above  about  100  m  show  a  strong  seasonal  signal.  They 
become  warmer  and  slightly  fresher  in  summer  as  a  result  of  the  positive  surface  heat 
flux  and  resulting  ice  melt.  In  fall,  they  cool  rapidly  and  become  saltier  as  they  mix  down 
into  the  remnant  Winter  Water.  Finally,  the  shelf  is  nearly  fully  covered  with  ice  in 
winter,  and  almost  totally  ice  free  in  summer. 
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Figure  4.1.1  Marguerite  Hay  region  of  the  western  Antarctic  Peninsula  shelf.  The  moorings  and 
weather  stations  are  shown  by  llieir  name  (e.g..  Cl  and  Dismal,  respectively).  The  thick  black  line 
represents  the  coast  and  approximate  position  of  the  permanent  ice  shelf,  and  the  5tll>-m  and  1000-m 
isobaths  arc  shown,  with  the  latter  being  the  location  of  the  shelf  break. 

The  presence  of  sea  ice  has  a  strong  influence  on  the  vertical  transport  of  heat  and 
momentum  between  the  atmosphere  and  the  ocean.  One  of  the  most  biologically 
productive  polynyas  on  the  continent  exists  in  Marguerite  Bay  (Anri go  and  van  Dijken, 
2003  ).  Furthermore,  the  presence  of  seasonal  sea  ice  is  believed  to  play  a  vital  role  in 
krill  ecosystems  throughout  Antarctica  (Quetin  and  Ross,  1991;  Quetin  et  al.,  1994; 
Siegel  and  Loeb,  1995).  Finally,  ice  production  has  been  found  to  be  the  main  control  on 
upper  ocean  stratification  in  Marguerite  Bay  (Meredith  et  al,  2004).  Therefore,  the 
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characterization  of  the  presence  of  sea  ice,  its  thickness  and  motion  in  Marguerite  Bay  are 
important  to  the  overall  Southern  Ocean  GLOBEC  goals. 

The  sea  ice  in  the  Marguerite  Bay  area  has  been  shown  to  be  strongly  influenced 
by  snow  cover.  Deep  snow  weighs  down  the  ice  and  causes  surface  flooding  as  sea  water 
percolates  up  through  the  ice.  This  leads  to  the  formation  of  snow-ice  which  makes  up  a 
significant  portion  of  the  total  ice  (Perovich  et  al,  2004).  This  is  typical  of  Antarctic  ice, 
which  tends  to  be  deformed  first-year  ice  with  significant  flooding  and  snow-ice 
formation.  The  existence  of  significant  snow-ice  without  accurate  precipitation  data 
prevents  the  estimation  of  brine  rejection  since  snow-ice  introduces  fresh  water  into  the 
sea  ice  budget.  Perovich  et  al  (2004)  also  found  evidence  of  significant  ridging  which 
suggests  the  importance  of  ice  dynamics  in  the  region.  Outside  of  the  Bay,  on  the  outer 
shelf,  Stammerjohn  et  al  (2003)  examined  a  winter  period  (1992)  when  sea-ice  advance 
and  retreat  were  both  anomalously  early.  They  found  that  the  dynamics,  as  opposed  to 
the  thermodynamics,  initiate,  and  thereby  dominate,  the  production  of  ice-edge 
anomalies. 

This  chapter  begins  with  a  brief  description  of  the  unique  data  set  available  from 
the  first  moored  time  series  of  ice  data  on  the  w  AP  shelf,  including  ice  concentration, 
draft  and  motion  for  two  winter  seasons  in  200 1  -2003.  Wc  refer  to  the  austral  winter  of 
2001-2002  as  the  ‘first  season’,  and  the  austral  winter  of  2002-2003  as  the  ‘second 
season’.  Then  some  relevant  aspects  of  the  observed  meteorology  are  presented.  We 
show'  time  series  of  ice  coverage  and  draft  lor  both  seasons.  We  start  with  a  comparison 
of  the  seasonal  cycle  between  the  outer  shelf  and  Marguerite  Bay  during  the  first  year. 

We  then  compare  the  ice  during  the  two  seasons  at  the  same  location  in  the  Bay  and 
examine  both  the  onset  and  break-up  of  the  seasonal  ice  cover.  These  comparisons  along 
with  historical  data  are  used  to  characterize  sea  ice  in  Marguerite  Bay.  Next,  we  examine 
the  motion  of  the  ice  and  how  it  relates  to  the  wind  and  ocean  forcing.  We  accomplish 
this  by  examining  progressive  vector  diagrams  and  calculating  the  full  linear  momentum 
balance.  Finally,  we  consider  the  energetic  near-inertial  band  of  motion  in  both  the  sea- 
ice  and  upper  ocean. 
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4.2.  Dataset 

The  Southern  Ocean  GLOBEC  program  collected  data  using  a  wide  variety  of 
instruments.  This  work  focuses  on  data  collected  by  moored  instruments,  satellites  and 
Automatic  Weather  Stations. 

4.2.1.  Mooring  Data 

A  total  of  nine  subsurface  moorings  were  deployed,  six  the  first  year  2001-2002 
(A I,  A2,  A3,  BI,  B2,  B3),  and  three  the  second  year  2002-2003  (Cl,  C2,  C3),  with  the 
locations  and  names  shown  in  Figure  4.1.1.  The  moorings  were  deployed  during  the  ice- 
free  season  and  recorded  data  during  the  formation  and  subsequent  breakup  of  the 
seasonal  sea  ice.  The  instrumentation  was  relatively  successful  in  enduring  the  ice  on  the 
wAP  shelf.  One  upward-looking  sonar  (ULS)  provides  data  on  the  evolution  of  the  ice 
draft  at  A2.  Up  ward- looking  broadband  acoustic  Doppler  current  profilers  ( ADCPs) 
successfully  measured  currents  and  ice  coverage,  draft  and  velocity  at  the  B2.  Cl  and  C3 
sites.  The  B2  and  C3  moorings  were  collocated  and  give  a  two-year  time  series,  while 
the  Cl  and  C3  mooring  were  separated  by  49  km  and  provide  a  look  at  the  spatial 
variability  within  the  Bay.  The  A2  ULS  allows  a  comparison  between  its  location  at  mid- 
shelf  and  the  B2  site  139  km  away  in  the  Bay  itself.  The  ice  data  were  derived  from  the 
ADCP  data  using  the  method  described  by  Hyatt  et  al  (2005)  (Chapter  3).  The  drafts 
have  an  uncertainty  of  ±  0.20  m  during  midwinter  and  ±  0.40  m  when  the  base  of  the 
mixed  layer  is  above  the  ADCP  (spring  and  fall).  Conductivity  and  temperature  sensors 
at  various  depths  provide  data  on  the  evolution  of  vertical  stratification.  Pressure  sensors 
were  used  to  observe  tides  and  to  remove  atmospheric  effects  when  estimating  ice  draft 
(time  mean  pressure  of  813,  254  and  236  dbar  at  B2,  C3  and  Cl,  respectively).  The  total 
water  depths  at  the  sites  were  808  and  432  at  B2/C3  and  Cl,  respectively.  The  reader  is 
referred  to  the  moored  array  technical  report  (Moffat  et  al,  2005)  for  more  details  about 
the  moored  data  collected  in  the  SO  GLOBEC  program 

4.2.2.  Satellite  Data 

Sea  ice  concentration  maps,  concentrations  at  A2  and  longer  term  data  (1978- 
2002)  were  extracted  from  the  Cavalieri  et  al.  (1997)  Nimbus-7  SMMR  and  DMSP 
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SSM/I  passive  microwave  data  provided  on  CD  by  the  National  Snow  and  Ice  Data 
Center  (NSIDC). 

4.2.3.  Automatic  Weather  Station  Data 

Two  Automatic  Weather  Stations  (AWSs)  were  deployed  in  Marguerite  Bay  in 
close  proximity  to  each  other  to  provide  redundancy  in  case  of  failure,  as  well  as  the 
ability  to  check  for  failure  by  intercomparison.  The  stations  were  located  on  small  rocky 
islets  named  Kirkwood  Island  and  Dismal  Island  and  had  sensor  heights  of  12  and  25  m 
above  mean  sea  level,  respectively,  and  were  separated  by  29  km  (Figure  4.1.1).  They 
transmitted  their  data  via  ARGOS  satellites  to  the  Antarctic  Meteorological  Research 
Center  at  the  University  of  Wisconsin,  where  the  data  were  processed  into  a  time  series  of 
atmospheric  pressure,  temperature,  relative  humidity,  and  wind  speed  and  direction  with 
an  averaging  interval  of  10  minutes.  The  data  were  then  edited  and  averaged  to  construct 
time  series  with  hourly  resolution.  The  Kirkwood  and  Dismal  atmospheric  pressure  data 
agreed  to  within  1  mb  after  adjustment  for  the  difference  in  sensor  height.  The  Dismal 
station  was  slightly  closer  to  the  moorings  (78  km,  12  km  and  38  km  from  A2,  Cl  and 
B2/C3,  respectively)  so  it  was  used  in  this  study.  The  COARE  2.6  algorithm  (Fairall  et 
al,  2003)  was  used  to  calculate  wind  stress.  The  reader  is  referred  to  Moffat  et  al,  (2005) 
for  a  detailed  description  of  the  AWS  instrumentation  and  data  processing,  and  to 
Beardsley  et  al.  (2006)  for  an  in-depth  analysis  of  the  local  meteorology. 

The  meteorological  variables  significant  to  the  sea  ice  are  the  wind  forcing  and  air 
temperature.  The  strong  wintertime  winds  tended  to  blow  towards  the  southwest,  out  of 
the  Bay  during  both  years,  with  light  and  variable  winds  in  between.  Taking  all  wind 
events  with  a  wind  stress  magnitude  greater  than  0.2  N/m~  and  blowing  towards  the 
southwestern  quadrant,  we  find  44  in  the  first  ice  season,  and  23  in  the  second  ice  season. 
These  winds  also  tended  to  be  relatively  warm,  averaging  -2.5°  C  the  first  season  and  - 
3.3°  C  the  second  season,  but  still  cold  enough  for  strong  air-sea  interaction.  The  average 
wind  stress  during  these  events  was  0.5  N/m;  both  years.  Thus,  the  first  observed  ice 
season  had  more  frequent  strong  forcing  events. 
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4.3.  Sea  ice  cover  and  draft. 


4. 3.1.  Austral  winter  2001-2002 

For  the  first  season  a  comparison  between  the  A2  and  B2  sites  was  possible.  The 
A2  mooring  was  placed  on  the  wAP  shelf  61  km  from  the  coast  and  139  km  away  from 
the  B2  mooring,  which  was  deployed  at  the  mouth  of  Marguerite  Bay,  in  the  axis  of  the 
deep  Marguerite  Trough  (Figure  4.1.1).  Figure  4.3. 1  a  shows  ice  drafts,  which  have  been 
boxcar  filtered  with  a  24-hour  window.  The  drafts  are  the  same  (within  experimental 
uncertainties)  between  July  and  early  October.  The  drafts  grow  steadily  to  1  m  by  the 
end  of  July,  nearly  2  m  by  the  end  of  August  and  over  2  m  by  the  end  of  September.  In 
October,  the  ice  continues  to  thicken  at  B2  but  starts  to  retreat  at  A2. 

The  2-minute  sampling  of  the  A2  ULS  shows  the  finer  details  of  the  ice  cover 
(gray  shaded  area  in  Figure  4.3.1a),  including  4  ice  keels  greater  than  10  m  (one  just 
missing  the  ULS  at  45  m),  which  are  presumably  icebergs.  Furthermore,  frequent  deep 
keels  greater  than  5  m  suggest  significant  ice  convergence  and  associated  ridging  at  the 
mid-shelf  site,  as  these  drafts  are  unattainable  with  simple  thermodynamic  ice  production. 
The  B2  site  also  experienced  deep  ice  drafts,  with  daily  averaged  ice  drafts  of  up  to  5  m 
observed  during  late  October/early  November. 
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Figure  43. 1  Daily  averaged  (a)  ice  draft  and  (h)  concentration  during  the  first  season  at  the  A2 
(black)  and  B2  (grey)  mooring  sites;  daily  averaged  (c)  ice  draft  and  (d)  concentration  during  the 
second  season  at  the  Cl  (black)  and  C3  (grey)  mooring  sites.  In  (a),  the  shaded  grey  area  shows  the 
range  of  the  2-minute  measurements  taken  by  the  ULS  at  A2.  The  large  variability  in  range  shown 
before  the  onset  of  ice  in  July  and  after  the  ice  retreat  in  October/Novemher  is  due  to  surface  waves. 
In  (b),  the  ice  concentrations  at  A2  site  were  measured  by  SSMI  satellite.  In  (b)  and  (d),  the  dotted 
lines  represent  ihe  weekly  climatology  of  ice  concentration.  Note  that  (a)-(d)  refer  to  top  to  bottom 
panels,  respectively. 

Figure  4.3.  lb  shows  the  daily  averaged  ice  concentration  at  the  A2  and  B2  sites, 
with  a  weekly-binned  climatology  derived  from  the  NSIDC  data  plotted  as  dotted  lines. 
The  ice  concentrations  at  the  A2  site  were  measured  by  SSMI  satellite,  as  the  ULS  does 
not  measure  ice  concentration  in  the  same  manner  as  the  ADCP  (Hyatt  et  al,  2005).  The 
ice  covered  the  mooring  sites  at  about  the  same  time,  but  the  ice  remained  for  about  2.5 
months  longer  in  the  Bay  near  B2.  The  meteorological  forcing  measured  by  the  AWSs 
within  the  Bay  is  not  necessarily  representative  of  that  on  the  shelf.  Therefore  we  cannot 
confidently  explain  the  shorter  ice  coverage  at  A2.  However,  based  on  the  long-term 
NSIDC  ice  concentration  data,  ice  typically  lasts  approximately  1  to  2  months  longer  at 
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B2  than  A2  (dotted  lines  in  Figure  4.3.1b).  At  B2,  the  ice  covered  the  mooring  site  in 
early  July  and  remained  until  mid-December.  The  ice  cover  in  late  December  and 
January  of  the  first  season  was  isolated  to  southern  Marguerite  Bay,  as  revealed  by  the 
NSIDC  satellite  data  (not  shown),  just  barely  reaching  the  mooring  site.  Although  ice 
was  not  present  at  the  mooring  site  for  December/January  of  the  previous  year  in  the 
NSIDC  data,  the  satellite  record  shows  it  frequently  extended  to  the  mooring  site  at  that 
time  of  year. 

Overall,  the  sites  are  similar  until  mid-October,  when  the  ice  retreats  from  the 
mid-shelf  site  (with  a  brief  return  a  few  weeks  later),  but  thickens  dramatically  at  B2. 

The  ice  retreats  from  B2  in  mid-December,  only  to  return  again  for  a  short  period.  This 
produces  a  climatologically  short  ice- free  season  in  the  Bay,  as  described  later  in  this 
chapter. 

4.3.2.  Austral  winter  2002-2003 

In  the  second  season,  two  moorings  (Cl  and  C3)  in  the  mouth  of  Marguerite  Bay 
provide  time  series  of  ice  cover  and  draft  with  a  49-km  separation.  Cl  was  further  into 
the  Bay,  almost  due  east  of  C3  (Figure  4.1.1).  Figure  4.3.1c  shows  ice  draft  at  the  two 
sites,  which  have  been  boxcar  filtered  with  a  24-hour  window,  and  Figure  4.3.  Id  shows 
the  daily  averaged  ice  concentrations  with  the  NS1  DC-based  weekly  climatology  shown 
by  the  dotted  lines. 

At  Cl.  the  ice  covered  the  mooring  site  in  mid-May  and  remained  until  late 
December,  with  daily  averaged  ice  drafts  of  up  to  4  m  observed  during  late  October/earfy 
November.  At  C3,  ice  also  covered  the  mooring  site  in  mid-May  and  remained  until  late 
December.  There,  daily  averaged  ice  drafts  of  2-3  m  were  observed  from  mid-August 
through  mid-October.  Around  October  12  the  ice  broke  up,  and  the  ice  draft  dropped 
dramatically  to  ~l  m  at  both  sites.  In  addition,  the  ice  concentration  dropped  somewhat 
at  both  sites,  but  still  remained  above  60%  for  most  of  the  remaining  ice  season. 

The  drafts  at  Cl  and  C3  are  the  same  (within  observational  uncertainties)  except 
for  a  period  in  August  with  5-6  m  drafts  at  the  C3  mooring  and  from  mid-September 
through  mid-October  when  the  ice  draft  at  C3  was  1-2  m  larger  than  at  Cl.  Ice  first 
appears  at  the  C3  mooring  3  days  after  the  Cl  mooring,  and  retreats  from  the  C3  site  7 
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days  before  Cl.  The  ice  cover  at  the  two  sites  has  a  correlation  coefficient  of  0.85,  and 
the  ice  draft  has  a  correlation  coefficient  of  0.67. 

A  polynya  appeared  for  a  short  period  at  the  Cl  site  around  October  12,  2002,  but 
did  not  extend  to  C3.  The  polynya  was  also  apparent  in  satellite  imagery  (Figure  4.3.2a) 
which  showed  a  short-lived  area  of  low  ice  cover.  The  moored  observations  clearly  show 
that  the  polynya  did  not  extend  to  the  C3  mooring  site,  despite  the  ambiguity  in  the 
NSIDC  satellite  data.  This  event  is  described  in  greater  detail  in  section  3.4. 
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of  low  ice  concentration  around  the  Cl  site,  which  did  not  extend  to  the  C3  site,  (b)  Weekly 
climatological  transect  of  ice  coverage  in  at  siv  points,  with  the  colors  corresponding  to  locations  of 
the  colored  dots  in  (a).  The  climatology  has  been  repeated  for  clarity  and  is  based  on  the  period  of 
1 978-2002. 
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4.3.3. Interannual  comparison  of  sea  ice  cover  and  draft 

The  collocation  of  the  B2  and  C3  moorings  provides  an  intcrannual  comparison  of 

the  sea  ice  cover  and  draft  independent  of  any  spatial  variability.  Both  years  have 
roughly  seven  months  of  nearly  complete  ice  cover,  but  the  onset  of  ice  formation  is 
about  two  months  earlier  in  2002-2003  than  in  the  previous  year  (Figure  4.3. 1  a  and 
4.3.1c).  The  early  ice  in  2002-2003  can  be  explained  by  the  large  negative  surface  heal 
flux  associated  with  cold  air  temperatures  which  came  earlier  than  normal  in  early  May. 
As  shown  in  Figure  4.3.3, 2003  had  the  coldest  May/June  on  record  (since  1976),  The 
climatology  is  based  on  temperature  recorded  at  the  British  Antarctic  Survey’s  Rothcra 
station  on  the  northern  shore  of  Marguerite  Bay,  and  is  a  good  proxy  for  temperature  in 
the  center  of  the  Bay  (Beardsley  ct  al.,  2006).  In  both  years,  the  ice  cover  gets  above 
60%  only  after  the  air  temperature  goes  below  -9°  C.  To  further  test  the  relationship 
between  the  first  cold  event  in  austral  autumn  and  the  onset  of  ice,  we  uti  lized  the  1 977- 
2003  air  temperature  data  from  Rothcra  station  and  the  NSIDC  ice  concentration  data. 

We  compared  the  first  date  of  -9.2°  C  in  a  3-day  filtered  temperature  record  and  found  it 
correlated  to  the  first  time  that  greater  than  70%  3-day  filtered  ice  concentration  occurred 
at  the  B2/C3  site  with  a  correlation  coefficient  of  0.8 1 .  The  standard  deviation  of  the 
mismatch  in  the  least-squares  linear  regression  is  1 8  days,  This  suggests  that  local 
surface  cooling  is  the  dominant  process  in  ice  formation  in  the  region. 


82 


.35  L-1 _ i _ _ 1- _ J _ i _ i _ i _ I _ i _ l _ -i- - 

Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov  Dec  Jan  Feb 

Figure  4.3.3  Air  temperatures  (24-hr  box-filtered)  measured  at  the  AWS  during  the  first  season 
(gray)  and  second  season  (black).  The  dotted  black  line  is  the  air  temperature  in  2001  recorded  at 
the  British  Antarctic  Survey  Rothera  station  before  the  AWS  was  deployed.  The  monthly 
climatology  based  on  measurements  at  Rothera  for  the  period  1976-2003  is  plotted  as  a  dashed  black 
line,  with  the  light  gray  shaded  area  showing  plus  and  minus  one  standard  deviation. 


A  climatology  of  ice  cover  along  a  transect  from  within  Marguerite  Bay  to 
beyond  the  shelf  break  provides  further  evidence  of  local  ice  formation  (Figure  4.3.2b). 
The  ice  on  the  shelf  forms  at  nearly  the  same  time  in  June,  with  the  off-shelf  points 
following  about  two  weeks  later.  The  more  gradual  disappearance  of  ice,  with  the 
inshore  points  lagging  the  offshore  ones,  suggests  that  the  ice  docs  not  melt  locally,  but 
rather  breaks  up  and  is  adveclcd  away.  A  final  piece  of  evidence  which  suggests  local 
formation  of  ice  in  the  Marguerite  Bay  region  in  the  second  season  is  found  in  the  net 
surface  heat  flux  (the  sum  of  shortwave,  longwave,  sensible  and  latent  fluxes)  calculated 
from  the  shipboard  meteorological  data  collected  on  the  RV1B  N.B.  Palmer  during  a 
survey  in  the  area  (Figure  4.3.4).  This  shows  that  in  the  time  preceding  the  onset  of  ice, 
the  net  surface  heat  flux  averaged  about  40  W/m~  until  May  3,  2002.  Then  there  were  six 
days  of  average  heat  loss  in  excess  of  1 80  W/m 2  at  which  point  the  ice  began  to  form, 
followed  by  six  more  days  with  an  average  heat  loss  of  260  W/m2.  If  we  assume  that  this 
heat  loss  went  solely  into  ice  production,  the  amount  of  ice  formed  during  the  May  3-9 
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period  would  be  ~  0.5  m  thick  and  during  the  May  3-16  period  would  be  ~  1.5  m  thick. 
These  large  heat  losses  coinciding  with  the  onset  of  sea  ice  imply  local  formation  rather 
than  advcction  into  the  region. 
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Figure  4.3.4  (a)  Air-sea  heat  flux  calculated  from  meteorological  sensors  onboard  the  Rl  'IB 
Palmer ,  which  was  surveying  the  area  and  (b)  sea  ice  concentration  at  C3  during  the  same  period  in 
2002. 


Both  years  have  very  thick  ice,  with  extended  periods  (including  one  continuous 
month)  with  drafts  over  3  m.  The  thick  ice  occurs  about  one  month  earlier  in  2002  than 
in  2001 ,  which  can  be  attributed  to  the  early  ice  season.  However,  the  month  of  highest 
ice  draft  had  1  -2  m  thicker  draft  in  2001 . 

The  combination  of  late  ice  in  Dcccmbcr/January  at  the  B2  site  and  early  ice 
observed  at  the  C3  site  made  for  a  relatively  short  ice-free  period  of  only  3  months  at  the 
B2/C3  site.  Historical  NSIDC  satellite  data  since  1978  show  only  four  shorter  ice-free 
seasons  in  24  years.  In  1980,  ice  remained  at  the  B2/C3  site  for  the  entire  summer. 
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4.3.4.0ctober  12,  2002  Polynya  Event 

On  October  12, 2002,  a  polynya  opened  up  at  C 1  which  did  not  extend  to  C3 
(Figure  4.3.2a).  A  similar  polynya  was  not  observed  in  the  NSIDC  satellite  data  in  200 1 . 
The  open  area  was  about  50  km  across  and  occurred  during  a  strong  westward  wind 
event.  This  wind  event,  however,  was  not  exceedingly  strong  for  the  region  ( Beardsley  ct 
ah,  2006).  The  meteorological  data  and  ice  conditions  surrounding  the  event  arc  shown 
in  Figure  4.3.5.  The  air-sea  heat  flux  is  the  sum  of  the  sensible  and  latent  flux  (Fairall  ct 
al,  2003)  weighted  by  the  open  water  fraction  (assuming  a  sea  surface  temperature  of  - 
1 .8°C)  and  is  meant  to  be  a  gross  estimate  to  show  that  there  is  a  large  heat  flux  from  the 
ocean  to  the  atmosphere.  Although  not  measured,  the  net  radiative  heat  flux  (long  and 
shortwave)  was  presumably  negative  and  possible  quite  large,  based  on  the  shipboard 
measurements  shown  in  Figure  4.3.4(a),  and  would  give  an  even  larger  negative  surface 
heat  flux. 
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Figure  4.3.5  (a)-(d)  AWS  measurements  surrounding  the  period  of  open  water  on  October  12,  2002. 
Sensible  and  latent  heat  flux  (c)  into  the  ocean  weighted  by  the  open  water  fraction  (see  text  for 
details),  and  ice  concentration  (0  and  draft  (g)  derived  from  the  ADCP  at  Cl.  Note  that  (a)-(K)  refer 
to  top  to  bottom  panels,  respectively. 


The  upper  ocean  at  the  Cl  site  responded  strongly  to  the  event.  Direct  exposure 
to  the  strong  westward  winds  initiated  a  strong  westward  current  of  15  cm/s  at  6  m  below 
the  bottom  of  the  ice,  which  is  quite  energetic  for  the  wintertime  in  the  region.  The  upper 
150  m  cooled  approximately  0.5°  C  in  these  two  days,  with  effects  observed  down  to  200 
m  (Figure  4.3.6).  The  strong  cooling  at  the  surface  was  presumably  caused  by  the  large 


86 


net  air-sea  heat  flux  with  the  cold  (-10°  to  -5°  C)  atmosphere  and  strong  winds.  This  was 
followed  on  October  14  by  a  sharp  warming  of  about  1°  C.  This  warming  could  be 
associated  with  up  welling  along  the  new  ice  edge  due  to  a  wind  burst  towards  the  north 
(e.g.  Rood  and  Obrien,  1983). 


Cl  moonng,  Oct  10-12  ice  opening 


2002.  The  depth  of  each  temperature  record  is  labeled  on  the  right. 

One  possible  explanation  is  that  the  ice  broke  along  the  edge  of  thicker  ridged  fast 
ice  within  the  Bay.  The  ice  draft  at  C3  declined  significantly  during  the  event,  further 
evidence  that  the  ice  broke  up.  Evidence  of  this  can  be  seen  in  two  MODIS  satellite 
images  of  the  wAP  shelf  on  (a)  October  24,  2001  and  (b)  September  6,  2002  (Figure 
4.3.7).  The  interior  of  the  Bay  contains  ice  which  appears  nearly  solid.  The  mouth  of  the 
Bay  and  beyond,  while  having  very  high  ice  concentrations,  contains  large  pieces  of  ice, 
with  areas  of  broken  ice  in  between,  which  could  lower  internal  ice  stresses  and  allow 
greater  ice  motion. 
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Figure  4J.7  MODIS  satellite  images  of  the  region  on  (a)  October  24,  2001  and  (b)  September  6,  2002 
and  the  approximate  locations  of  the  moorings. 


4.4.  Ice  Motion 

lec  motion  was  measured  using  an  ADCP  at  the  B2,  C I  and  C3  sites.  Some 
statistics  of  this  motion  arc  given  in  Table  I .  The  statistics  do  not  vary  greatly  from  year 
to  year  and  site  to  site.  The  average  distance  traveled  in  one  day  by  the  ice  ranged  from 
9.6  km  at  the  C3  site  to  10.5  km  at  the  B2  site,  with  an  average  icc  speed  ranging  from 
10.5  cm/s  at  Cl  to  12.9  cm/s  at  B2.  This  motion  was  rarely  unidirectional,  so  that  the 
time-averaged  net  distance  traveled  per  day  was  6.6  and  7.6  km  at  the  Cl  and  C3  sites, 
and  7.8  km  at  the  B 2  site. 
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Mooring 

Distance  traveled  in  a 
day  (km) 

Average  Daily 
Displacement  (km) 

Mean  Daily-averaged 
speed  (cm/s) 

Days 
of  icc 

cover 

>  70% 

Mean 

Sid 

Min 

Max 

Mean 

Std 

Min 

Max 

Mean 

Std 

Min 

Max 

B2 

10.5 

6.0 

1,6 

30.6 

7.8 

5.7 

0.2 

28.5 

12.9 

7.4 

1,9 

38,1 

122 

Cl 

8.8 

8.2 

1.5 

44.3 

6.6 

7.5 

0.2 

39.2 

10.5 

9.8 

1.7 

53.0 

189 

C3 

9.6 

7.0 

1.2 

36,0 

7.6 

6.4 

0.1 

35.5 

1 1.6 

8.3 

1-5 

42J 

185 

Tabic  1:  Icc  velocity  statistics. 


Another  way  to  present  the  motion  of  the  icc  and  upper  ocean  is  to  use  a 
progressive  vector  diagram  (PVD).  This  integrates  the  Eulcrian  velocity  in  time 

assuming  that  it  is  spatially  homogeneous,  i.c.  x  =  j  udt,  to  give  a  synthetic  particle 

Ji 

trajectory.  During  the  first  season,  the  more  frequent  occurrence  of  strong  southwestward 
winds  caused  a  net  southward  icc  motion  at  B2  (Figure  4.4. la).  At  C3  during  the  second 
season,  during  the  period  of  nearly  solid  ice  pack  (June  2 1 , 2002  through  December  20, 
2002)  the  icc  had  almost  no  net  motion,  with  the  displacements  remaining  in  a  roughly 
1 00  km  by  1 00  km  box  (Figure  4.4. 1  b).  The  southward  wind  events  push  the  icc 
southward,  then  the  slow  but  steady  northward  ocean  current  returns  the  ice  to  near  its 
position  prior  to  the  wind  event.  This  suggests  that  the  same  icc  remains  in  the  C3  area 
for  6  months  during  the  second  season. 
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Figure  4.4.1  Progressive  vector  diagrams  of  (a)  B2  (b)  C'3  and  (c)  Cl  with  ice  shown  in  blue  and  the 
upper  ocean  shown  in  red.  The  time  of  isolated  very  low  (<10%)  ice  concentration  at  Cl  is  identified 
by  the  dashed  line.  The  integrations  begin  at  (0,0). 


The  nearby  Cl  site  exhibits  entirely  different  behavior  (Figure  4.4.1  c).  There,  the 
ocean  has  a  net  southward  drift,  in  the  same  direction  as  the  strong  winds.  This  suggests 
a  strong  southwestward  ice  drift  over  the  winter,  implying  that  the  same  ice  docs  not 
remain  at  this  site  for  the  winter. 

Inner  rotary  coherence  spectra  (Emery  and  Thomson,  200 1 )  of  the  ice  at  C 1  and 
C3  for  the  time  series  shown  in  Figure  4.4.1b  and  c  show  that  the  near-inertial  motions 
arc  coherent  in  the  ice.  The  left-turning  rotary  coherence  squared  in  the  frequency  band 
from  12.2!  to  13.75  hours,  centered  at  12.89  hours,  (the  inertial  period  at  69  S  is  12.81 
hours)  is  0.59  with  zero  phase  lag.  The  right-turning  rotary  coherence  for  the  same  band 
has  low  coherence,  as  do  higher  frequencies  in  both  directions. 
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4.5.  Momentum  Balance  in  the  Sea  Ice 


A  momentum  balance  vertically  integrated  from  the  bottom  to  the  top  of  the  ice 
can  be  estimated  from  the  data  available.  Neglecting  nonlinear  advection,  we  represent 
this  balance  as 


where  p.  is  the  ice  density,  taken  to  be  constant  at  920  kg/m3,  assuming  the  ice  is  in 
hydrostatic  balance,  the  ice  thickness  ht  is  the  measured  draft  divided  by  the  density  ratio 
of  ice  to  seawater,  u,  is  the  ice  velocity,  f  is  the  Coriolis  parameter,  is  the  sea 
surface  slope,  A  is  the  fractional  area  ice  coverage,  r ,  and  r„  arc  the  air  and  water 

stresses  on  the  ice,  respectively,  and  Fh  is  the  ice  dynamics  term  (Hiblcr,  1979).  Using 
geostrophy,  the  surface  slope  can  be  expressed  in  terms  of  the  ocean  gcostrophic  flow 


Assuming  no  gcostrophic  shear ,  we  used  the  velocities  at  a 


g 


deep  current  meter  (400  m,  232  m  and  247  m  at  B2,  Cl  and  C3,  respectively)  as  a  proxy 
for  ua,  which  is  accurate  to  within  1 0%  at  C3,  the  only  mooring  where  shear  was 

measured.  The  water  stress  ru,  is  calculated  using  a  quadratic  drag  law  with  a  drag 

coefficient  of  0.015  (based  on  Table  1  in  Shirasawa  and  Ingram,  1991).  No  turning 
angles  were  used  in  the  air  or  water  stresses  due  to  the  close  proximity  of  the 
measurements  to  the  ice.  Next,  we  assume  complete  ice  cover  and  estimate  the 
divergence  of  the  internal  ice  stresses  as  the  residual  sum  of  the  directly  estimated  terms, 

or  A  *  I,  Fh  =  V  •  <y .  Therefore,  we  only  estimate  the  balance  when  the  measured  ice 

concentration  exceeds  70%.  This  leaves  us  with  an  expression  for  the  momentum 
balance  in  the  ice,  where  the  residual  of  the  measured  terms  is  the  divergence  of  the 
internal  ice  stresses. 


Since  the  residual  of  the  measured  terms  includes  both  the  internal  ice  stresses  as  well  as 
the  sum  of  the  errors  in  the  estimated  terms,  this  residual  estimate  should  be  used 
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cautiously.  We  refer  to  the  second  term  on  the  left  hand  side  as  the  Coriolis/tilt  term.  The 
time  series  used  to  constnict  the  momentum  balance  terms  were  low-pass  filtered  with  a 
33-hour  window. 

The  momentum  terms  derived  at  B2  show  a  primary  balance  between  wind  stress 
and  internal  ice  stresses,  with  ocean  stresses  taking  a  primary  role  at  times,  particularly 
early  in  the  season  (Figures  4.5. 1  and  4.5.2).  In  July  and  early-August,  strong  ocean 
stresses  on  the  ice  are  balanced  primarily  by  internal  ice  stresses.  Then,  beginning  in 
mid-August,  there  were  many  strong  southward  wind  bursts  (Beardsley  et  al.,  2006), 
which  are  generally  balanced  by  internal  ice  stresses.  Overall,  wind,  ocean  and  ice 
stresses  form  the  dominant  balance,  the  Coriolis/tilt  term  plays  a  minor  role  at  times  and 
the  acceleration  term  is  not  significant. 
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Figure  4.5.1  Momentum  terms  in  ice  at  the  B2  site.  The  variance  of  each  term  is  shown  in  the  upper 
right  of  each  axis. 
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Figure  4.5.2  Monthly  averaged  and  overall  mean  momentum  balance  at  112. 


The  following  season  at  C3,  the  primary  balance  is  between  wind  stress  and 
internal  ice  stresses,  with  ocean  stresses  taking  a  primary  role  at  times  (Figures  4.5.3  and 
4.5.4).  In  2002  there  were  less  frequent  southward  wind  bursts,  and  all  of  the  terms  in  the 
balance  are  generally  smaller.  Again,  earlier  in  the  season  the  ocean  stress  is  more 
important.  Overall,  wind  and  ice  stresses  form  the  dominant  balance. 
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Residual  Coriolis  Temporal  Ocean  Stress  Wind  Stress 


Momentum  terms  at  C3 
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2002  2003 

Figure  4.5.3  Momentum  terms  in  ice  at  the  C3  site.  The  variance  of  each  term  is  shown  in  the  upper 
right  of  each  axis. 
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Figure  4.5.4  Monthly  averaged  and  overall  mean  momentum  balance  at  C3. 


At  the  nearby  C 1  site,  we  once  again  see  a  primary  balance  between  wind  stress 
and  internal  ice  stresses,  with  ocean  stresses  taking  a  primary  role  at  times  (Figures  4.5.5 
and  4.5.6).  In  particular,  we  see  this  early  in  the  season,  and  then  again  in  October  and 
November.  As  at  the  C3  site,  wind,  ocean  and  ice  stresses  form  the  dominant  balance. 

The  ocean  stress  was  more  important  at  Cl  than  C3.  At  both  Cl  and  C3  the  Coriolis/tilt 
term  made  a  small  contribution  and  the  acceleration  term  was  essentially  negligible.  The 
variance  of  each  term,  an  indicator  of  its  importance,  is  shown  in  the  upper  right  comer  of 
each  axis  for  all  terms  at  all  sites. 
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Residual  Coriolis  Temporal  Ocean  Stress  Wind  Stress 


Momentum  terms  at  Cl 


Figure  4.5.5  Momentum  terms  in  ice  at  the  Cl  site.  The  variance  of  each  term  is  shown  in  the  upper 
right  of  each  axis. 
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Figure  4.5.6  Monthly  averaged  and  overall  mean  momentum  balance  at  C3. 


4.6.  Near-  Inertial  Oscillations  in  Sea  Ice  and  the 
Upper  Ocean 

Thick  (~  1  m)  sea  ice  acts  an  effective  insulator,  reducing  air-sea  heat  flux  by  one 
to  two  orders  of  magnitude  from  open  water  values.  Therefore,  the  overall  air-sea  heat 
flux  budget  can  be  controlled  by  small  areas  (as  little  as  1%)  of  open  water  or  thin  ice  in 
the  icepack.  High  frequency  motions,  such  as  near-inertial  or  tidal  oscillations,  are  often 
spatially  inhomogeneous,  with  alternating  convergence  and  divergence.  This  differential 
ice  motion  can  open  and  close  leads  on  time  scales  shorter  than  the  thermodynamic  time 
scale  for  freezing,  thus  providing  insufficient  ice  strength  to  prevent  relative  ice  motion: 
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the  new  ice  is  crushed  and  deposited  around  floe  edges.  Heil  et  a).,  (2001 )  realized  an 
increase  in  ice  thickness  of  about  13%  over  a  5-month  model  simulation  by  inclusion  of 
moderate  high  frequency  (subdaity)  deformation  dynamics  in  a  dynamic/thermodynamic 
ice  model.  Heil  and  Hibler  (2002)  included  “inertial  embedding”  (McPhee,  1978)  into 
the  dynamics  of  the  oceanic  boundary  layer  of  a  Arctic  sea  ice  and  found  not  only 
improved  reproduction  of  observed  ice  motion,  but  an  increase  in  winter  ice  mass  of 
about  20%.  Furthermore,  shear  at  the  base  of  the  mixed  layer  and  convection  during  ice 
production  associated  with  high  frequency  ice  motion  could  enhance  oceanic  vertical 
mixing.  This  motivates  the  characterization  of  observed  near-inertial  motions  in  the  sea 
ice  and  upper  ocean. 

Power  density  spectra  of  ice  and  ocean  velocities  at  B2,  Cl  and  C3  (Figure  4.6. 1 , 
dashed  lines  show  total  spectra)  show  peaks  at  a  period  of  about  12-13  hours.  There  is 
more  energy  in  this  band  in  the  ice  than  in  the  water,  C3  had  the  largest  peak,  with  Cl 
nearly  as  large.  The  peak  at  lower  frequencies  is  the  synoptic  3-4  day  weather  band.  The 
first  year  (B2)  had  a  less  pronounced  spectral  peak  near  the  inertial  frequency,  despite  the 
more  frequent  occurrence  of  strong  wind  forcing  events.  This  apparent  contradiction 
could  possibly  be  explained  by  the  fact  that  near-inertial  coupling  requires  anti -clockwise 
rotation  in  the  wind  stress  vector,  as  well  as  strong  forcing. 
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Rotary  and  total  spectra  at  B2  of  ice  (black)  and  water  (gray) 


Figure  4.6.1  Rotary  and  total  spectra  at  (upper  panel)  B2  (middle  panel)  Cl  and  (lower  panel)  C3  for 
the  same  time  periods  as  shown  in  Figure  4.4.1.  The  ice  is  shown  in  black  and  upper  ocean  in  gray, 
with  solid  lines  showing  the  anticlockwise  component  and  dashed  lines  the  clockwise  component.  The 
short  vertical  lines  with  circles  show  the  95%  confidence  interval  at  three  points  for  the  ice  (black) 
and  ocean  (gray).  The  solid  black  vertical  line  shows  the  inertial  frequency  and  the  dashed  black 
vertical  line  shows  the  M2  tidal  frequency. 

At  these  latitudes,  the  near-inertial  band  and  the  M2  semidiurnal  tide  have  very 
similar  periods  of  12.92  and  12.42  hours,  respectively.  Here  vve  show  that  the  energy  is 
indeed  near- inertial.  The  rotary  spectra  at  C3  show  that  the  spectral  energy  peak  is  almost 
entirely  in  the  anti-clockwise  direction  in  both  the  ice  and  water,  consistent  with  near- 
inertial  motion  (Figure  4.6.1).  Furthermore,  near-inertial  motion  is  wind-forced  and 
therefore  intermittent,  while  the  tide  should  be  present  throughout  the  time  scries  and 
have  a  constant  phase.  Figure  4.6.2b  shows  that  during  times  of  high  measured  kinetic 
energy,  the  kinetic  energy  accounted  for  by  counter-clockwise  complex  demodulation  at 
the  inertial  period  is  also  high  and  not  phase  locked  (Figure  4.6.2a).  These  energetic 
periods  have  rms  amplitudes  up  to  -30  cm/s.  A  tidal  analysis  (Pawlowicz  et  al,  2002)  of 
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the  same  period  shows  the  dominant  semidiurnal  component  M2  having  current 
amplitudes  of  only  ~2  cm/s.  Similar  results  also  hold  at  B2  and  Cl  (not  shown  for 
brevity).  Therefore  we  conclude  that  the  motion  is  due  to  near- inertial  oscillations. 


(a) 

CCW  demodulation  for  ice  at  c3 


(b) 


600- 

500- 
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Figure  4.6.2  Counterclockwise  complex  demodulation  of  ice  velocity  at  C3  at  the  near-inertial 
frequency  showing  (a)  phase  and  (b)  measured  kinetic  energy  in  gray  shading  and  kinetic  energy' 
accounted  for  by  the  complex  demodulation  in  black. 

An  example  of  strong  near-inertial  oscillations  in  the  second  year  on  October  28- 
3 1 , 2002  is  shown  in  Figure  4.6.3.  There  is  a  strong  wind  event  on  October  26-28,  2002 
with  an  anti-clockwise  wind  shift  followed  abruptly  by  calm  conditions,  causing 
oscillations  in  both  the  sea  ice  and  upper  ocean.  The  kinetic  energy  per  unit  volume  in 
the  near-inertial  band  is  greater  in  the  sea  ice  than  in  the  upper  ocean.  The  internal  ice 
stresses  are  low  during  this  period,  with  about  75%  of  the  kinetic  energy  in  the  ice 
reproduced  with  near-inertial  complex  demodulation.  Then  on  October  3 1 ,  a  strong 
internal  ice  stress  (some  non-local  process)  causes  northward  ice  velocity  and  is  partially 
balanced  by  ocean  stress  and  the  oscillations  cease.  The  high-pass  filtered  momentum 
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balance  during  this  period  using  the  method  described  in  section  4.5  (not  shown)  does  not 
show  a  simple  inertial  wave  balance  between  the  temporal  change  and  coriolis  terms. 
There  are  also  significant  contributions  at  high  frequencies  from  the  wind  stress  and 
internal  ice  stresses. 


Figure  4.6.3  An  example  of  a  near-f  event  at  the  C3  mooring. 


The  effects  of  the  initial  wind  burst  which  caused  the  ncar-inertial  oscillations  can 

o 

be  seen  down  to  200  m  (figure  4.6.4).  The  heat  content  pc p  jTdz  for  each  mooring  is 

calculated  using  linear  interpolation  of  the  temperature  data  and  assuming  that  the  surface 
temperature  is  at  the  freezing  point  (Figure  4.6.4b).  A  possible  scenario  is  as  follows: 
the  Cl  site  had  60-90%  ice  cover  during  this  period,  so  leads  could  presumably  open  and 
close  with  the  inertial  period.  We  see  strong  cooling  (loss  of  heat  content)  and  mixing 
during  the  initial  wind  burst  on  October  27  which  extends  down  to  200-m.  Then  the 
near-inertial  motions  from  October  28  to  November  1  possibly  mixed  warm  water  back 
up  to  107-m  and  132-m  depths,  but  the  heat  content  only  rises  slightly  compared  to  the 
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cooling  event.  At  the  C3  mooring  we  do  not  see  the  strong  loss  in  heat  content,  yet  see 
strong  high  frequency  oscillations  in  the  temperature  records  down  to  over  200-m  during 
the  time  of  strong  near-inertial  oscillations.  The  heat  content  at  both  sites  does  rise 
during  October  29-31,  a  period  of  strong  near-inertial  oscillations,  suggesting  a  flux  of 
heat  from  below  the  lowest  temperature  measurement.  The  relative  roles  that  the  initial 
wind  burst  and  subsequent  near-inertial  shears  play  in  vertical  mixing  remains  an  open 


question  and  the  subject  of  on-going  work. 
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Figure  4.6,4  lime  series  of  temperature  at  the  (a)  Cl  and  (c)  C3  moorings  before,  during  and  after 
the  wind  burst  and  near-f  oscillations  in  October,  2002  at  various  depths.  The  heat  content  (b)  for 
Cl  (black)  and  C3  (grey)  is  shown.  The  depth  of  each  temperature  record  is  labeled  on  the  right* 
Note  that  (a)-(c)  refer  to  top  to  bottom  panels,  respectively. 


4.7.  Discussion  and  Conclusions 


We  have  presented  two  years  of  sea  ice  data  which  cover  the  winter  icc  seasons  of 
2001-2003  within  Marguerite  Bay  and  on  the  wAP  shelf.  The  strong  wind  events  in 
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Marguerite  Bay  are  primarily  in  a  southwestward  direction,  and  occurred  more  frequently 
in  the  first  season. 

In  the  first  season,  we  found  that  the  mid-shelf  and  Marguerite  Bay  sites  are 
similar  until  mid-October,  2001,  when  the  ice  retreats  from  mid-shelf  and  thickens  in  the 
Bay,  In  the  second  season,  the  two  sites  within  the  Bay  have  similar  and  highly 
correlated  ice  concentrations  and  drafts.  The  primary  difference  between  the  two  sites 
begins  in  mid-October,  2002,  when  the  ice  apparently  broke  up  and  left  open  water  for  a 
number  of  days  at  Cl  but  not  C'3.  At  this  time,  the  ice  draft  dropped  dramatically  to  -1  m 
at  both  sites,  and  the  ice  concentration  also  dropped  somewhat  from  less  than  70%  to  less 
than  40%,  but  still  remained  above  60%  for  most  of  the  remaining  ice  season.  An 
interannual  comparison  at  the  B2/C3  site  shows  that  the  onset  of  ice  formation  occurred 
about  two  months  earlier  in  the  second  season.  This  is  due  to  unseasonably  cold  air 
temperatures  in  May  and  June,  2002.  the  coldest  on  record  since  meteorological 
measurements  were  started  in  1976  at  Rothera  Station  in  northern  Marguerite  Bay. 

Our  analysis  suggests  local  formation  of  ice.  The  timing  of  the  formation  of  at 
least  70%  ice  cover  is  correlated  to  the  timing  of  -9.2°  C  weather  with  a  correlation 
coefficient  of  0.81  over  the  26  years  of  historical  data  available.  Furthermore,  the  thick 
ice  drafts  suggest  ridging  and  the  importance  of  ice  dynamics  in  the  region. 

A  polynya  on  October  12,  2002  demonstrates  how  highly  compacted  ice  in  the 
Bay  breaks  apart,  allowing  easier  advection  and  strong  near- inertial  responses  to  wind 
events.  Both  years  show  significant  events  happening  in  October  involving  a  sudden 
noteworthy  reduction  in  ice  concentration  and/or  draft. 

Ice  speeds  were  typically  around  10  cm/s,  with  the  ice  typically  traveling  about  3 
km  each  day,  with  a  time-mean  net  daily  displacement  of  200-500  m.  This  indicates  that 
there  is  considerable  ice  kinetic  energy  in  all  frequency  bands  from  sub-daily  to  seasonal 
time  scales.  At  B2/C3,  the  intermittent  strong  southwestward  wind  events  act  in 
conjunction  with  a  slow  but  steady  northward  ocean  current  to  keep  ice  from  leaving  the 
area  for  most  of  the  second  winter  season  when  the  winds  were  less  energetic.  At  C 1 , 
however,  the  current  is  southward,  and  ice  moves  quickly  through  the  area. 

Wc  calculate  the  subtidal  momentum  balance  in  the  ice  and  find  that  the  temporal 
acceleration  term  is  essentially  negligible,  the  Coriolis/tilt  term  is  small,  and  the  dominant 
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balance  is  between  the  ocean  stress,  wind  stress,  and  internal  ice  stresses.  At  the 
beginning  and  end  of  the  ice  season,  when  the  ice  concentration  is  slightly  lower,  the 
ocean  stress  becomes  more  important  and  the  internal  ice  stresses  less  so.  This  has  also 
been  found  to  be  true  in  both  a  numerical  model  of  the  Arctic  Ocean  (Steele  et  al,  1997) 
and  observations  made  from  ice  buoys  in  the  Weddell  Sea  (Uotila  et  al,  2000). 

The  experimental  uncertainties  in  the  measured  variables  (e.g.,  wind  speed,  ice 
thickness,  ice  and  water  velocity)  used  to  compute  the  terms  in  the  momentum  balance 
are  generally  small,  less  than  1%,  indicating  that  the  residual  term  contains  little 
experimental  error  due  to  the  measurements.  As  a  result,  uncertainty  in  the  residual  term 
is  primarily  sensitive  to  uncertainties  in  the  parameterizations  used  to  estimate  the  air-ice 
(air-ocean)  stress  and  ice-ocean  stress,  which  can  be  significantly  larger  than  1%. 

In  order  to  investigate  the  sensitivity  of  the  low-frequency  momentum  balance  to 
the  ice-ocean  drag  coefficient,  we  varied  Cj  from  0.010  to  0.030,  a  somewhat  expanded 
range  of  measured  values  for  rough  first  year  ice  (Shirasawa  and  Ingram,  1991)  and 
contains  C<j  =  0.01 5,  the  value  used  throughout  this  thesis.  First,  recall  that  the  primary 
balance  at  all  moorings  is  between  the  wind  stress  and  the  residual  term,  with  a  minor 
contribution  from  the  ice-ocean  drag  at  times.  Therefore,  if  we  let  Cj  tend  towards  zero, 
the  balance  should  tend  towards  wind  stress  and  the  internal  ice  stress,  which  is  in  fact 
the  case  with  an  unrealistically  low  Cd  of  0.005. 

The  mean  magnitude  of  the  residual  term  only  increases  as  Cd  increases.  The 
increased  ice-ocean  stress  is  balanced  not  only  by  the  wind  stress,  but  also  by  the  residual 
term.  Upon  inspection  of  the  vectors  at  various  times,  this  can  be  attributed  to  the  fact 
that  the  ice-ocean  stress  does  not  purely  oppose  the  wind  stress,  which  would  be  the  case 
in  a  motionless  (or  very  slow  moving)  ocean  beneath  wind-forced  ice.  Furthermore, 
there  is  no  value  of  Cd  that  creates  a  balance  between  the  wind  and  ocean  stress  and 
eliminates  the  residual  term.  The  residual  term  not  only  balances  the  wind  forcing,  but 
also  the  ocean  forcing.  It  should  be  noted  that  there  is  no  physical  reason  to  believe  Cd  is 
even  constant  in  time,  since  the  ice  could  become  rougher  or  smoother  during  the  8- 
nionth  ice  season. 

The  calculated  residual  term  is  indeed  sensitive  to  the  chosen  value  of  Cd,  with  the 
mean  of  the  absolute  value  of  the  residual  increasing  significantly  with  a  doubling  of  the 
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drag  coefficient.  If  one  were  to  have  any  reason  to  believe  that  the  ice-ocean  drag  were 
higher,  then  the  balance  would  be  a  larger  residual  term  balancing  both  the  wind  and  the 
ice-ocean  stresses.  For  these  reasons,  we  conclude  that  while  the  stress  parameterizations 
used  here  contribute  some  uncertainty  in  the  momentum  balance,  the  interpretation  of  the 
residual  term  as  internal  ice  stress  is  robust. 

We  observed  strong  semi-diurnal  oscillations  in  both  the  sea  ice  and  upper  ocean. 
Power  density  spectra  and  rotary  spectra  show  that  these  are  not  tidal  but  near-inertial 
oscillations.  An  analysis  of  individual  wind  events  elucidates  the  coupled  ice-ocean 
response  in  this  energetic  band.  A  typical  example  at  the  end  of  October,  2002  shows  a 
pulse  of  wind  stress  which  spins  up  near-inertial  oscillations  in  both  the  sea  ice  and  upper 
ocean.  These  oscillations  persist  for  3  days,  until  an  increase  in  internal  ice  stress  forces 
northward  motion  in  the  ice.  The  cause  of  the  strong  internal  ice  stresses  shows  the 
importance  of  non-local  processes,  which  are  unknown  in  this  case. 

Future  work  will  investigate  the  ocean  mixing  caused  by  these  wind  and  near- 
inertial  events.  Specifically,  how  can  these  events  move  heat  and  nutrients  up  from 
deeper  water?  The  periodic  opening  and  closing  of  leads  will  cause  ice  production  and 
thus  convection.  The  stress  at  the  ice-ocean  interface  creates  turbulent  kinetic  energy 
which  also  acts  to  mix  heat  up  from  below.  Finally,  near-inertial  shear  at  the  base  of  the 
mixed  layer  could  also  cause  vertical  mixing  and  is  now  thought  to  be  the  primary  source 
of  small-scale  mixing  on  the  wAP  shelf  (Howard  et  al,  2004). 
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5.  Chapter  5 


One-dimensional  modeling  of  vertical  mixing  in  a 
seasonally  ice-covered,  stratified,  high-latitude  coastal 
ocean 

5.1.  Introduction 

Air- sea  heat  flux  can  be  one  to  two  orders  of  magnitude  larger  over  open  water 
than  over  thick  (~  l  m)  sea  ice.  Therefore,  the  exchanges  across  small  areas  of  open 
water  or  thin  ice  in  the  icepack  (as  little  as  1%)  can  dominate  the  overall  air-sea  heat  flux 
budget.  High  frequency  motions,  such  as  near-inertial  oscillations  or  tidal  motions,  can 
lead  to  differential  ice  motion.  This  can  open  and  close  leads  on  time  scales  shorter  than 
the  thermodynamic  time  scale  to  freeze  them  solidly  enough  to  provide  sufficient  icc 
strength  to  inhibit  relative  ice  motion.  Consequently,  the  new  ice  formed  in  leads  is 
crushed  and  deposited  around  floe  edges.  Inclusion  of  moderately  high  frequency 
(subdaily)  deformation  dynamics  in  the  forcing  of  a  thermodynamic  ice  model  results  in 
an  increase  in  ice  thickness  of  about  13%  over  a  5-month  modet  simulation  (Heil  et  al., 
2001 ).  Including  an  integrated  momentum  equation  for  the  ice  and  upper  ocean  rather 
than  the  ice  alone,  or  “inertial  embedding”  (McPhee,  1978),  into  the  dynamics  of  the  ice 
increased  winter  ice  mass  by  about  20%  (Heil  and  Hibler,  2002).  In  addition  to  its  role  in 
augmenting  ice  formation,  high  frequency  ice  motion  and  lead  formation  will  enhance 
oceanic  mixing  by  near-inertial  shear  at  the  base  of  the  mixed  layer  and  convection 
during  ice  production  and  brine  rejection. 

As  shown  in  Chapter  4,  power  density  spectra  of  ice  and  ocean  velocities  at  the 
Marguerite  Bay  moorings  show  peaks  at  the  inertial  frequency,  with  more  energy  in  this 
band  in  the  ice  than  in  the  water.  These  energetic  motions  have  rms  amplitudes  up  to  -30 
cm/s.  Most  of  the  periods  of  high  kinetic  energy  at  the  inertial  frequency  occurred  during 
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periods  of  lower  ice  concentration  (Figure  5.1.1),  supporting  the  notion  that  a  looser  ice 
pack  can  sustain  near-inertial  oscillations.  This  motivates  the  investigation  of  oceanic 
near-inertial  motions  in  the  presence  of  sea  ice. 

(a) 

CCW  demod  for  ice  at  c3  red  observed,  b!ue  from  demod 


(b) 
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Figure  5.1.1  (a)  Upper  panel  shows  ice  kinetic  energy  (U!)  observed  (red)  and  that  accounted  for  by 
near-inertial  oscillations  (blue)  as  calculated  with  a  complex  demodulation.  The  lower  panel  (h) 
shows  observed  ice  concentrations  at  the  C3  mooring.  Similar  results  were  found  at  the  nearby  Cl 
mooring  site  (not  shown). 


Deep-keeled  icebergs  precluded  moored  instrumentation  in  the  upper  100  m  of 
ocean  in  the  SO  GLOBEC  field  program.  Therefore  we  arc  unable  to  observe  important 
mixing  processes  with  the  moorings.  For  example,  we  cannot  resolve  the  depth  of  the 
mixed  layer.  The  long-wave  radiative  cooling  has  been  shown  to  be  an  important  term  in 
the  net  air-sea  heat  flux  (Beardsley,  ct  al.,  2006).  The  radiative  heat  fluxes  were 
measured  in-situ  only  with  shipboard  instruments,  not  on  the  Automatic  Weather 
Stations,  and  there  were  no  ships  in  this  area  during  the  spring  break-up  when  strong  ice- 
covered  near-inertial  oscillations  in  both  the  sea  ice  and  upper  ocean  were  observed. 
Furthermore,  CTD  casts  and  hull-mounted  ADCPs  provide  important  initial  conditions 
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for  any  mixing  study.  This  prompts  the  use  of  an  one -dimensional  surface  mixed  layer 
model  to  study  the  relative  importance  of  vertical  mixing  processes  in  an  idealized 
seasonally  ice-covered,  stratified,  high-latitude  coastal  ocean. 

Here  I  examine  two  potentially  important  ocean  mixing  processes,  namely  shear 
associated  with  a  wind  event  and  subsequent  near-inertial  oscillations  and  convection  due 
to  brine  rejection  associated  with  freezing  as  leads  open  and  close  due  to  the  near-inertial 
motions  (Figure  3.1.1).  Shear  associated  with  near-inertial  oscillations  has  been  identified 
as  a  potentially  important  process,  especially  when  compared  with  double  diffusion 
(Howard  ct  al.,  2004). 
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Figure  5.1.2  Schematic  of  mixing  processes  considered  in  the  PYVP  model  with  ice. 


It  is  important  to  examine  these  processes  individually.  Therefore  the  analysis 
begins  with  a  base  ease  with  no  ice  and  a  wind  event  which  generates  inertial  oscillations. 
Then  the  consequences  of  an  ice  cover  arc  explored,  as  a  function  of  both  the  thickness 
and  roughness  (drag  coefficient)  oflhc  ice.  The  initial  and  surface  boundary  conditions 
arc  idealized  based  on  the  western  Antarctic  Peninsula  (wAP)  shelf  environment. 
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Subsequently,  an  air-sea  heat  flux  modulated  at  the  inertial  frequency  is  included 
to  simulate  the  phenomenon  of  leads  opening  and  closing  with  divergent  high  frequency 
ice  motion.  The  feedback  associated  with  the  upward  flux  of  heat  and  ice  production  is 
discussed,  as  well  as  the  slowing  of  the  mixed  layer  due  to  entrainment  of  deeper 
motionless  fluid.  Additional  effects  of  rotation  of  the  wind  stress  vector  and  reflection  of 
a  barotropic  wave  from  the  coastal  boundary  arc  explored.  Finally,  the  implications  of 
the  results  on  mixing  in  Marguerite  Bay  arc  discussed. 

5.2.  The  Model:  PWP  with  ice. 


The  mixing  model  employed  is  based  on  that  of  Price,  Weller  and  Pinkcl  ( 1 986) 
(PWP)  with  some  modifications  made  as  noted  below, 

5.2. 1.  The  basic  PWP  model 


The  basic  PWrP  vertical  mixing  model  allows  for  mixing  due  to  static  instability 
(density  inversions)  and  shear  instability,  and  allows  for  free  inertial  oscillations.  The 
basic  horizontal  momentum,  salt  and  heat  balance  equations  arc  as  follows  (with  positive 
z  downward  and  downward  surface  wind  stress,  heat  and  salt  fluxes  positive  into  the 
ocean): 
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with  the  following  stability  criteria: 
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The  model  is  initialized  with  specified  temperature  ( T ),  salinity  ( S )  and  velocity 
( u )  profiles.  Velocities  rotate  freely  at  the  inertial  frequency  of  /  (equation  5.1),  which 
equals  2 n  i  12.74  hrs  at  70°  S.  Prescribed  time  series  of  surface  fluxes  of  momentum 
( r ),  heat  ( Qr )  and  salt  ( Qs )  arc  applied  as  the  surface  flux  boundary  conditions  for 
equations  5. 1 ,  5.2  and  5.3. 

At  each  time  step,  the  model  first  applies  the  temperature  and  salinity  surface 
boundary  conditions  (described  in  detail  in  section  5.2.3)  in  the  top  depth  cell.  Then  the 
density  profile  is  checked  for  static  stability  starting  from  the  surface  cell  downward 
(Equation  5.4)  and  the  water  column  is  mixed,  including  momentum,  until  stability  is 
met.  This  creates  a  homogenized  surface  layer,  or  slab,  of  depth  h .  Then  the  mode! 
applies  the  surface  wind  stress  by  adding  the  momentum  to  the  “surface  mixed  layer”, 
defined  by  a  change  in  density,  Ap ,  of  5  x  10"4  kg/m1  from  the  surface  value.  One 
consequence  of  this  definition  is  that  the  “surface  mixed  layer”  may  contain  some  weak 
stratification  and  shear.  The  bulk  Richardson  number  Rih  stability  criterion  (equation 
5.5)  is  then  checked  and  cells  arc  entrained  into  the  slab  layer  until  the  profile  is  stable. 
The  differences  A p  and  Ah  arc  taken  between  the  slab  layer  and  the  level  just  beneath. 
This  simulates  entrainment  into  the  mixed  layer  due  to  shear  instability. 

The  first  two  criteria  (equations  5.4  and  5.5)  often  produce  a  sharp  jump  in  u  and 
p  at  the  base  of  the  slab  layer.  The  final  mixing  criterion  (equation  5.6)  uses  a  gradient 
Richardson  number  RiK  stability  criterion  to  simulate  mixing  due  to  small-scale  shear 
instability.  A  profile  of  Ri  is  made  and  any  adjacent  cells  which  do  not  meet  the 

stability  criterion  have  their  temperature,  salinity  and  momentum  mixed  together.  Since 
this  generally  produces  other  cell  pairs  which  do  not  meet  the  RiK  stability  criterion,  the 

process  iterates  until  the  entire  profile  is  stable.  This  process  bears  some  resemblance  to 
a  diffusive  mixing  process  in  that  it  spatially  smoothes  the  profile.  The  region  between 
the  base  of  the  slab  mixed  layer  and  the  depth  to  which  the  shear  mixing  penetrates  is 
known  as  the  transition  layer  (Schudlich  and  Price,  1 992)  (Figure  5.2. 1 ). 
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>0.25 


(5.6) 
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Figure  5*2  J  Illustration  of  the  defined  slab  layer,  model-defined  misted  layer  and  depth  of  maximum 
mixing  and  (zoomed  in  inset)  the  sequence  of  the  model's  mixing  scheme  at  t=  1 3  hrs.  The  model 
“mixed  layer"  contains  a  slab  layer  over  a  region  of  very  weak  (A p  of  5  x  HT*  kg/nrr)  stratification. 

The  inset  shows  salinity  profiles  at  intermediate  stages  in  the  mixing  scheme.  The  green  line  denotes 
the  initial  profile  in  the  time  step,  with  added  brine  in  the  surface  cell.  The  red  lane  indicates  the 
profile  after  this  added  brine  has  been  mixed  down  far  enough  to  satisfy  the  static  stability  and  Rih 
criterion  (equations  5.4  and  5.5),  producing  a  slab-like  mixed  layer  with  a  sharp  gradient  at  its  base. 
The  blue  line  shows  the  final  salinity  profile,  after  shear  flow  stability  (equation  5.6)  has  been 
achieved  hv  fluxing  salt,  heat  and  momentum  upwards  and  smoothing  the  gradient.  Note  that  the 
water  column  is  salt  stratified,  with  a  weak  destabilizing  thermal  gradient  (not  show  n). 


The  density  p  is  calculated  from  the  model  temperature  and  salinity  using  the 
fully  nonlinear  equation  of  state  (Fofonoff  and  Millard,  1 983).  There  is  a  completely 
mixed  layer  resulting  from  criteria  5.4  and  5.5  which  is  treated  as  a  slab  layer  wrilh 
infinite  vertical  diffusivity  and  viscosity.  After  the  gradient  Richardson  number  criterion 
(equation  5.6)  is  applied,  a  final  “surface  mixed  layer"  is  calculated  for  use  in  plotting  the 
results.  This  final  “surface  mixed  layer"  is  defined  by  a  Apof  5  x  10^  kg/m’  from  the 
surface  value,  identical  to  the  definition  used  for  adding  momentum  from  wind  stress. 
This  mixed  layer  is  technically  not  a  fully  mixed  slab  layer,  so  the  actual  slab  layer 
thickness  is  also  shown  in  many  figures.  The  choice  of  A p  (within  reasonable  bounds) 
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has  a  negligible  effect  on  the  modeled  stratification  and  velocity.  The  simulations  were 
run  at  high  resolution,  with  a  grid  size,  Ar ,  of  0.5  m,  a  A/ of  120  s  and  zero  background 
vertical  diffusivity  for  scalars  and  momentum. 

Vertical  fluxes  of  heat  and  salt  arc  calculated  after  the  simulation  has  been 
completed  using  equations  5.7  and  5.8,  which  integrate  from  the  bottom  of  the  domain  up 
to  the  surface  and  arc  identical  in  form  to  equation  6.1  in  PWP: 


<5-7> 

<5-8> 

The  water  column  in  the  Marguerite  Bay  region  is  salt  stratified.  Using  the 
definitions  in  McDougall  (1987),  the  change  in  density  due  to  salt  /3A S  is  one  order  of 

magnitude  larger  than  that  due  to  temperature  aAT ,  taking  the  differences  in  S  and 
T  between  the  top  and  bottom  of  the  domain.  Therefore,  vertical  fluxes  of  heat  have  a 
much  smaller  impact  on  the  stratification  than  those  of  salt.  This  is  discussed  further  in 
section  5.3.4. 


5.2.2.  PWP  with  ice 


Simplified  ice  dynamics  were  added  to  the  PWP  model  to  investigate  the  role  of 
ice  in  vertical  mixing  on  a  stratified  continental  shelf.  Unlike  standard  PWP  applications. 


the  wind  first  acts  on  the  ice,  with  the  air-ice  stress  r„.  calculated  using  the  velocity 

difference  between  the  two  and,  for  simplicity,  assuming  neutral  stability  (Large  and 
Pond,  1981).  The  ice  motion  is  modeled  as  free  drift  and  allows  for  inertial  oscillations: 


dul 

~St~ 


(5.9) 


It  is  assumed  that  the  ice  covers  a  large  enough  percentage  of  the  sea  surface  that 
the  wind  stress  acts  only  on  the  ice.  The  ice-ocean  stress  rw  is  calculated  using  a 
quadratic  drag  law: 
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(5.10) 


u  -u 


with  Cj=15  x  10  3  (Shirasawa  and  Ingram,  1991 )  (unless  otherwise  noted),  p„thc  density 
of  seawater  ( 1 024  kg/m  )  and  pt  the  density  of  sea  ice  (920  kg/nr).  The  ice-ocean  stress 


zia  is  used  as  the  surface  ocean  boundary  condition  r  in  the  PWP  ocean  model 
(equation  5.1). 

The  latent  and  sensible  heal  fluxes  arc  calculated  by  applying  bulk  formula 
(Fairall  ct  al.,  1996)  to  the  simulated  sea  surface  temperature  and  salinity  and  a  constant 
air  temperature  (-20  °C)  and  relative  humidity  (70%)  based  on  observations  in  Marguerite 
Bay  (Beardsley  ct  ah,  2006).  After  the  wind  event,  a  background  wind  of  5  m/s  is  used  in 
the  bulk  formula  to  calculate  heat  fluxes,  based  on  the  observed  light  and  variable  winds 
in  Marguerite  Bay,  but  is  not  used  in  the  calculation  of  the  surface  stress  to  simplify 
interpretation  of  the  model  output.  The  radiative  heat  fluxes  arc  assumed  constant,  with 
values  based  on  shipboard  measurements.  The  net  radiative  iicat  Dux  is  set  at  100  W/m 2 
(loss),  under  the  assumption  that  the  shortwave  heat  flux  is  small  during  the  period  of 
time  of  interest,  and  the  flux  is  primarily  long-wave  heat  loss.  Assuming  no  heat  flux 
through  the  ice,  the  total  heat  flux  is  then  sealed  by  the  open  water  fraction  A  ,  the 
calculation  of  which  is  described  in  section  5.3.3. 

The  ice  model  docs  not  contain  dynamics  such  as  ridging  and  breaking,  a 
horizontal ly-lwo-dimcnsional  problem  impossible  to  include  into  the  one-dimensional 
PWP  model.  The  phenomenon  oflcads  opening  and  closing  due  to  convergence  and 
divergence  of  the  ice  drift  is  incorporated  by  allowing  the  open  water  fraction  A  to  be 
time-dependent. 

The  simple  thermodynamic  ice  model  is  achieved  by  creating  ice  when  the 
ocean's  surface  temperature  Tiurf  drops  below  the  freezing  point  of  seawater  at  the 

surface  salinity  S  .  in  the  simulation.  The  excess  heat  is  removed  by  freezing  and  ice  is 
created  at  a  rate  of  Ft  at  a  salinity  of  Sice  which  is  set  at  4  PSU  (Perovich,  2004) 

(equation  5. 1 1 ).  The  basic  equations  can  be  expressed  as 

p,LlF,=p.crT*''J"'  te  (5.11) 
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(5.12) 


—  =  ^(5^,-5*,)/ Az 

A  t 

with  cp  the  specific  heat  of  seawater  (4183.3  J/kg/°C)  and  Lf  is  the  latent  heat  of  fusion 

of  sea  ice  (333  kJ/kg,  Wadhams,  2000)).  F(  is  then  used  to  calculate  the  salt  flux  into  the 

ocean  at  the  surface  (equation  5.12).  This  model  ignores  changes  in  ice  temperature  and 
the  melting  of  ice  and  assumes  that  thin  ice  in  the  leads  has  no  effect  on  the  air-sea  heal 
flux,  which  is  reasonable  since  a  maximum  of~I  cm  of  ice  is  formed  in  the  model  during 
an  inertial  period.  Note  that  more  elaborate  thermodynamic  ice  models  exist  (c.g. 

Maykut  and  Untcrstcincr,  1 97 1 ),  however  we  believe  the  simple  ice  production  model 
used  in  this  Chapter  captures  the  essential  feature  of  ice  production  on  short  time  scales 
and  subsequent  brine  rejection  and  that  iccmclt  in  these  situations  (especially  the  cold 
upper  ocean)  would  be  negligible  (M.  McPhcc,  pres.  comm.). 

The  air-sea  heat  flux  allows  heat  to  vent  out  from  the  ocean.  This  venting  is 
thought  to  be  important  as  heat  brought  up  from  deeper  waters  by  mixing,  which  would 
otherwise  melt  ice  and  form  a  thin  fresh  surface  layer  which  would  isolate  the  wind  stress 
from  the  deeper  waters  (McPhcc,  pers.  comm.),  is  instead  released  to  the  atmosphere. 

The  crushing  of  newly  formed  thin  ice  in  the  leads  is  assumed  to  have  no  effect  on  the 
ice’s  momentum  (equation  5.9). 

The  model  is  thought  to  be  a  spatial  average  over  an  area  consisting  of  both  open 
water  and  ice.  Essentially,  horizontal  diffusivity  over  this  area  is  thought  to  be  infinite. 
No  plume  dynamics  arc  included,  but  could  be  important.  Plumes  could  potentially  not 
mix  completely  as  they  carry  denser  water  downwards  and  then  overshoot  the  level 
predicted  by  the  PWP  model,  causing  further  mixing.  Furthermore,  a  spatial  pattern  of 
convccting  plumes  beneath  leads  could  create  a  pattern  of  internal  waves  on  the 
pycnocltnc,  heaving  upwards  between  leads  and  potentially  causing  feedbacks  by  melting 
ice  (M.  Levine,  pers.  comm.).  This  problem  is  inherently  three-dimensional  and  beyond 
the  scope  of  this  study. 

5.2.3.  Initial  and  Boundary  Conditions 

All  model  runs  were  initialized  with  an  idealized  stratification  which  consisted  of 
a  uniform  65-m  deep  mixed  layer  overlying  a  linear  stratification  down  to  100  m  (Figure 
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5.2.2a).  The  domain  was  limited  to  the  upper  100  m  because  the  maximum  depth  of 
mixing  (defined  as  the  deepest  point  which  has  been  influenced  by  the  mixing  model)  did 
not  penetrate  below  that  level  in  any  of  the  simulations.  The  stratification  was  chosen  to 
mimic  a  CTD  cast  taken  on  July  29,  2001  in  the  region,  at  66°  15.27’  S,  70°  20.81’  W 
(Figure  5.2.2b).  This  cast  was  chosen  as  it  contains  the  essential  features  characterizing 
the  western  Antarctic  Peninsula  shelf  regions,  namely  a  deep  (~65  m)  mixed  layer  of  cold 
(-1.83  °C),  fresh  (33.91  PSU)  Antarctic  Surface  Water  overlaying  a  region  of  sharp 
stratification  and  a  lower  layer  of  Upper  Circumpolar  Deep  Water.  A  linearized  version 
of  the  stratification  was  used  to  simplify  and  clarify  the  interpretation  of  the  model  results 
and  remove  the  small-scale  vertical  structure  in  the  observed  profile.  It  is  noted  that  the 
details  of  the  model  results  arc  sensitive  to  the  initial  temperature  and  salinity  profiles, 
but  the  general  conclusions  regarding  mixing  arc  robust.  The  initial  water  and  ice 
velocities  arc  zero,  unless  otherwise  noted. 
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Figure  5.2.2  (a)  Initial  model  temperature  (solid  black),  salinity  (solid  red)  and  potential  density 
(solid  green)  profiles  used  in  the  simulations  overlaying  the  observed  values  (dotted  lines),  (b)  The 
full  CTD  cast. 

The  model’s  surface  wind  was  based  on  typical  events  observed  on  the  wAP  shell' 
(Beardsley  ct  al.,  2006):  a  unidirectional  (eastward)  wind  event  lasting  24  hrs  with  a 
Gaussian  shape  and  a  peak  intensity  varying  between  5  and  25  m/s,  depending  on  the 
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simulation  (15  m/s  in  Figure  5.2.3).  As  noted  earlier,  the  wind  stress  on  the  ice  was 
calculated  using  the  neutral  quadratic  drag  formula  of  Large  and  Pond  (1981).  The  ice- 
ocean  stress  was  calculated  using  a  quadratic  drag  law  (equation  5. 10)  assuming  rough 
first-year  ice  with  the  ice-ocean  drag  coefficient  set  at  15  x  10  '  unless  otherwise  noted. 
Additional  simulations  were  done  with  a  rotating  wind,  as  well  as  with  non-zero  air-sea 
heat  fluxes. 
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Figure  5.2.3  Wind  forcing  for  idealized  runs,  with  this  example  having  a  maximum  wind  of  15  m/s. 

5.3.  Experiments 

5.3.1.  Ice  versus  no  ice 

An  initial  set  of  experiments  compared  the  ocean  response  to  a  wind  event  under 
ice-free  conditions  to  those  with  total  ice  cover.  Both  surface  heat  and  salt  fluxes  were  set 
to  zero  and  thus  there  was  no  static  instability  in  the  water  column. 

The  initial  wind  event  causes  all  of  the  observed  mixing,  with  subsequent 
persistent  inertial  oscillations  stable  to  shear  instability.  The  presence  of  ice  makes 
neither  a  large  difference  in  the  modeled  response  of  temperature,  salinity  and  velocity  to 
the  wind  event  (Figures  5.3.1, 5.3.2  and  5.3.3),  nor  has  a  great  impact  on  the  overall 
structure  of  the  vertical  heat  and  salt  fluxes  (Figures  5.3.4  and  5.3.5).  While  these  results 
are  for  a  wind  event  with  a  maximum  velocity  of  15  m/s  (Figure  5.2.3),  ice  had  little 
impact  on  the  model  results  for  wind  events  ranging  in  strength  between  5  and  25  m/s. 
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Figure  5.3.1  Depth-time  contour  plots  of  the  temperature  response  to  a  wind  event  with  a  maximum 
speed  ofl5  m/s  and  (a)  no  ice  and  (b)  1  m  ice.  White  line  denotes  the  model-defined  mixed  layer 
depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of  the 
slab  layer. 
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Figure  5.3.2  Depth-time  contour  plots  of  the  salinity  response  to  a  wind  event  with  a  maximum  speed 
of  15  m/s  and  (a)  no  icc  and  (h)  1  m  ice.  White  Line  denotes  the  model-defined  mixed  layer  depth,  the 
red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of  the  slab  layer. 


122 


(a) 


Max  Wirid  =  1 5  m/s  Ice  =  0  fn 


10  20  30  4Q  50  60  70 

Time  thrsj 


(b) 


Max  Wind  *  1 5  m/s.  Ice  ■  I  m 


10  20  30  40  50  60  70 

Time  sh*s> 


Figure  5.3*3  Depth-time  contour  plots  of  the  velocity  response  to  a  wind  event  with  a  maximum  speed 
of  15  m/s  and  (a)  no  ice  and  (b)  1  m  ice*  White  line  denotes  the  model-defined  mixed  layer  depth,  the 
red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of  the  slab  layer* 


123 


(a) 


Max  Wind  *  15  m/s  Ice  *  0  m 


Time  ihrsi 


(b) 


Max  Wind  =  IS  m/s  ice  a  t  m 


Tim#  (hfsi 


Figure  5 .3,4  Pepth-timc  contour  plots  of  the  vertical  heat  flux  response  to  a  wind  event  with  a 
maximum  speed  of  15  m/s  and  (a)  no  ice  and  (b)  1  m  ice.  White  line  denotes  the  model-defined  mixed 
layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of 
the  slah  layer. 
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Figure  5.3.5  Depth-time  contour  plots  of  the  vertical  salt  flux  response  to  a  wind  event  with  a 
maximum  speed  of  15  m/s  and  (a)  no  ice  and  (b)  l  m  ice*  White  line  denotes  the  model-defined  mixed 
layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  hasc  of 
the  slab  layer. 
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Varying  the  ice  thickness  between  0  and  2  m  had  little  effect  on  neither  the  final 
stratification  (Figures  5.3.6,  5.3.7  and  5.3.8)  nor  inertial  oscillation  speed  (Figure  5.3.9) 
for  a  wind  burst  as  shown  in  Figure  5.2.3.  The  small  effect  of  ice  thickness  is  due  only  to 
the  added  mass  of  ice,  a  relatively  small  effect  when  compared  to  the  total  mass  of  the 
icc/mixcd  layer  system  (-  1 .5%  for  1  m  ice). 


Effect  of  tee,  Max  Wind=1 5 


Figure  5.3.6  The  effect  of  various  ice  thicknesses  on  the  temperature  profile.  The  insets  are  zoomed 
in  on  the  change  in  the  mixed  layer  depth  (left)  and  maximum  depth  of  mixing  (right). 
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Effect  of  Ice,  Max  W»nd=15 
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Figure  5.3.7  The  effect  of  various  ice  thicknesses  on  the  salinity  profile.  The  insets  are  zoomed  in  on 
the  change  in  the  mixed  layer  depth  (left)  and  maximum  depth  of  mixing  (right). 


Effect  of  Ice,  Max  Wmd=15 


Figure  53.8  The  effect  of  various  ice  thicknesses  on  the  density  profile.  The  insets  arc  zoomed  in  on 
(he  change  in  the  mixed  layer  depth  (left)  and  maximum  depth  of  mixing  (right). 
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Effect  of  Ice,  Max  Wmd=15 
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Figure  53,9  The  effect  of  various  ice  thicknesses  on  the  current  speed  profile. 

The  observation  that  the  initial  wind  event  causes  all  of  the  observed  mixing  is  consistent 
with  Plucddcmann  and  Farrar  (2006),  who  examined  a  variety  of  icc-frcc  datasets  and 
found  that  the  majority  of  strong  mixing  events  arc  very  short,  less  than  1/10  of  an 
inertial  period,  with  the  longest  events  being  at  most  1/2  an  inertial  period. 


5.3,2.  Sensitivity  to  the  ice-ocean  drag  coefficient,  Cd 

The  next  set  of  simulations  explored  the  effect  of  varying  the  icc-occan  drag 

coefficient,  Cj.  To  do  this,  Cj  was  varied  from  0.000 1  to  0.030,  the  range  of  values 
estimated  through  observational  methods  such  as  eddy  correlation  in  a  variety  of  regions 
in  both  the  Northern  and  Southern  Hemispheres  (Shirasawa  and  Ingram,  1991).  The 
large  range  in  values  reflects  variations  in  the  roughness  of  the  bottom  of  the  ice  (Iocs. 
The  maximum  wind  speed  was  15  m/s  and  the  ice  thickness  was  1  m  for  all  runs  in  this 
section. 

For  the  set  of  initial  and  boundary  conditions  tested  here,  changes  in  Cj  had  little 
effect  on  the  simulation  for  values  greater  than  0.00 1  The  ice  velocity  becomes  sensitive 
to  the  choice  of  Cd  with  values  lower  than  0.001 ,  even  having  a  different  final  phase  of 
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ihc  inertial  oscillations  compared  with  runs  at  higher  values  (Figure  5.3. 10a).  In 
particular,  the  distance  covered  by  the  ice  during  the  initial  wind  burst  is  greatly 
increased  for  the  very  low  values  of  Cj.  One  should  note,  however,  that  observational ly 
based  values  of  C<j  less  than  0.001  have  not  been  observed  for  first-year  drifting  ice 
(Shirisawa  and  Ingram,  1991 )  such  as  that  found  in  Marguerite  Bay,  which  is  relatively 
rough  (Pcrovich  ct  at.,  2004). 

The  response  of  the  upper  ocean  mixed  layer  current  to  changes  in  the  drag 
coefficient  is  small  (Figure  5.3.10b),  w'ith  all  flows  simulated  with  a  Q  >  0.0005  within 
+/*  3%  of  the  open  water  value.  This  insensitivity  is  due  to  low  values  of  the  ice-ocean 
quadratic  drag  smoothing  the  ice-ocean  stress  in  time,  thus  weakening  the  final 
magnitude  of  the  inertial  oscillations.  More  abrupt  changes  in  surface  stress  lead  to 
stronger  inertial  oscillations.  Cases  with  low  Cd  have  lower  ice-ocean  stress  and 
therefore  allow  higher  ice  velocities  and  lower  ocean  velocities.  This  creates  a  smoothing 
in  time  as  the  ocean  takes  longer  to  match  velocity  with  the  ice. 

Changes  in  the  ice-ocean  drag  coefficient  act  to  distribute  the  momentum 
differently,  but  the  total  momentum  of  the  system  remains  largely  unchanged.  Any  small 
changes  in  the  total  momentum  come  from  the  dependence  of  the  air-ice  stress  on  the  ice 
velocity,  so  when  the  ice  moves  faster  (lower  Cj)  and  the  velocity  difference  between  the 
ice  and  air  is  smaller,  the  total  imparted  momentum  is  decreased. 
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Figure  5,3,10  Northward  (solid  lines)  and  eastward  (dashed  lines)  (a)  ice  velocity  and  (b)  ocean  mixed 
layer  velocity  (m/s)  for  runs  with  different  ice-ocean  drag  coefficients.  Values  in  the  legend  indicate 
the  ice-ocean  drag*l0\  with  4fre’  in  (b)  referring  to  the  open  water  case. 
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Examination  of  the  final  temperature  profile’s  response  to  changes  in  the  drag 
coefficient  reveals  deeper  maximum  depth  of  mixing  and  shallower  mixed  layer  as  Cd  is 
increased  (Figure  5.3. 1 1 ).  Again,  for  reasonable  values  of  Cd  (greater  than  0.00 1 ),  the 
response  docs  not  vary  greatly  between  runs. 

One  conclusion  from  these  experiments  is  that  as  Cd  increases,  the  response  in  the 
model  tends  towards  the  open  water  ease  and  maximum  mixing  for  a  given  wind  event. 
Furthermore,  the  modeled  response  of  the  upper  ocean  is  not  sensitive  to  the  range  of 
values  reasonable  for  the  rough  deformed  first-year  ice  found  on  the  wAP  shelf  (Pcrovich 
ct  al.,  2004),  i.c.  Cd>0.001  (Shirisawa  and  Ingram,  1991). 


Ice-Ocean  Drag  Coefficient  *  1e-3  Effects,  Temp 


Figure  5.3.1 1  Final  temperature  profile  of  runs  with  different  values  of  Cd.  Values  in  the  legend 
indicate  the  ice-ocean  drag*  1(>\  with  ‘Ini1  referring  to  the  initial  profile. 


5.3.3.  Lead  opening  and  closing  and  air-sea  heat  fluxes 

A  second  set  of  simulations  was  performed  in  which  the  effects  of  ice  leads  were 
examined.  As  presented  earlier,  this  was  accomplished  by  allowing  a  surface  air-sea  heat 
flux  through  leads  that  open  and  close  due  to  inertial  divergent  motions  in  the  ice  pack. 
The  hypothesis  being  tested  is  that  periodic  ice  formation  and  subsequent  brine  rejection 
and  static  instability  and  mixing  will  influence  the  stratification  and  result  in  enhanced 
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shear  instability  during  and  after  the  initial  wind  burst.  All  runs  in  this  section  had  1  -m 
thick  tee. 

As  described  above,  the  surface  heat  flux  is  modulated  by  the  open  water  fraction 
A  ,  the  fraction  of  the  surface  area  that  is  ice-free  over  a  spatial  average.  The  open  water 
fraction  is  set  ad  hoc  proportional  to  the  eastward  ice  motion  u,(t) , 


A -A. 


2U; 


with  /lm!K  depending  on  the  individual  simulation  and  u( , 


taken  from  the  ease  with  no  lead  opening  and  closing  but  the  same  wind  forcing  (Figure 
5.3.12a).  An  alternative  method  for  calculating  A  could  be  to  use  *,(/)  instead  of  m,(0  • 
However,  integrating  «/,(/)  in  time  to  get  x((/)  simply  changes  the  phase  of  the  periodic 
signal  by  n! 4 .  Since  these  simulations  arc  idealized  such  that  the  north-south  and  east- 
west  directions  arc  interchangeable,  the  use  of  t/,(r)  or  .*■,(?)  in  calculating  A  docs  not 
change  the  results. 

When  the  surface  water  is  cooled  below  freezing,  ice  is  formed  and  the  equivalent 
amount  of  brine  (equations  5.1 1  and  5. 12)  is  injected  into  the  top  layer  of  the  ocean 
model  causing  static  instability  and  mixing.  The  simulated  sea  surface  temperature  and 
salinity  and  a  constant  air  temperature  (-20  °C)  and  relative  humidity  (70%)  arc  input  into 
bulk  formulae  (Fairall  ct  al.,  1996)  to  calculate  gross  latent  and  sensible  heat  fluxes 
(Figure  5.3.12b  blue  and  green  lines)  which  arc  then  scaled  by  the  open  water  fraction  A  . 
In  order  to  simulate  the  observed  light  and  variable  winds  in  Marguerite  Bay  (Beardsley 
ct  al,  2006),  a  background  wind  of  5  m/s  is  used  in  the  bulk  formulae  to  calculate  heat 
fluxes  after  the  initial  wind  event  but  is  not  used  in  the  calculation  of  the  surface  stress. 
This  causes  the  non-zero  sensible  and  latent  fluxes  after  the  initial  24  hours.  The 
radiative  heat  fluxes  arc  assumed  constant  at  100  W/m:  (loss)  (Figure  5.3. 1 2b  red  line) 
before  applying  the  open  water  fraction  A  .  The  total  heat  flux  (sum  of  sensible,  latent 
and  radiative)  is  then  sealed  by  the  open  water  fraction  A  to  get  the  net  heat  flux  applied 
in  the  model  (Figure  5.3.12b  black  line).  For  example,  with  the  maximum  open  water 
Amm  set  to  0.5  and  a  maximum  wind  of  15  m/s,  the  net  surface  cooling  oscillates  between 
20  and  180  W/m2. 
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Figure  5.3.12  Time  series  of  (a)  applied  open  water  fraction  A  with  a  prescribed  maximum  Amja  of 

0.5  and  (b)  resulting  air-sea  heat  fluxes  consisting  of  sensible  (blue),  latent  (green),  net  radial  ive  (red) 
components.  The  ncl  (black)  heat  flux  is  the  sum  of  the  sensible,  latent  and  radiative  fluxes 
multiplied  by  the  open  water  fraction  A  . 

The  active  air-sea  heat  flux  not  only  allows  for  icc  production  and  destabilizing 
brine  rejection,  the  resulting  vertical  mixing  vents  stabilizing  heat  from  deeper  water  to 
the  mixed  layer  and  ultimately  the  atmosphere.  This  eliminates  the  warming  in  the 
region  from  roughly  30-  to  65-m  depth  seen  in  Figure  5.3.6  and  acts  to  deepen  the  mixed 
layer. 

Comparison  of  the  evolution  of  the  temperature  profile  with  the  forcing  shown  in 
Figure  5.3.12  (Figure  5.3.13)  with  the  identical  case  without  air-sea  heat  fluxes  (Figure 
5.3.1)  shows  both  a  deeper  mixed  layer  (71  m  vs.  50  m)  and  maximum  depth  of  mixing 
(74  m  vs.  69  m). 
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Figure  5.3.13  Depth-time  contour  plot  of  the  temperature  response  to  a  wind  event  with  a  maximum 
speed  of  15  m/s  and  50%  maximum  open  water  fraction.  White  line  denotes  the  model-defined 
mixed  layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the 
base  of  the  slab  layer. 

The  combined  effect  of  static  and  shear  instability  can  be  seen  in  the  vertical  heat 
flux  (Figure  5.3. 14a).  The  wind  event  causes  the  initial  period  of  shear  mixing  as  in  the 
case  with  no  leads  (Figure  5.3.4),  but  it  is  combined  with  the  initial  surface 
cooling/brining  (Figure  5.3.12)  to  flux  more  heat  than  wind  alone.  The  pulsed  brining 
causes  static  instability  which  changes  the  stratification.  This  new  stratification  is 
unstable  with  the  existing  background  shear  (equation  5.6)  and  mixes  further.  In  the  case 
without  strong  wind  forcing  and  thus  no  strong  inertial  shear,  the  mixed  layer  depth  is 
identical  to  the  case  with  both  wind  and  leads,  but  the  transition  layer  has  zero  thickness 
(Figure  5.3.14b). 
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Figure  5*3,14  Depth-time  contour  plot  of  the  vertical  heat  flux  response  to  a  wind  event  with  a 
maximum  speed  of  (a)  15  m/s  and  (b)  5  m/s,  both  with  50%  maximum  open  water  fraction.  White 
line  denotes  the  model-defined  mixed  layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing 
and  the  green  line  indicates  the  hasc  of  the  slab  layer. 
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The  relative  roles  of  the  mixing  processes  contained  in  the  PWP  model  can  be 
examined  by  saving  the  temperature  and  salinity  profiles  partway  through  each  time  step. 
First,  output  is  saved  after  the  static  instability  criterion  has  been  met,  then  again  after  the 
bulk  Richardson  number  criterion  has  been  met,  and  then  after  the  gradient  Richardson 
number  criterion  has  been  satisfied.  Fluxes  can  then  be  calculated  for  each  process  in  the 
model,  as  shown  in  Figure  5.3.15  for  a  simulation  with  a  15  m/s  maximum  wind  and  50% 
maximum  open  water  fraction. 

The  PWP  model  first  checks  for  static  stability,  starting  at  the  surface,  and 
entrains  cells  into  the  mixed  layer  until  the  profile  is  stable  (Figure  5.3.15a).  The 
addition  of  salt  at  the  surface  has  caused  the  flux  of  heat  and  salt  upwards  into  the  mixed 
layer.  The  PWP  model  then  checks  the  Rib  criterion,  which  was  not  exceeded  for  the 
simulation  forced  with  a  15  m/s  maximum  wind  and  50%  maximum  open  water  fraction. 
Entrainment  into  the  mixed  layer  due  to  shear  (Rib  criterion  mixing)  did  occur  in  runs 
with  stronger  forcing,  and  played  a  significant  role  in  the  simulation  in  Section  5.4. 
Finally,  the  Rig  criterion  is  applied  and  the  model  iterates  to  make  the  flow  in  the 
transition  layer  stable  (Figure  5.3.15b). 

The  region  in  the  ‘mixed  layer’  which  shows  Rig  fluxes,  especially  between  hours 
10  and  20,  is  the  very  weakly  stratified  part  in  the  lower  portion  of  the  mixed  layer  which 
is  not  a  true  fully-mixed  slab.  The  density  difference  between  the  base  of  the  slab  and  the 
base  of  the  mixed  layer  is  of  the  order  10'5  kg/m1. 
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Figure  53.15  Depth-time  contour  plot  of  the  (a)  static  instability  and  <b)  shear  instability  vertical 
heat  flux  response  to  a  wind  event  with  a  maximum  speed  of  15  m/s  and  50%  maximum  open  water 
fraction.  White  line  denotes  the  model-defined  mixed  layer  depth,  the  red  line  marks  the  maximum 
depth  of  mixing  and  the  green  line  indicates  the  base  of  the  slab  layer. 
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The  base  ten  logarithm  of  the  inverse  of  Rig  is  a  convenient  way  of  showing  the 
potential  for  shear  instability  (Figure  5.3.16).  The  higher  the  value,  the  closer  to 
instability,  with  the  criterion  for  mixing  logio(  1/0.25)  =  0.60.  The  profile  is  rearranged 
until  the  criterion  is  met  and  no  further,  leaving  Rig  nearly  critical  (<1,  near  the  neutral 
value  of  0.25)  everywhere  in  the  transition  layer.  Therefore,  when  the  top  of  this  region 
is  changed  due  to  some  surface  forcing,  it  becomes  unstable  and  mixes  with  the  layer 
below  it.  This  nearly  critical  cell  then  becomes  unstable  and  must  mix  even  further 
down.  This  process  continues  down  to  the  maximum  depth  of  mixing  where  a  previously 
undisturbed  cell  may  need  to  be  mixed  upwards  to  satisfy  the  Rig  criterion. 

The  results  were  qualitatively  the  same  for  runs  with  different  forcing.  The  wind 
burst  causes  shear  instability  combined  with  static  instability  from  brine  rejection.  Then 
subsequent  pulsed  static  instability  deepens  the  mixed  layer  over  time.  The  transition 
layer  is  nearly  critical  to  shear  instability,  and  undergoes  frequent  rearrangement  to 
satisfy  the  stability  criterion,  particularly  whenever  another  previously  undisturbed  model 
level  is  entrained. 
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Figure  5.3.16  Depth-time  contour  plot  of  Log,0  of  the  inverse  gradient  Richardson  number  response 
to  a  wind  event  with  a  maximum  speed  of  15  m/s  and  50%  maximum  open  water  fraction.  White  line 
denotes  the  model-defined  mixed  layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing  and 
the  green  line  indicates  the  base  of  the  slab  layer. 
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The  combination  of  shear  and  static  instability  also  transports  more  passive  tracer 
up  to  the  mixed  layer  than  either  process  by  itself,  as  seen  in  the  time  series  of 
concentration  of  a  passive  tracer  in  the  model’s  surface  mixed  layer  for  runs  with  and 
without  leads  and  wind  (Figure  5.3. 1 7).  The  initial  tracer  concentration  profile 
(traccr/m3)  was  prescribed  as  identical  to  the  salinity  profile  minus  the  surface  value,  i.e. 
large  concentrations  below  the  haloclinc.  This  initial  tracer  profile  simulates  the  salty 
and  nutrient-rich  UCDW  underlying  the  nutrient-poor  surface  waters  on  the  wAP  shelf. 
There  is  no  flux  of  tracer  at  the  air-sea  or  ice-sea  interface,  which  makes  the  final  profile 
different  from  the  salinity  profile.  The  ease  without  leads  docs  not  allow  any  tracer  to  the 
surface.  The  combination  of  wind  and  leads  yields  a  concentration  of  passive  tracer  in 
the  mixed  layer  which  is  more  than  double  that  caused  by  leads  alone.  These 
experiments  demonstrate  the  critical  role  leads  play  in  fluxing  tracer  to  the  underside  of 
the  ice. 


The  flux  of  tracer  is  relevant  to  biological  processes  which  are  limited  by  nutrient 
fluxes  from  deeper  waters.  It  should  be  noted  that  those  same  processes  which  allow  for 
greater  upward  tracer  flux  also  allow  heat  to  vent  to  the  atmosphere,  and  thus  keep  the 
mixed  layer  at  or  near  the  freezing  point. 
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Figure  5.3.17  Time  series  of  passive  tracer  concentration  in  the  surface  slab  layer  for  runs  with  and 
without  leads  and  wind.  The  “leads”  and  “no  leads"  cases  have  an  \mAl  of  0.5  and  0,  respectively. 
The  “wind"  cases  have  a  maximum  wind  of  15  m/s  and  the  “light  wind"  case  has  a  maximum  wind  of 
5  m/s. 


139 


Finally,  simulations  were  repeated  with  the  open  water  fraction  constant  at  half 
the  amplitude  of  the  oscillating  open  water  fraction  of  a  comparable  run  (i.c.  Figure 
5.3. 12a)  to  examine  the  effect  of  steady  rather  than  pulsed  brining.  This  results  in  a 
nearly  identical  final  profile  of  temperature,  salinity  and  velocity.  In  other  words,  the 
mixing  averaged  over  the  simulation  period  was  the  same.  This  indicates  that  the  specific 
timing  of  the  brining  is  not  important.  However,  without  the  inertial  oscillations  and  high 
frequency  ice  motion,  there  would  be  no  mechanism  to  keep  leads  open  to  continue  the 
strong  heat  fluxes  and  brine  rejection. 

5.3.4.  Feedbacks  from  vertical  heat  flux 

Taking  the  deltas  as  temporal  changes  in  mixed  layer  properties  during  the 
simulations  with  leads  opening  and  closing,  PAS  is  three  orders  of  magnitude  larger 
thanaAF .  In  other  words,  vertical  fluxes  of  heat  have  negligible  impact  on  the  density 
of  the  mixed  layer,  and  thus  the  overall  stratification.  This  is  due  not  only  to  the  values 
of  p  and  a  in  this  region  of  T  —  S  space,  but  also  because  the  mixed  layer  temperature 
remains  at  the  freezing  point  as  excess  heat  is  vented  out  through  open  leads  and  ice  is 
formed  once  the  temperature  reaches  the  freezing  point,  producing  a  very  small  AT  in 
the  mixed  layer  during  the  simulation. 

However,  the  vertical  heat  flux  docs  impact  the  formation  of  ice  and  hence  the 
amount  of  brine  rejected  which  has  a  large  impact  on  the  final  stratification.  This  can  be 
seen  by  taking  the  total  heat  flux  specified  at  the  surface  (Figure  5.3. 12b,  black  line)  and 
calculating  the  amount  of  ice  it  produces  (taking  into  account  depression  of  the  freezing 
point  by  the  added  brine)  and  comparing  with  the  full  model.  The  15  m/s  maximum 
wind,  50%  4™,  ease  without  mixing  produces  18%  more  ice.  This  difference  is  due  to 

the  fact  that  in  the  full  model,  the  surface  heat  flux  must  first  cool  the  mixed  layer  to  the 
freezing  point  before  creating  brine.  The  excess  heat  in  the  mixed  layer  has  been  fluxed 
up  from  deeper  waters. 

The  vertical  flux  of  heat  has  little  direct  effect  on  the  stratification.  However,  it 
docs  have  a  large  effect  on  ice  and  brine  production,  which  plays  a  large  role  in  the 
stratification,  especially  in  the  presence  of  shear.  This  suppression  of  ice  and  brine 


140 


formation  causes  less  mixing,  which  in  turn  causes  less  vertical  heat  flux,  or  a  negative 
feedback.  The  actual  change  in  ice  production  will  depend  on  the  initial  profiles  of 
temperature  and  salinity,  but  is  noted  as  significant  when  compared  with  the  modeled 
20%  increase  in  ice  production  associated  with  high  frequency  ice  motion  {Hcil  and 
Hiblcr,  2002). 

5.3.5.Slowing  of  mixed  layer  due  to  entrainment 

The  high  frequency  differential  ice  motion  associated  with  the  mixing  processes 
outlined  above  is  reported  to  be  responsible  for  a  20%  increase  in  ice  volume  in  a 
simulation  of  the  Arctic  (Hcil  and  Hiblcr,  2002).  However,  Hcil  and  Hiblcr  (2002)  did 
not  take  into  account  any  oceanic  mixing  and  subsequent  effects  on  the  velocity  of  the 
mixed  laycr/ice  system.  A  change  in  mixed  layer  velocity  due  to  entrainment  could  have 
an  effect  on  a  large-scale  ice  model,  which  uses  the  velocity  of  the  mixed  laycr/ice 
system  (and  its  spatial  gradients)  in  the  calculation  of  ice  dynamics  such  as  ridging  and 
breaking.  Here  I  compare  the  final  mixed  layer  inertial  current  amplitude  for  the  case 
with  no  leads  with  the  case  with  leads.  The  decrease  in  final  mixed  layer  inertial  current 
amplitude  is  due  to  entrainment  of  deeper  fluid  which  is  at  rest. 

Figure  5.3.1 8  shows  that  the  mixed  layer  slows  by  up  to  10%  compared  to  the 
case  with  no  open  water.  Furthermore,  one  can  sec  that  the  combined  effect  of  shear  and 
static  instability  causes  the  largest  relative  slowing  of  the  mixed  layer. 
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Figure  5,3.18  Slowing  of  mixed  layer  due  to  lead  convection  and  shear  processes.  The  slowing  is 
shown  as  a  percentage  decrease  compared  to  the  case  of  no  leads  for  each  wind  speed. 

The  specific  results  shown  in  Figure  5.3.18  arc  sensitive  to  the  initial  and 
boundary  conditions.  The  main  point  is  to  show  that  if  one  chooses  to  include  the 
mechanism  of  high-frequency  ice  motion  in  creating  a  significant  icc  volume,  one  must 
also  consider  entrainment  which  slows  the  mixed  layer,  and  hence  the  ice. 

5.3.6.  Effects  of  a  Rotating  Wind 

The  oceanic  mixed  layer  currents  respond  primarily  to  the  component  of  wind 
stress  rotating  at  the  inertial  frequency  {counterclockwise  in  the  Southern  Hemisphere) 
(c.g.  D’Asaro,  1985;  Pollard  and  Millard,  1970).  Ail  simulations  presented  thus  far  have 
been  for  a  purely  eastward  wind  event.  In  this  section  the  effects  of  a  rotating  wind  in  the 
presence  of  icc  are  examined. 

In  Figure  5.3.19  the  red  dots  show  the  rotation  versus  wind  stress  recorded  at  the 
Dismal  Island  Automatic  Weather  Station  during  2001  and  2002.  The  blue  dots  indicate 
October,  November  and  December,  2002,  the  period  with  lower  observed  icc 
concentration  and  more  inertial  energy.  The  strong  winds  (greater  than  0.2  N/nr)  were 
observed  to  rotate  very  little,  mostly  between  plus  and  minus  10  deg/hr,  far  from  the 
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inertial  frequency  (28  deg/hr).  The  black  dotted  line  in  Figure  5.3. 19  shows  the  winds  for 
October  26, 2002,  the  strong  forcing  event  which  produced  the  observed  near-inertial 
oscillations  described  in  the  previous  chapter.  For  reference,  the  wind  stress  from  a  10 
m/s  wind  is  0.14  N/m2  and  a  20  m/s  wind  yields  a  stress  of  0.87  N/m:,  assuming  neutral 
stability  (Large  and  Pond,  1981). 

Despite  the  fact  that  strong  winds  were  not  observed  to  have  strong  rotation  in  the 
wAP  region,  the  basic  PWP  with  ice  model  was  run  with  a  rotating  wind  event  at 
different  frequencies.  A  wind  burst  of  10  m/s  and  no  rotation  generates  small  inertia! 
motions  of  3  cm/s  (cyan  curve  in  Figure  5.3.19).  A  wind  of  10  m/s  rotating  at  the  inertial 
frequency  docs  not  generate  inertial  oscillations  as  strong  as  a  non-rotating  wind  of  20 
m/s.  A  20  m/s  wind  rotating  at  the  inertial  frequency  drives  just  over  double  the  inertial 
currents  as  the  non-rotating  case.  Conversely,  a  wind  rotating  at  the  inertial  frequency  in 
the  opposite  sense  to /generates  virtually  no  inertial  current.  Note  that  true  resonance 
(and  unrealistically  high  mixed  layer  velocities)  at/can  only  occur  in  an  undamped  slab 
model.  Any  model  with  damping,  in  this  ease  in  the  form  of  vertical  mixing  of 
momentum,  has  a  restriction  on  mixed  layer  velocity. 
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Figure  53 A 9  Observed  rotation  versus  wind  stress  {red  dots),  with  blue  denotes  denoting  the  period 
of  lower  ice  concentration  in  October  through  December,  2002.  The  black  dotted  line  shows  October 
26,  2002.  The  vertical  dashed  black  line  shows  the  inertial  frequency.  The  green  curve  shows  the 
modeled  final  ocean  mixed  layer  surface  velocity  (leftmost  y-axis)  when  forced  by  a  2ft  m/s  wind 
rotating  at  the  frequency  denoted  on  the  %-axis.  The  cyan  curve  shows  the  results  for  a  rotating  10 
m/s  wind. 


When  the  observed  wind  stress  vector  in  the  wAP  region  was  strong  enough  to 
generate  notable  near-inertial  motions,  it  did  not  rotate  enough  to  make  a  great  difference 
to  the  overall  conclusions  in  this  study.  The  inertial  period  is  shorter  than  the  timescale 
for  rotation  of  strong  winds  in  this  region.  This,  however,  is  not  always  the  case  at  lower 
latitudes  with  longer  inertial  periods  or  other  high-latitude  regions  with  faster  rotating 
winds. 


5.4.  Application  of  idealized  study  to  wAP  region 

Assessing  the  role  of  near-inertial  oscillations/high  frequency  ice  motion  and 
production  in  mixing  the  waters  of  the  wAP  region  presents  significant  challenges.  First, 
the  air-sea  heat  fluxes  arc  not  well  constrained,  given  that  the  net  radiative  flux  plays  a 
significant  role  (Beardsley  ct  al.,  2006)  and  was  only  measured  by  shipboard 
instrumentation.  Ships  were  not  present  in  the  region  for  much  of  the  study  period,  in 
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particular  October  through  December  when  high  near-inertial  activity  was  observed. 
Therefore,  the  surface  boundary  condition  is  not  well  known.  Second,  the  lack  of  initial 
conditions  from  CTD  casts,  not  only  to  initialize  the  computations,  but  also  to  compare 
with  the  results,  further  precludes  any  detailed  simulations.  Furthermore,  the  mixed  layer 
depth  is  unknown  for  most  of  the  year  since  there  were  no  temperature  or  salinity 
measurements  above  104  m.  Finally,  the  wind  stress  has  a  significant  portion  of  its 
energy  put  into  the  internal  ice  stress  term  during  the  strong  wind  event  as  described  in 
Chapter  4,  the  modeling  of  which  is  beyond  the  scope  of  this  study. 

Despite  these  obstacles,  in  this  section  I  make  a  rough  estimate  of  the  potential 
mixing  caused  by  shear  and  static  instabilities.  The  results  of  the  idealized  experiments 
reveal  several  key  features  which  allow  for  simplification  of  the  modeling.  First,  while 
static  and  shear  instabilities  together  create  more  mixing  than  either  alone,  the  timing  of 
the  buoyancy  and  wind  stress  forcing  is  not  critical  at  higher  (inertial)  frequencies. 
Therefore  one  can  apply  a  heat  flux  on  the  daily-averaged  open  water  fraction  and  make 
the  assumption  that  high  frequency  motions  arc  keeping  the  leads  open.  Next,  in  order  to 
work  around  the  complicated  and  often  strong  observed  internal  ice  stresses,  i  apply  the 
observed  ice  motion  as  the  surface  boundary  condition  and  calculate  the  ice-ocean  stress 
directly  from  that  (observed  winds  arc  not  used).  I  linearly  interpolate  the  available 
moored  temperature  and  salinity  measurements,  and  assume  an  initial  mixed  layer  of  65 
m  with  the  same  temperature  and  salinity  properties  of  the  CTD  cast  shown  in  Figure 
5.2.2.  Again,  the  net  radiative  air-sea  heat  flux  is  assumed  to  be  a  constant  1 00  W/m2 
loss. 

I  concentrate  on  the  inertial  event  described  in  Chapter  4  at  the  C3  site,  October 
25  through  October  31,  2002.  The  end  of  the  run  is  chosen  just  before  the  “non-local” 
forcing  due  to  divergence  of  internal  ice  stresses  occurs.  This  October  26  storm  is  a 
strong  event,  with  slightly  weaker  but  noteworthy  inertial  current  events  occurring 
roughly  on  the  order  of  once  a  month  (Figure  5.1  .la).  Figure  5.4. 1  shows  the  boundary 
conditions  used  to  force  the  run.  There  is  a  large  net  heat  flux  during  the  wind  event 
followed  by  a  period  of  relatively  warm  air  (but  still  at  or  below  the  sea  surface 
temperature)  and  heat  flux  comprised  primarily  of  the  specified  100  W/m2  radiative  loss. 
Warm  air  is  typically  associated  with  the  strong  “fohn”  winds  in  Marguerite  Bay 
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(Beardsley  ct  ah*  2006).  There  is  an  initial  southward  velocity  in  the  ice,  followed  by 
near-inertial  oscillations*  The  time  step  in  this  simulation  is  15  minutes.  The  simulation 
start  time  is  October  25,  2002: 
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Figure  5.4,1  Time  series  of  (a)  total  surface  heat  flux  before  multiplying  by  the  open  water  fraction, 
ice  concentration  (l-A),  and  the  resulting  net  heat  flux,  (b)  Ice  velocity. 

The  simulation  is  not  meant  to  reproduce  the  observations,  but  rather  to  give  some 
insight  into  a  plausible  mixing  scenario  in  Marguerite  Bay.  The  modeled  temperature, 
salinity  and  velocity  fields  arc  qualitatively  similar  to  those  for  the  idealized  experiments, 
with  the  initial  wind  event  causing  large  changes,  with  smaller  changes  during  the  inertial 
oscillations  (Figure  5.4.2).  The  initial  wind  event  causes  current  velocities  of  about  1 2 
cm/s,  with  the  subsequent  oscillations  of  6  to  8  cm/s,  similar  to  what  was  observed. 

About  4.5  cm  of  ice  was  formed  in  the  model  during  the  first  one  and  a  half  days,  writh  no 
subsequent  ice  formation  and  brine  rejection. 
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Figure  5.4.2  Depth-time  contour  plot  of  the  (a)  temperature,  (b)  salinity  and  (c)  speed  response  to  a 
simulation  based  on  an  observed  event  in  October,  2002.  White  line  denotes  the  model-defined  mixed 
layer  depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of 
the  slab  layer. 


The  model  vertical  heat  flux  is  shown  in  Figure  5.4.3.  The  maximum  depth  of 
mixing  descends  quickly  during  the  large  surface  forcing  (both  stress  and  density),  then 
remains  somewhat  steady,  consistent  with  the  findings  earlier  in  this  chapter.  After  the 
initial  wind  event,  the  mixed  layer  depth  moves  up  and  down,  with  the  main  pycnoclinc 
remaining  at  the  maximum  depth  of  mixing.  The  movement  of  the  mixed  layer  depth  is 
due  to  a  weak  stratification  in  conjunction  with  weak  shear  smoothing  the  gradient 


148 


Total  Flux 


Figure  5.4,3  Depth-time  contour  plot  of  the  vertical  heat  flux  response  to  a  simulation  based  on  an 
observed  event  in  October,  2002,  White  line  denotes  the  model-defined  mixed  layer  depth,  the  red 
line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of  the  slab  layer. 
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Figure  5.4.4  Depth-time  contour  plots  of  the  contributions  to  the  vertical  heal  flux  response  from  (a) 
static  (b)  bulk  Richardson  number  and  (c)  gradient  Richardson  number  instability  in  a  simulation 
based  on  an  observed  event  in  October,  2002.  While  line  denotes  the  model-defined  mixed  layer 
depth,  the  red  line  marks  the  maximum  depth  of  mixing  and  the  green  line  indicates  the  base  of  the 
slab  layer. 

The  flux  resulting  from  the  individual  mixing  mechanisms  can  be  separated  for 
analysis  (Figure  5.4.4).  Before  the  wind  burst  on  October  26,  there  is  no  shear  mixing, 
and  static  instability  due  to  brine  rejection  causes  a  deepening  of  the  mixed  layer  from  65 
m  to  70  m.  Then  (he  icc  formation  and  brine  rejection  continue  for  an  additional  halfday 
(first  half  of  October  26),  the  first  half  of  the  initial  wind  burst  and  sea  icc  movement. 
During  this  -1/2  day  period  of  strong  forcing,  the  addition  of  salt  and  momentum  cause 
entrainment  into  the  mixed  layer  by  bulk  Richardson  number  mixing.  The  associated 
changes  in  the  profiles  of  density  and  velocity  at  the  top  of  the  transition  layer  triggers 
gradient  Richardson  number  mixing,  which  deepens  the  maximum  depth  of  mixing  to  90 
m.  The  heat  fluxed  upward  from  the  deeper  waters  is  released  to  the  atmosphere  until 
midday  October  26,  when  the  surface  heat  flux  is  reduced  (Figure  5.4. 1 )  due  to  higher  air 
temperature  and  humidity,  yet  the  strong  surface  stress  continues  for  another  half  day. 
This  results  in  a  warming  of  the  surface  layer  to  about  -0. 1  °C  above  the  freezing  point. 
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After  the  strong  forcing  cvent.there  is  both  weak  surface  cooling  and  weak  shear. 
The  vertical  bands  of  flux  arc  due  to  the  surface  being  relatively  warm  after  the  strong 
wind  event  and  then  subsequent  weak  surface  cooling  adding  enough  density  without  ice 
formation  to  convcct  partway  down  before  the  upward  flux  of  salt  stabilizes  the  density 
profile.  The  static  instability  docs  not  penetrate  into  the  deep  reservoir  of  heat,  while  the 
weak  shear  creates  some  mixing  to  smooth  the  gradients.  There  is  little  change  in  the 
temperature  and  salinity  profiles  after  the  strong  forcing  (Figure  5.4.2a,b). 

Overall,  the  simulation  reflects  the  salient  features  shown  in  the  idealized 
experiments.  The  initial  wind  burst  performs  the  majority  of  the  mixing  with  subsequent 
oscillations  stable  to  shear  instability.  The  surface  heat  flux  and  brining  keeps  the  surface 
mixed  layer  cold  and  changes  the  stratification,  which  causes  further  shear  instability. 
However,  the  strong  surface  stress  persists  longer  than  the  cooling,  causing  warming  of 
the  surface  waters  above  freezing.  Thus,  the  duration  of  the  strong  surface  heat  flux  and 
stress  arc  important.  After  the  strong  event,  mixing  (deepening  maximum  depth  of 
mixing)  docs  not  penetrate  significantly  into  the  deep  reservoir  of  heat  and  salt. 

Therefore  it  is  the  initial  wind  event  which  is  of  greatest  interest  to  studies  of  mixing  in 
Marguerite  Bay,  especially  since  the  wind  events  lend  to  be  followed  by  relatively  warm 
(-5  to  0  °C)  air  (Beardsley  ct  al.,  2006). 

The  upward  heat  flux  evaluated  at  the  base  of  the  mixed  layer  is  quite  large  for  a 
short  period,  with  a  peak  value  of  ~  900  W/m:  and  an  average  of  385  W/m:  for  October 
26.  The  average  heat  loss  for  October  27-30  is  1 7  W/m:,  w  hich  can  be  compared  to  the 
1-2  W/m’  estimated  from  an  analysis  of  observed  near-inertial  oscillations  on  the  wAP 
shelf  in  Chapter  2.  Furthermore,  low  values  of  heat  flux  (0-5  W/m2)  have  been  found 
with  a  CTD-mountcd  microstructurc  measurements,  but  these  data  arc  still  being 
processed  and  analyzed  so  it  is  not  clear  if  they  were  taken  during  strong  atmospheric 
forcing  events  (L.  Padman,  pers.  comm.). 

The  time-averaged  vertical  heat  flux  for  the  entire  6-day  simulation  is  79  W/ntf. 
Making  the  rough  assumption  of  one  large  event  every  month  during  the  ice  season,  (or 
equivalently  one  day  of  385  W/m"  and  29  days  of  17  W/m  )  gives  an  average  flux  of -29 
W/m:  across  the  base  of  the  mixed  layer.  This  can  be  compared  with  estimates  of  an 
average  vertical  heat  flux  across  the  pycnoclinc  based  on  wAP  hydrography  of  between  5 
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and  10  W/m"  (Klinck,  1998;  Smith  ct  al.,  1999).  Note  that  while  an  annual  average  is 
calculated,  the  conditions  under  which  the  6-day  simulation  was  made  do  not  occur 
throughout  the  year,  lee  covers  the  wAP  shelf  about  seven  months  a  year  (Chapter  4), 
and  through  half  of  that  interval  strong  observed  near-inertial  motions  arc  not  observed 
(Figure  5.1.1).  Overall,  this  work  supports  the  idea  of  a  weak  average  vertical  heat  flux 
across  the  pycnocline  in  Marguerite  Bay  made  up  of  episodic  strong  events  which  may  be 
difficult  to  capture  obscrvationally.  In  this  rough  calculation,  nearly  half  of  the  total  flux 
across  the  base  of  the  pycnocline  during  a  month  occurs  during  a  one-day  event. 
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Figure  5.4.5  Time  series  of  downward  vertical  heat  flux  at  the  model  mixed  layer  depth  (black  line) 
and  the  individual  contributions  from  static  (blue),  bulk  Richardson  number  (red)  and  gradient 
Richardson  number  (green)  instability  in  a  simulation  based  on  an  observed  event  in  October,  2002. 


5.5.  Summary  and  Discussion 

In  this  chapter  I  have  used  the  PWP  one-dimensional  upper  ocean  mode!  to 
examine  vertical  mixing  beneath  sea  ice,  with  an  emphasis  on  inertial  oscillations  and 
high  frequency  ice  motion.  In  the  parameter  space  tested,  nearly  all  mixing  occurs  within 
the  first  inertial  period,  consistent  with  Plucddcmann  and  Farrar  (2006).  The  presence  of 
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icc  acts  to  smooth  the  surface  stress  on  the  ocean  in  time,  but  docs  not  reduce  mixing  for 
reasonable  values  of  the  ice-ocean  drag  coefficient  (0.015).  The  response  of  the  mode! 
begins  to  diverge  from  the  open  water  ease  with  drag  coefficients  less  than  0.001 . 

Overall  the  presence  of  ice  has  little  effect  on  the  flows  and  mixing. 

The  model  results  revealed  an  interesting  synergy  between  wind  and  buoyancy 
forcing  in  ocean  mixing  beneath  ice  in  inertial  events.  High  frequency  icc  motion, 
particularly  leads  opening  and  closing,  had  a  great  effect.  The  active  surface  heat  flux 
allows  for  heat  loss  and  subsequent  icc  formation.  The  resulting  brine  rejection  leads  to 
periodic  changes  in  the  stratification,  causing  subsequent  shear  instability  and  more 
mixing  than  that  caused  by  air-sea  heat  flux  atone.  The  leads  also  allow  atmospheric 
venting  of  heat  mixed  up  from  the  deeper  water. 

The  oscillatory  opening  and  closing  of  leads  introduces  an  average  open  water 
fraction  that  allows  for  heat  loss  to  the  atmosphere.  It  is  not  the  details  of  the  oscillations, 
but  the  open  water  they  create,  which  has  a  great  effect.  Furthermore,  entrainment  of 
slower  fluid  into  the  mixed  layer  decreases  the  mixed  layer  speed  by  up  to  1 0%.  This 
could  affect  icc  models  such  as  those  employed  in  Hcil  and  Hiblcr  (2002). 

Based  on  a  simplified  simulation  of  an  observed  event,  the  modeled  vertical  heat 
flux  across  the  base  of  the  mixed  layer  during  a  strong  forcing  event  is  quite  large,  -385 
W/m2  averaged  over  one  day.  The  remainder  of  the  simulation  was  characterized  by 
weak  fluxes,  with  the  6-day  average  value  over  the  entire  simulation  coming  to  79  W/m2. 
This  translates  into  an  estimated  overall  average  heat  flux  of  29  W/m',  assuming  roughly 
one  strong  event  per  month  with  the  remainder  of  the  month  having  weak  fluxes.  The 
heat  gained  by  the  mixed  layer  from  the  UCDW  was  completely  fluxed  to  the  atmosphere 
by  the  end  of  the  6-day  simulation  (see  Appendix  5  A. 2) 

A  deep  flow  180°  out  of  phase  with  the  mixed  layer  flow  was  observed  in  the 
shipboard  ADCP  data  (Chapter  2),  but  not  simulated  in  the  one-dimensional  PWP  model. 
This  deep  flow  is  probably  due  to  a  barotropic  wave  reflecting  from  the  coastal  boundary. 
This  deep  flow  docs  not  appear  to  have  a  great  effect  on  mixing  (see  Appendix  5 A.  1 ). 

The  structure  in  the  moored  time  series  temperature  data  (last  Figure  from 
Chapter  4)  suggest  something  more  is  happening  in  Marguerite  Bay.  One  possibility  is 
dense  brine  plumes  penetrating  the  pycnoclinc  and  generating  internal  waves.  This  could 
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potentially  be  investigated  in  a  three  dimensional  high  resolution  simulation,  but  is 
beyond  the  scope  of  this  study.  Another  possibility  is  that  large  wind  events  set  deep- 
keeled  icebergs  in  motion  that  drives  mixing  and  radiate  internal  waves. 

It  should  be  noted  that  the  divergent  high  frequency  ice  motion  can  occur  in  the 
absence  of  a  coast.  Other  causes  for  such  divergent  motions  could  be  a  horizontally 
inhomogeneous  wind  forcing  field  or  a  horizontally  inhomogeneous  ice  strength  Held 
(Hcil  and  Hibler,  2002),  In  fact,  Hcil  and  Hiblcr  (2002)  found  that  even  with  a 
temporally-smoothed  wind  forcing  field,  inertial  oscillations  in  the  ice  occurred  due  to 
the  nonlinearity  of  the  sea  ice  model.  The  ice  model  can  produce  abrupt  ice  motions 
when  strength  thresholds  arc  exceeded,  and  thus  creating  temporally  abrupt  and  spatially 
inhomogeneous  ice-ocean  stresses. 


APPENDIX  5.A.1  Coastal  Boundary  Effects 

The  modeled  response  of  an  ice-covered  ocean  thus  far  has  not  included  the  flow 
observed  in  the  lower  layer  which  is  1 80°  out  of  phase  with  the  upper  layer  flow  (Chapter 
2).  This  lower  layer  flow  is  most  likely  the  result  of  a  barotropic  wave  reflecting  off  of 
the  coastal  boundary  (Pettigrew,  1981 ;  Kundu  ct  al.,  1983).  This  lower  flow  is  predicted 

to  begin  at  time  t  -  xi  -J~gH  ,  where  x  is  the  distance  to  the  coast,  H  is  the  total  water 

depth,  and  t  =  0  is  the  initiation  of  the  surface  wind  forcing.  This  is  about  1 8  minutes  for 
the  mooring  sites  in  Marguerite  Bay,  and  therefore  may  be  considered  instantaneous. 
Furthermore,  the  observation  that  the  inertial  band  ice  velocities  are  greater  than  the 
water  velocities  (Chapter  4),  may  also  be  due  to  these  coastal  boundary  effects.  Since 
this  flow  is  part  of  a  fundamentally  three-dimensional  process,  the  one-dimensional  PWP 
model  can  reproduce  neither  the  lower  layer  flow  nor  the  inertial  band  velocities  which 
arc  stronger  in  ice  than  the  water. 

Despite  this  obstacle,  an  effort  was  made  to  investigate  the  potential  role  of  the 
opposite  phased  lower  layer  flow  by  using  an  ad  hoc  approach  of  simulating  oscillations 
after  the  wind  event  has  occurred.  This  was  accomplished  by  taking  the  mixed  layer 
velocity  resulting  from  a  certain  strength  wind  event  in  a  run  with  no  leads  (  /tnwJ(  —  0 ), 
then  using  the  resulting  surface  current  speed  to  scale  the  mode  I  structure  of  the 
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strati  Heat  ion  profile.  The  modal  structure  was  made  using  the  full  observed  stratification 
(not  the  idealized  linear  profile)  shown  in  Figure  5.2.2b.  This  sealed  mode  1  velocity 
profile  was  then  used  as  the  initial  condition  in  a  run  with  no  wind  stress  forcing,  and  a 
specified  varying  open  water  fraction  as  described  earlier.  The  mode  1  velocity  profile 
was  chosen  because  it  resembles  the  observed  velocity  profile  presented  in  Chapter  2, 
Figure  4. 

The  resulting  heat  flux  is  nearly  identical  to  a  simulation  with  leads  and  no  wind 
forcing  at  all  (Figure  5.3. 14).  The  modal  velocity  structure  docs  not  provide  enough 
shear  to  alter  the  mixing  when  compared  with  a  simulation  which  begins  at  rest.  This 
result  is  robust  for  a  range  of  velocities  and  open  water  fractions. 

The  modal  structure  has  a  gradient  Richardson  number  profile  which  is  finite 
throughout  the  water  column  below  the  mixed  layer  (where  Rig  is  undefined),  i.c.  where 
there  is  non-zero  shear  and  stratification  (blue  line  in  Figure  5. 5 A.  1 ).  The  model  runs 
which  begin  at  rest  spin  up  inertial  oscillations  which  have  a  Rig  which  is  undefined  in 
both  the  mixed  layer  and  infinite  below  the  maximum  depth  of  mixing,  where  there  is  no 
shear  (red  line  in  Figure  5.5A.1).  However,  in  the  transition  layer,  Rig  is  critical,  so  any 
destabilizing  change  to  the  profile  results  in  mixing. 

The  barely  stable  profiles  of  temperature,  salinity  and  velocity  which  result  from 
the  initial  wind  event  arc  vital  to  proper  simulation  of  response  to  subsequent  forcing. 
Therefore,  the  coastal  boundary  effects  must  be  simulated,  including  the  actual  wind 
event,  with  a  true  higher-dimensional  model  to  investigate  the  details  of  vertical  mixing 
beneath  high  frequency  sea  ice  motion. 
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Figure  5.5A.1  Profile  of  Ioglu(I/  Rig)  at  the  end  of  two  simulations,  the  red  having  no  initial  modal 
velocity  structure,  the  blue  with  an  initial  velocity  structure. 


APPENDIX  5.A.2  Mixed  Layer  Heat  Balance 


In  order  to  investigate  the  fate  of  heat  fluxed  from  the  UCDW  into  the  AASW  in 
the  simulation  described  in  section  5.4,  one  can  also  calculate  a  depth-integrated  heal 
balance  for  the  mixed  layer  (c.g.  McPhadcn,  1982),  ignoring  lateral  effects: 


where  QHF{ 0)  and  Q,IF(~h)  arc  positive  upwards  and  T  is  the  average  temperature  in  the 
mixed  layer  of  thickness  h ,  Q„F( 0)  represents  only  the  heat  flux  at  the  surface  after  any 
latent  heat  of  fusion  in  ice  production  has  been  removed.  On  the  second  day,  when  the 
surface  stress  is  strongest,  upward  heat  flux  through  the  base  of  the  mixed  layer  is 
primarily  balanced  by  a  wanning,  with  some  loss  of  heat  to  the  atmosphere  (Figure 
5.5A.2).  For  the  remainder  of  the  simulation,  the  balance  is  maintained  by  a  reduction  in 
mixed  layer  temperature  as  heat  is  lost  to  the  atmosphere. 
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Misfed  layer  Heat  Balance 


Figure  5.5A.2  Terms  in  the  mixed  layer  heat  balance  in  a  simulation  based  on  an  observed  event  in 
October,  2002.  The  sign  convention  has  the  terms  sum  to  zero  for  clarity. 

Integrating  these  terms  in  time  (Figure  5.5A.3)  shows  that  the  heat  content  in  the 
mixed  layer  returns  to  its  original  value  at  the  end  of  the  6-day  simulation.  The  heat 
gained  by  the  mixed  layer  from  the  deeper  waters  during  the  wind  event  is  gradually  lost 
to  the  atmosphere  during  the  following  four  days  of  cooling.  This  suggests  the 
importance  of  the  cooling  between  strong  events  to  maintain  the  wintertime  temperature 
profile.  Heat  from  the  UC'DW  effectively  moves  through  the  mixed  layer  to  the 
atmosphere.  The  heat  lost  from  the  UCDW  must  be  replenished  through  some  process  not 
included  in  the  P WP  model,  such  as  lateral  advection  of  UCDW  from  offshore  onto  the 
shelf. 


158 


x  -jo4  Integrated  Mixed  Layer  Heat  Balance 


Figure  5.5A.3  Time-integrated  terms  in  the  mixed  layer  heat  balance  in  a  simulation  based  on  an 
observed  event  in  October,  2002.  The  sign  convention  has  the  terms  sum  to  zero  for  clarity. 
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6.  Chapter  6 


Summary  and  Discussion 

6.1.  Summary 

Subsurface  intrusions  of  Upper  Circumpolar  Deep  Water  (UCDW)  onto  the 
western  Antarctic  Peninsula  (wAP)  shelf  are  believed  to  provide  much  of  the  heat,  salt 
and  nutrient  flux  into  the  region  (Smith  et  al,  1999).  In  order  to  investigate  the  possible 
mechanisms  responsible  for  this  flux,  moored  oceanographic  and  automatic  weather 
station  data  were  collected  within  Marguerite  Bay  on  the  shelf  as  part  of  the  Southern 
Ocean  Global  Ocean  Ecosystems  Dynamics  (SO  GLOBEC)  program.  The  subsurface 
moorings  were  deployed  in  the  ice-free  austral  summer  and  observed  the  onset,  draft  and 
motion  of  sea  ice  during  the  winter  ice  seasons  of  2001-2003,  in  addition  to  oceanic 
temperature,  salinity,  pressure  and  velocity  at  various  depths.  In  this  thesis  I  develop  a 
novel  method  for  observing  sea  ice  with  an  upward-looking  acoustic  Doppler  current 
profiler  (ADCP)  and  use  it  to  characterize  the  sea  ice  within  Marguerite  Bay  on  the  wAP 
shelf;  then  I  investigate  the  possible  processes  responsible  for  vertical  fluxes  of  heat,  salt 
and  nutrients  in  an  ice-covered  ocean. 

In  Chapter  2,  initial  estimates  were  made  of  upward  diapycnal  fluxes  of  heat,  salt 
and  nutrients  from  the  UCDW  to  the  surface  mixed  layer  due  to  two  small-scale 
processes,  double  diffusive  convection  (DDC)  and  shear  instability.  The  estimates  were 
based  on  data  from  conductivity-temperature-depth  (CTD)  and  vessel-mounted  ADCP 
data  collected  during  three  research  cruises  between  April  and  September  2001 .  Near- 
inertial  motions  generated  by  wind  stress  provided  most  of  the  shear  across  the  sharp 
pycnocline  at  the  base  of  the  mixed  layer.  The  mean  vertical  diffusivity  associated  with 
shear  instability  was  estimated  at  <1  x  10  '  mV,  corresponding  to  a  heat  flux  into  the 
base  of  the  mixed  layer  of  <2  W  m'2.  While  prior  work  suggests  that  DDC  could  provide 
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upward  heat  fluxes  of  order  10  W  m 2  on  the  wAP  shelf  (Marguerite  Bay  north  to  the  tip 
of  the  Antarctic  Peninsula)  (Smith  and  Klinck,  2002),  the  analyses  in  Chapter  2  indicate 
almost  no  contribution  to  diapycnal  fluxes  from  DDC  (<0.2-0.4  W  m' ).  A  cursory 
analysis  of  one  observed  wind  event  prompted  the  suggestion  that  mixing  on  the  wAP 
shelf  could  be  dominated  by  episodic  events  not  captured  in  the  hydrographic  survey 
data. 

In  Chapter  3,  a  methodology  was  presented  for  determining  ice  coverage,  draft 
and  velocity  from  data  from  a  moored  upward-looking  ADCP.  The  technique  was 
successfully  applied  to  data  collected  on  the  wAP  shelf  during  SO  GLOBEC.  Ice 
detection  was  accomplished  using  windowed  variances  of  ADCP  vertical  velocity  and 
error  velocity,  with  the  daily-average  ice  concentration  estimated  as  the  percentage  of 
time  ice  in  which  was  detected  during  a  given  day.  The  ADCP-derived  ice  concentration 
compared  well  with  a  satellite-derived  estimate  (Cavalieri  et  al„  1 997).  Furthermore,  the 
ADCP-derived  ice  concentration  data  improved  on  observation  of  ice  concentrations  on 
small  spatial  scales,  observing  a  localized  polynya  over  the  Cl  mooring  site  which 
appeared  as  a  broad  area  of  low  ice  concentration  in  the  satellite-derived  ice 
concentration  data  (Chapter  4).  Acoustic  signal  correlation  and  backscatter  intensity 
were  found  to  be  poor  indicators  of  the  presence  of  ice  at  this  site.  Ice  draft  was 
estimated  using  a  combination  of  ADCP  backscatter,  atmospheric  and  oceanic  pressure 
data,  and  information  about  the  stratification.  This  estimate  required  corrections  to  the 
ADCP-derived  range  for  instrument  tilt  and  sound  speed  profile.  Uncertainties  of  ±  0.20 
m  during  midwinter  and  ±  0.40  m  when  the  base  of  the  mixed  layer  was  above  the  ADCP 
for  ice  draft  were  estimated  based  on  (a)  a  Monte  Carlo  simulation,  (b)  uncertainty  in  the 
sound  speed  correction  and  (c)  performance  of  the  zero-draft  estimate  during  times  of 
known  open  water.  Ice  velocity  was  taken  as  the  ADCP-measurcd  horizontal  velocity  in 
the  depth  bin  specified  by  the  range  estimate.  This  method  is  currently  in  use  by  other 
researchers  in  both  the  Arctic  (F.  Straneo,  C.  Lee)  and  Antarctic  (M.  Meredith). 

Presented  in  Chapter  4  are  time  series,  spanning  the  winter  ice  seasons  of  2001- 
2003,  of  sea  ice  concentration,  thickness  and  motion  within  Marguerite  Bay  that  were 
estimated  using  the  methodology  described  in  Chapter  3.  Both  years  have  roughly  seven 
months  of  nearly  complete  ice  cover,  but  the  onset  of  ice  formation  was  about  two 
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months  earlier  in  2002-2003  than  in  the  prior  year  due  to  one  of  the  coldest  May/June 
periods  on  record  since  1976.  Satellite  data  show  that  the  onset  of  ice  at  the  mouth  of 
Marguerite  Bay  correlates  well  with  the  first  instance  of  -9.2°  C  air  temperatures  during 
the  period  of  1978-2002.  Sizeable  ice  draft  of  3  m  or  more  lasted  for  weeks  during  both 
observed  winters,  with  an  annual  mean  sub-inertial  ice  velocity  of  ~10  cm  s'1  during  the 
ice-covered  period.  A  linear  momentum  balance  shows  the  importance  of  internal  ice 
stresses  that  are  necessary  to  explain  the  observed  motion  of  the  ice  pack.  These  internal 
ice  stresses  are  consistent  with  the  observed  high  concentrations  and  deep  drafts  of  the  ice 
and  the  semi-enclosed  nature  of  Marguerite  Bay.  Tidal  motions  are  weak,  and  most  of 
the  kinetic  energy  in  the  ice  and  upper  ocean  is  in  the  near-inertial  band. 

In  Chapter  5  the  PWP  (Price  et  al,  1986)  one-dimensional  vertical  mixed  layer 
model  was  adapted  for  investigation  of  mixing  beneath  an  ice-covered  ocean. 

Simulations  using  an  idealized  storm  forcing  of  a  stratification  representative  of  winter  on 
the  wAP  shelf  revealed  the  relative  roles  of  wind  mixing  and  convection  due  to  brine 
rejection  in  vertical  mixing.  When  reasonable  ice-ocean  drag  coefficients  were  used 
(0.015),  the  presence  of  ice  did  not  greatly  change  the  modeled  response  of  the  water 
column  when  compared  with  the  open  water  case.  Shear  mixing  driven  by  the  initial 
wind  event,  not  the  subsequent  shear  associated  with  inertial  oscillations,  caused  the 
mixing  and  subsequent  deepening  of  the  mixed  layer. 

The  model  results  indicate  that  the  combined  effects  of  shear  and  static  instability 
can  be  significant  over  short  time  scales.  The  model  forced  by  observed  atmospheric  and 
ice  conditions  on  the  wAP  shelf  produces  large  vertical  fluxes  for  short  periods,  a  -385 
W/nr  one-day  averaged  heat  tlux  across  the  base  of  the  mixed  layer  during  a  storm  event. 
However,  when  averaged  over  the  entire  6-day  simulation,  the  tlux  reduces  to  a  more 
modest  79  W/m:  across  the  base  of  the  mixed  layer.  Making  the  rough  assumption  of 
one  storm  event  per  month  with  the  remainder  of  the  month  having  a  weak  flux  ( 1 7 
W/nf )  one  estimates  a  monthly  average  tlux  of  29  W/nrr  across  the  base  of  the  mixed 
layer. 

The  small  (17  W/ni:)  fluxes  estimated  during  the  period  of  inertial  oscillations 
following  the  storm  event  are  consistent  with  the  estimate  made  in  Chapter  2  of  ~5  W  m  * 
over  the  sharp  pycnocline  during  a  4-day  period  of  inertial  oscillations  in  August,  2001. 
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Furthermore,  the  results  in  Chapter  5  support  the  notion  proposed  in  Chapter  2  that 
intermittent  but  relatively  energetic  events  may  be  the  dominant  contributors  to  mean 
turbulent  diapycnal  fluxes. 

6.2.  Discussion 

Prior  to  the  SO  GLOBEC  observational  program  the  average  upward  heat  flux 
through  the  pycnocline  was  thought  to  be  about  10  W  m  2  (Klinck,  1998;  Smith  et  al„ 
1999),  with  a  significant  contribution  from  DDC  based  on  the  need  for  different 
diffusivities  of  heat  and  salt  to  close  the  annual  budgets  (Smith  and  Klinck,  2002).  These 
initial  studies  calculated  a  net  heat  loss  from  the  upper  ocean  based  on  budgets  in  a  one¬ 
dimensional  model  and  observed  hydrography.  These  estimates  rely  on  the  net 
differences  of  large  air-sea  heat  fluxes  over  the  annual  cycle,  themselves  containing  non¬ 
trivial  uncertainty.  In  Chapter  2,  an  analysis  of  observations  made  in  the  first  year  of  SO 
GLOBEC  could  not  find  evidence  of  notable  DDC  fluxes  {<0.2  -  0.4  W  m  :)  and  very 
weak  shear  mixing  (~  1  W  m  “).  However,  nutrient  data  collected  subsequently  as  part  of 
the  SO  GLOBEC  program  reinforced  the  idea  that  there  must  be  some  upward  material 
transport  from  the  UCDW  to  the  Antarctic  Surface  Water  (AASW)  ( Y.  Serebrennikova, 
pers.  comm.). 

The  analysis  in  Chapter  5  suggests  that  the  vertical  mixing  during  times  of  ice 
coverage  in  Marguerite  Bay  is  driven  by  episodic  wind  events.  Thus,  the  shipboard 
observations  used  in  the  analysis  in  Chapter  2  may  have  missed  significant  mixing 
events.  In  fact,  the  simulations  in  Chapter  5  showed  weak  fluxes  associated  with  near- 
inertial  oscillations  after  the  wind  event,  consistent  with  Chapter  2. 

Despite  the  fact  that  their  five-year  PWP  simulation  did  not  reproduce  the 
observed  stratification,  Meredith  et  al.  (2004)  found  ice  production  to  be  the  main  control 
on  upper  ocean  stratification.  That  work,  which  focuses  on  interannual  variability  in  the 
northern  arm  of  Marguerite  Bay,  may  have  missed  key  one-dimensional  vertical 
processes,  but  instead  attributed  the  model-data  discrepancy  to  lateral  processes.  The 
Rothera  Time  Series  program  measured  stratification  approximately  weekly,  but  did  not 
make  velocity  measurements,  vital  to  any  estimate  of  shear  mixing.  However,  careful 
inspection  of  the  stratification  before  and  after  strong  wind  events  during  times  of  lower 
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ice  concentrations  (more  freely  drifting  ice)  could  potentially  give  further  insight  into  the 
vertical  mixing  mechanisms  at  work  in  the  region. 

This  thesis  suggests  that  the  important  factors  to  consider  in  vertical  mixing 
beneath  sea  ice  are  not  only  the  frequency  and  magnitude  of  strong  wind  events,  but  the 
ice  conditions  at  the  time  of  the  events.  In  particular,  mixing  should  occur  under  ice 
conditions  favorable  for  the  transfer  of  momentum  from  the  atmosphere  through  the  ice 
to  the  upper  ocean,  i.e.  less  of  the  wind  stress  balanced  by  internal  ice  stresses.  Lower 
ice  concentrations  (<~80%)  (Chapter  5)  and  larger  distances  from  coastal  boundaries 
(-400  km)  (Thorndike  and  Colony,  1982)  are  generally  conditions  which  favor  a  sea  ice 

momentum  balance  closer  to  the  free  drift  solution  (i.e.,  V  •  cr  =  0 ,  see  Chapter  4).  In 
fact,  the  first  set  of  simulations  in  Chapter  5  showed  that  the  air-sea  momentum  flux 
through  free  drift  sea  ice  is  nearly  identical  to  that  of  an  ice-free  ocean  (not  yet  taking 
into  account  the  importance  of  sea  ice  effects  on  air-sea  heat  flux).  Higher  ice 
concentrations  and  associated  internal  ice  stresses  not  only  insulated  the  wAP  shelf  upper 
ocean  from  air-sea  heat  flux,  but  also  from  momentum  flux,  which  was  shown  in  Chapter 
5  to  be  the  primary  factor  in  mixing. 

6.3.  Implications  and  Future  Work 

Vertical  fluxes  of  momentum,  heat  and  salt  beneath  an  ice-covered  ocean  remains 
an  open  topic  in  physical  oceanography.  Researchers  now  recognize  that  in  order  to 
understand  and  simulate  sea  ice,  the  surface  ocean  must  be  taken  into  account  (e.g.  Heil 
and  Hibler,  2002).  Furthermore,  observations  continue  to  highlight  the  dominant  role  of 
internal  ice  stresses  in  balancing  the  wind  stress  (e.g.  Chapter  4;  Uotila  et  al,  2000).  Ice 
dynamics  is  fundamentally  a  spatially  two-dimensional  problem,  and  therefore  very 
difficult  to  measure. 

Exciting  research  is  currently  underway  to  observe  the  mixing  processes 
discussed  in  Chapter  5.  A  field  program  known  as  1CORTAS  (Ice  Covered  Ocean 
Response  to  Atmospheric  Systems)  is  being  organized  to  observe  and  quantify  mixing 
associated  with  high-frequency  ice  motion  and  formation  in  the  Laptev  Sea  in  the  Arctic 
Ocean  by  deploying  an  array  of  moorings  similar  to  those  used  in  the  SO  GLOBEC  field 
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program  (H.  Simmons,  pers.  comm.).  This  program  has  initial  plans  to  use  the  methods 
described  in  Chapter  3. 

Moored  observations  of  temperature  and  salinity  in  the  upper  ocean  beneath  sea 
ice  (including  those  made  in  SO  GLOBEC)  often  encounter  problems  with  drifting  deep- 
keeled  icebergs  (i.e.  keel  depths  >  ~50  m).  There  have  been  significant  technological 
advances  to  address  this,  such  as  transmitting  upper  ocean  temperature  and  salinity  data 
to  a  deeper  logger  via  an  acoustic  modem.  The  upper  instrumentation  is  attached  to  the 
deeper  portion  of  the  mooring  via  breakable  links,  so  if  an  iceberg  moved  across  the 
mooring  site,  the  instruments  may  be  lost,  but  the  data  and  the  rest  of  the  mooring  may 
still  be  recovered  (C.  Lee,  pers.  comm.).  There  have  also  been  advances  in  winched 
systems  which  profile  when  no  iceberg  is  detected,  but  remain  deep  when  an  iceberg  is 
drifting  through  the  mooring  site  (F.  Straneo,  pers.  comm.).  Measurements  of  the  upper 
ocean  stratification  will  not  only  provide  the  measured  sound  speed  profile  useful  for  the 
methodology  described  in  Chapter  3,  but  will  provide  a  basis  for  comparison  with 
simulations  such  as  those  in  Chapter  5.  In  particular,  the  observed  response  of  the  upper 
ocean  to  a  storm  event  and  any  subsequent  high-frequency  ice  motion  and  brine 
production.  A  favorable  model-data  comparison  (or  lack  thereof)  could  point  towards  the 
importance  of  physics  not  included  in  the  model,  such  as  lateral  advcction. 

Finally,  a  full  three-dimensional  model  of  sea  ice  and  the  upper  layers  of  the 
ocean  in  an  idealized  marginal  sea  with  idealized  forcing  could  possibly  illuminate  some 
of  the  important  processes  in  mixing  beneath  sea  ice.  In  particular,  the  divergence  and 
convergence  of  high-frequency  ice  motion  would  provide  a  more  realistic  time  series  of 
the  spatially  inhomogeneous  open  water  fraction,  rather  than  the  idealized  open  water 
fraction  used  in  Chapter  5.  Modern  sea  ice  models,  such  as  that  of  Zhang  and  Hiblcr 
(1997)  are  capable  of  reproducing  realistic  lead  structures  (J.  Zhang,  pers.  comm.).  Run 
at  high  enough  spatial  resolution,  the  model  could  simulate  buoyancy  flux  at  distinct 
points  on  the  surface,  rather  than  an  average  flux  distributed  over  an  area  containing  leads 
and  ice  as  envisioned  in  Chapter  5.  This  could  potentially  produce  convective  plumes 
and  give  insight  into  any  important  effects  on  the  three-dimensional  structure  of  the 
stratification  and  mixing.  The  results  from  Chapter  5  provide  an  excellent  starting  point 
for  a  full  three-dimensional  modeling  study. 
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